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This  report  Is  divided  Into  two  separate  sections,  Part  I  and 
Part  II,  which  correspond  to  the  major  contractual  tasks. 

^  Part  I  covers  the  work  on  *  Laser  Induced  Surface  Process! ng*. 

Part  II  covers  the  work  on  "Laser  Blow-Off  Ion  Implantation". 

Each  part  has  separate  section  numbering,  references  and 
appendices. 
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PART  I  -  LASER  INDUCED  SURFACE  PROCESSING 


SUMMARY  OF  FINAL  REPORT,  PART  I 


'pfitT  I 

The  objective  of  Task  4  of  this  program  was  to  evaluate  and 
characterize  laser-induced  processing  of  semiconductor  materials  and 
devices. .The  work  described  here  was  performed  under  Contract  No. 

MOA  903-78^>0284  during  the  period  of  June  1978  to  July  1979.  The 
results  of  thfs  work  were  also  presented  In  two  semi-annual  reports,  and 
the  salient  partsNqnd  conclusions  are  presented  here.  In  particular, 
reprints  of  two  publl^ed  papers  from  this  work  are  Included  as 
appendices  A  and  B  (also  as  References  1  and  2).  Appendix  C  contains  an 
unpublished  derivation  of  the  contribution  of  free  electrons  to  optical 
absorption  In  semiconductors. 

The  TRW  contribution  to  laser-induced  processing  of  semi-conductor 

surfaces  consisted  of  the  demonstration  of  the  effect  of  lattice  damage 

on  the  laser- Induced  heating,1  and  the  development  of  annealing 

2 

techniques  for  lightly  damaged  materials.  This  new  technique,  pioneered 
at  TRW,  consists  of  combining  a  laser  anneal  with  a  short  low  temperature 
thermal  anneal.  It  exploits  the  unique  ability  of  laser  processing  to 
dissolve  lattice  defects  without  providing  sufficient  energy  for 
melting  or  regrowth.  A  detailed  discussion  of  the  effect  of  damage  on 
the  laser-induced  heating  Is  given  In  Section  2.1  and  Appendix  A. 

In  our  new  technique,  the  sample  Is  first  Irradiated  by  a  pulsed  laser. 
This  removes  lattice  defects.  This  defect-free  material  Is  then  subjected 
to  a  moderate  thermal  anneal.  This  results  In  full  electrical  activation 
without  significant  change  in  the  Impurity  distribution.  Thus,  by 
combining  laser  irradiation  with  short  low  temperature  thermal  anneal.  It 
Is  possible  to  both  eliminate  defects  and  electrically  activate  lightly 
damaged  semi-conductors.  Since  melting  Is  avoided,  the  dopant  distribution 
remains  as  Implanted.  The  results  of  these  experiments  are  summarized  In 
Section  2.2  and  Appendix  B.  Preliminary  experiments  have  also  shown  that 
this  technique  can  be  applied  to  high  fluence  Implants. 
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1.0  INTRODUCTION 


In  our  work  at  TRW,  we  have  Investigated  the  Important  difference 
between  annealing  of  heavily  and  lightly  damaged  materials.  Our  original 
results  confirmed  findings  recently  reported  In  the  literature:  for  lightly 
damaged  samples,  the  crystal  regrowth  was  poor,  and  only  partial  electrical 
activation  was  obtained,  though  the  displacement  damage  was  removed.  In 
contrast,  for  high  fluence  samples,  excellent  activation  was  obtained. 

The  mechanism  for  the  Q-swItched  laser  annealing  Is  believed  to 

Involve  melting  of  the  amorphlzed  surface  followed  by  liquid  phase 
3  4 

regrowth.  ’  Here,  the  defects  are  almost  completely  removed  and  a  full 
activation  of  dopants  at  concentrations  far  exceeding  the  equilibrium  limit 

5 

of  solid  solubility  has  been  reported.  This  method  has  Its  biggest 

potential  In  applications  where  slight  changes  In  the  Implanted  profiles 

can  be  tolerated.  This  Is  because  melting  results  In  some  Impurity 

redistribution.  For  ruby  annealed  silicon  samples,  typical  profile 

1  5 

widths  are  around  0.3  ym.  * 

Our  early  experiments  lead  us  to  thoroughly  Investigate  the  effect 
of  damage  on  the  laser-induced  heating.  For  the  theoretical  effort,  we 
have  developed  a  model1  which  explains  the  observed  differences  between 
annealing  of  heavily  and  lightly  damaged  surfaces.  This  Is  done  by  taking 
Into  account  diffusion  of  the  photoexclted  carriers.  This  Is  discussed 

In  Section  2.1  and  Appendix  A.  In  our  experimental  effort,  we  have 

2 

subsequently  developed  a  novel  technique  which  permits  laser  annealing 
of  lightly  damaged  materials.  This  has  the  very  attractive  feature  of 
maintaining  the  as-implanted  Impurity  profiles.  This  Is  discussed  In 
Section  2.2  and  Appendix  B.  Preliminary  experiments  have  also  shown  that 
this  technique  can  be  applied  to  high  fluence  Implants. 
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2.0  RESULTS 


2.1  The  Effect  of  Damage  on  Laser  Heating.  The  conversion  of 
optical  to  thermal  energy  during  laser  annealing  consists  of  two  major 
processes : 

•  Absorption  of  laser  light 

•  Conversion  of  absorbed  optical  energy  to  heat. 

In  strongly  absorbing  semiconductors,  the  dominant  absorption  mechanism 
at  frequencies  higher  than  the  bandgap  frequency  Is  Interband  transitions. 
The  absorbed  energy  Is  released  to  heat  through  the  thermal Izatlon  of  the 
photoexclted  carriers,  l.e.,  heat  Is  distributed  Into  a  depth  comparable 
to  carrier  diffusion  length.  In  good  quality  materials,  the  carrier 
diffusion  lengths  are  longer  than  the  optical  absorption  depth.  Thus,  the 
heating  volume  Is  enlarged,  and  for  the  same  laser  power,  the  average 
temperature  Is  lower  than  In  heavily  damaged  material. 

The  carriers  recombine  preferentially  at  defect  sites,  producing 
localized  thermal  spikes.  If  the  density  of  the  nonradlatlve  recombination 
centers  Is  high,  the  average  temperature  Is  high.  This  Is  the  case  In 
heavily  damaged  materials.  At  these  elevated  temperatures  epitaxial 
regrowth  (liquid  phase)  and  Impurity  activation  occur.  If  the  density 
Is  low,  the  average  temperature  Is  low.  This  Is  the  case  In  lightly  damaged 
materials.  Here  the  application  of  laser  fluxes  equivalent  to  that  used  In 
high  fluence  Implants  produces  no  melting  In  low  fluence  Implants.  No 
regrowth  and  little  or  no  dopant  activation  occur  even  though  the  defects 
dissolve  due  to  the  preferential  supply  of  energy  at  the  defect  sites. 

Thus,  It  is  possible  In  lightly  damaged  materials  to  remove  the 
displacement  damage  even  without  providing  sufficient  average  energy  to 
activate  Implanted  Impurities. 
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Our  detailed  thermal  analysis  supports  this  picture.  We  have 
developed  a  model^  which  explains  the  observed  differences  between 
annealing  of  lightly  and  heavily  damaged  materials  by  taking  Into  account 
the  diffusion  of  the  photoexcited  carriers.  Previous  thermal  models  were 
concerned  with  heavily  damaged  materials,  In  that  they  assumed  that  the 
absorbed  optical  energy  is  converted  to  heat  In  the  Immediate  vicinity 
of  the  absorbed  photons.  In  our  work,  the  time  evolution  of  temperature 
profiles  was  calculated  for  pulsed  ruby  and  Nd:YAG  illumination  of  silicon 
and  ruby  Illumination  of  GaAs.  The  material,  optical,  and  electronic 
properties  were  varied,  to  represent  varied  material  quality.  Surface 
recombination  of  the  carriers  was  also  taken  into  account.  The  effect 

of  the  damage  on  the  heating  was  found  to  be  much  more  pronounced  for 

Nd:YAG  than  for  ruby  irradiation  because  a  larger  portion  of  the  ruby  laser 

energy  (excess  above  bandgap)  Is  given  rapidly  to  the  lattice  through 

quasl-thermallzatlon  Inside  the  semiconductor  energy  bands.  Thermal 
spikes  were  formed  as  the  result  of  the  preferential  recombination  of  the 
carriers  at  the  surface.  For  a  sufficiently  fast  surface  recombination, 
the  temperature  profile  may  develop  a  secondary  maximum  in  the  bulk. 

2.2  Laser  Annealing  of  Low  Fluence  Ion- Imp! anted  Silicon.  The  most 

successful  applications  of  Q-swItched  lasers  have  been  for  annealing  of 

3 

high  fluence  Implants.  The  method  Involves  a  pulsed  laser  such  as  a 
Q-switched  Nd:YAG  or  ruby  laser,  with  a  typical  pulse  length  between  10 
and  200  nsec,  and  energy  density  between  .5  and  3  J/cm  .  The  annealing 

mechanism  is  believed  to  be  melting  of  the  amorphlzed  surface  followed  by 

4  5 

liquid  phase  regrowth.  *  Here  there  Is  almost  complete  removal  of  defects 

and  a  full  activation  of  dopant  concentrations  far  exceeding  the  limit  of 
3  5 

solid  solubility.  '  However,  this  method  can  be  only  applied  where  slight 

changes  In  Implanted  profile  can  be  tolerated.  This  Is  because  melting 

results  In  some  Impurity  redistribution.  Typical  profile  widths  In  ruby 

3  6 

laser  annealed  silicon  are  around  0.3  ym.  * 
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The  second  laser  anneal  method  employs  a  cw  laser  operated  in  a 

scanning  mode.  For  example,  a  7W  cw  argon  laser  beam  scanned  at 

2.76  cm/sec  and  producing  lines  of  22  pm  width  was  used  to  anneal  to  a 
o  7 

1000  A  depth  in  as-implanted  silicon.  The  recrystallization  is  through 
solid  state  regrowth,  because  of  the  low  powers  and  comparatively  long 
anneal  times.  The  dopant  profiles  are  not  changed  by  cw  laser  annealing. 

The  as-implanted  impurity  profiles  can  be  maintained  even  in  pulsed 
laser  annealing  by  using  techniques  developed  at  TRW.  These  techniques 
consist  of  combined  laser  irradiation  and  low  temperature  thermal  anneal. 

In  our  work  at  TRW  we  have  systematically  investigated  the  differences 

between  annealing  of  lightly  damaged  materials  and  heavily  damaged  materials. 

2 

We  have  applied  ruby  laser  pulses  of  1.5  to  2.5  J  cm  to  silicon  surfaces 

oriented  in  [100]  and  [111]  directions  and  implanted  with  lO1^  to  1015  B  cm'2 
14  14  -2 

and  10  to  5  x  10  B  cm  .  Boron  implantations  were  at  50  keV, 

phosphorous  at  80  keV.  It  was  found  that  laser  energy  densities  used  for 

high  fluence  implants  did  not  produce  melting  or  regrowth  in  the  low 

fluence  implants.  Although  the  electrical  activation  was  poor,  the  lattice 

defects  (e.g.,  dislocation  loops)  were  removed.  We  have  concluded  that  the 

average  temperature  induced  in  the  lightly  damaged  material  was  substantially 

lower  than  the  temperature  induced  in  the  heavily  damaged  material.  This 

permitted  only  a  small  fraction  of  the  Implanted  ions  to  migrate  to 

substitutional  positions.  The  damage  was  removed  through  the 

preferential  heating  at  the  defect  sites.1  However,  when  a  sequential 

moderate  thermal  anneal,  10  min  at  1000°C,  was  applied  to  these  samples, 

a  full  electrical  activation,  and  without  significant  impurity  redistribution, 

was  obtained.  Thus,  by  combining  a  laser  anneal  with  a  short  thermal 

anneal,  it  was  possible  to  achieve  defect  elimination  as  well  as  complete 

2 

electrical  activation.  This  was  the  first  reported  successful  use  of 
pulsed  laser  radiation  for  annealing  of  low  fluence  implants. 
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Preliminary  experiments  have  shown  that  we  can  also  apply  these 

techniques  to  high  fluence  Implants.  One  type  of  Implantation,  that  of 

14  2 

5  x  10  *  P  cm"*  Into  a  [111]  surface,  resulted  in  amorphism.  Our  early 
experiments  as  well  as  others  in  the  literature  showed  that  Irradiation 
of  amorphized  surfaces  by  a  Q-swItched  ruby  laser  results  in  melting. 

Thus,  In  order  to  avoid  melting,  the  damaged  region  should  first  be 
recrystallized.  This  was  done  by  introducing  a  moderate  thermal  anneal  before 
the  Irradiation.  To  date,  the  best  results  were  obtained  with  a  two-step 
thermal  anneal,  30  min  at  550°C,  and  10  min  at  1000°C.  This  caused  the 
amorphous  layer  to  regrow  into  single  crystallinity  by  solid  phase  epitaxy. 
The  subsequent  irradiation  eliminated  defects  without  melting.  Not 
only  did  we  eliminate  the  defects,  but  by  avoiding  melting,  we  also 
maintained  the  desirable  as-implanted  impurity  profiles. 

The  demonstrated  techniques  have  an  enormous  potential  for 
practical  applications.  Low  fluence  implants  account  for  more  than 
50  percent  of  the  potential  implant  uses  since  all  bases  for  junction 
transistors  fall  Into  this  category.  The  ability  to  anneal  high  fluence 
Implants  without  dopant  redistribution  is  important  for  the  production 
of  shallow  emitters,  which  are  essential  for  an  entire  category  of 
devices,  such  as  microwave  transistors.  The  ability  to  anneal  material 
without  changing  the  Implanted  profiles  would  permit  full  exploitation 
of  the  Ion-Implant  unique  capability  of  creating  narrow  accurately 
controlled  impurity  profiles. 
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The  effect  of  carrier  diffusion  on  laser  heating  of  lightly  damaged 
semiconductors 

Jaroslava  2.  Wilcox 

TR  W  Defense  and  Space  Systems  Group.  Redondo  Beach.  California  90278 

(Received  10  September  1979;  accepted  for  publication  6  December  1979) 

It  has  been  observed  that  the  surface  temperature  obtained  during  laser  heating  of  lightly 
damaged  material  is  generally  considerably  lower  than  the  temperature  obtained  in  heavily 
damaged  material.  We  develop  a  model  which  explains  this  observation  by  taking  into  acount  the 
diffusion  of  the  photoexcited  carriers.  The  time  evolution  of  temperature  profiles  was  calculated 
for  pulsed  ruby  and  Nd  :  YAG  illumination  of  silicon  and  ruby  illumination  of  GaAs.  The 
material,  optical,  and  electronic  properties  were  varied,  according  to  the  material  quality.  Surface 
recombination  of  the  carriers  was  also  taken  into  account.  For  the  same  laser  power,  the  induced 
temperature  was  lower  and  the  heating  times  were  longer  in  good-quality  materials  than  in 
heavily  damaged  materials.  This  is  because  of  longer  carrier  difusion  lengths  in  good-quality 
material.  The  effect  of  the  damage  on  the  heating  was  found  to  be  much  more  pronounced  for 
Nd  :  YAG  than  for  ruby  irradiation  because  a  larger  portion  of  the  ruby  laser  energy  (excess 
above  band  gap)  is  given  rapidly  to  the  lattice  through  quasithermalization  inside  the 
semiconductor  energy  bands.  Thermal  spikes  are  formed  as  the  result  of  the  preferential 
recombination  of  the  carriers  at  the  surface.  For  a  sufficiently  fast  surface  recombination,  the 
temperature  profile  may  develop  a  secondary  maximum  in  the  material  interior. 

PACS  numbers:  79.20.Ds,  66.70.  +  f,  72.20. Jv,  8l  .40.Ef 


I.  INTRODUCTION 

There  is  considerable  interest  in  the  use  of  lasers  for 
processing  of  semiconductor  surfaces.1  The  majority  of  ex¬ 
periments  have  been  in  ^-switched  ruby  and  Nd  :  YAG  la¬ 
ser  annealing  of  comparatively  high-dose  implantations.  In 
these  experiments,  the  resulting  activated  dopant  concentra¬ 
tions  were  often  in  excess  of  the  limits  of  equilibrium  solid 
solubility,  and  the  shapes  of  the  ion-implanted  profiles  were 
changed  as  a  result  of  the  illumination.'  This  led  to  the  sug¬ 
gestion  that  the  near-surface  regions  of  the  irradiated  mate¬ 
rial  actually  melted  and  underwent  a  type  of  liquid  epitaxy  in 
these  experiments.  The  suggested  model  is  strongly  support¬ 
ed  by  numerical  results  of  several  model  calculations.3 

Heavily  damaged  material  has  been  successfully  an¬ 
nealed  also  with  scanned  cw  Kr  and  Ar  ion  lasers*  No  sig¬ 
nificant  dopant  redistribution  was  observed  after  the  anneal¬ 
ing,  indicating  that  the  regrowth  was  through  solid-state 
epitaxy. 

Laser  annealing  of  lightly  damaged  material  has  en¬ 
joyed  much  less  attention.  The  available  data  indicate  a 
strong  sensitivity  to  the  dose  of  the  implanted  impurities. 
Excellent  single-crystal  regrowth  was  obtained  in 
Nd  :  YAG  annealing  of  high-dose  samples,  while  in  lower- 
dose  samples  the  regrowth  was  poor.’  Furthermore,  laser 
powers  which  melted  implanted  silicon  did  not  affect  the 
unimplanted  regions  of  the  material.6  For  ruby  laser  anneal¬ 
ing,  carrier  densities  for  lower-dose  implants  were  essential¬ 
ly  the  same  as  those  obtained  by  the  thermal  annealing,  while 
in  the  higher-dose  range  the  carrier  density  was  much  high¬ 
er.3  The  use  of  0-switched  ruby  laser  for  annealing  of  low - 
fluence  silicon  was  studied  in  greater  detail  in  Ref.  8.  The 
fluences  studied  ranged  from  10M  to  lO”  B/cnr  and  10M  to 
5  x  10u  P/cm’.  The  laser  energies  were  1-2.5  J  cm  :.  No 
melting  and  only  partial  electrical  activation  occurred  as  a 


result  of  the  irradiation.  However,  the  implant  defects  were 
fully  dissolved."  Full  electrical  activation  was  obtained  (at  2 
J  cm1)  only  for  the  1 1 1 1 1  surface  amorphized  by  an  implant 
of  5  x  I0N  P/cm2.  For  nonamorphized  surfaces,  powers  of 
2.5  J  cm  -’  produced  only  partial  activation;  these  surfaces 
could  be  activated  fully  by  applying  a  sequential  thermal 
anneal.  In  a  related  series  of  experiments,  we  have  applied 
ruby  radiation  also  to  silicon  which  had  been  subjected  to 
neutron  bombardment."  Again,  no  electrical  activation  oc¬ 
curred  unless  the  wafers  were  additionally  damaged  by  Nc 
implant.  However,  the  surfaces  showed  melt  spots  in  the 
vicinity  of  defects,  indicating  that  the  laser  energy  was  re¬ 
leased  primarily  at  the  defects. 

The  present  work  addresses  the  difference  between  la¬ 
ser-induced  heating  of  heavily  and  lightly  damaged  materi¬ 
als.  It  presents  the  results  of  the  thermal  calculations  which 
take  into  account  the  diffusion  of  the  photoexcited  carriers. 
Previous  thermal  models  have  been  concerned  with  heavily 
damaged  materials,  in  that  they  assumed  that  the  absorbed 
optical  energy  is  converted  to  heat  in  the  immediate  vicinity 
of  the  absorbed  photons.  However,  in  strongly  absorbing 
semiconductors,  the  primary  effect  of  the  absorbed  radiation 
is  to  generate  carriers.  Heat  is  generated  when  the  carriers 
recombine,  i.e.,  after  they  have  diffused  a  distance  on  the 
order  of  carier  diffusion  length.  In  good-quality  materials, 
the  diffusion  lengths  may  be  longer  than  the  optical  absorp¬ 
tion  depth.  Hence,  the  carrier  diffusion  enlarges  the  heating 
volume,  and,  for  the  same  laser  power,  causes  the  surface 
temperature  in  lightly  damaged  material  to  be  lower  than  in 
heavily  damaged  material.  Furthermore,  a  large  number  of 
carriers  will  reach,  and  recombine  at,  the  surface,  making  it 
qualitatively  feasible  to  melt  surface  defects  in  lightly  dam¬ 
aged  materials. 

The  conversion  of  optical  energy  to  heat  consists  of  two 
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major  step*:  (i)  absorption  of  laser  light,  and  (ii)  conversion 
of  absorbed  optical  energy  to  heat.  We  consider  strongly  se¬ 
miconductors.  In  strongly  absorbing  semiconductors,  the 
dominant10  absorption  mechanism  at  frequencies  higher 
than  the  band-gap  frequency  is  interband  transitions.  The 
conversion  of  absorbed  energy  to  heat  is  accomplished 
through  thermalization  of  photoexcited  carriers.  Once  the 
conversion  has  been  accomplished,  the  temperature  distri¬ 
bution  is  determined  by  solving  the  thermal  diffusion  equa¬ 
tion  in  the  semiconductor.  At  issue  is  the  proper  representa¬ 
tion  of  the  source  term  in  the  heat  equation. 

In  the  model  used  in  this  work,  we  categorize  the  ther¬ 
malization  as  consisting  essentially  of  two  steps,  operative 
either  in  parallel  or  in  sequence:  relaxation  of  the  electron 
"excess"  energy  (above  the  band  gap),  and  relaxation  of  the 
band-gap  portion  of  the  energy.  The  categorization  is  conve¬ 
nient  because  the  electron  redistribution  inside  the  conduc¬ 
tion  band  proceeds  generally  at  a  different  rate  than  the  elec¬ 
tron-hole  recombination  across  the  semiconductor  band 
gap.  The  band-gap  portion  of  the  electron  energy  is  released 
to  heat  at  the  rate  ~£,  Jar where  r  is  the  carrier  lifetime, 
£„  is  the  semiconductor  band  gap,  and  dn  =  An{x,t )  is  con¬ 
centration  of  the  photoexcited  carriers,  dn  is  calculated  by 
solving  the  electron  diffusion  equation.  Thus,  the  band-gap 
portion  of  the  electron  energy  is  released  into  a  volume 
bounded  by  a  depth  equal  to  max(L,  a' '),  where  L  ~(0r)'"\ 
D  is  the  carrier  diffusivity,  and  a  is  the  optical  absorption 
coefficient. 

The  conversion  of  the  excess  portion  of  the  electron 
energy  is  assumed  in  this  work  to  be  instantaneous.  This  is 
because  the  characteristic  “quasithermalization”  frequency 
is  on  the  order  of  the  lattice  vibrational  frequency./.  as  10' :- 
10°  sec  ',  which  is  much  faster  than  any  other  rate  in  the 
problem." 

The  presence  of  defects  and  interfaces  also  affects  the 
heating  profiles.  Interfaces  affect  heating  mainly  because  of 
surface  recombination  of  the  photoexcited  carriers.  The  ef¬ 
fect  of  the  surface  will  be  taken  into  account  in  this  work 
through  the  surface  recombination  velocity  u, . 

Bulk  defects  may  be  categorized  as  localized  or  distrib¬ 
uted.  Both  affect  optical  absorption  as  well  as  carrier  diffu¬ 
sion.  The  influence  of  the  distributed  defects  will  be  taken 
into  account  in  that  the  values  of  the  optical  absorption  coef¬ 
ficient,  carrier  lifetime,  and  diffusivity  will  be  different  for 
heavily  damaged  than  for  good-quality  materials.  However, 
the  crystal  will  be  considered  homogeneous,  and  its  proper¬ 
ties  independent  of  time.  In  this  sense,  the  values  of  the  mate¬ 
rial  parameters  used  are  some  “average”  values  of  the  pa¬ 
rameters  during  the  processing.  The  main  reason  for  this  is 
mathematical  simplicity.  The  localized  bulk  defects  will  not 
be  treated  here  explicitly.  However,  their  effect  on  the  heat¬ 
ing  should  be  similar  to  that  of  the  surface  defects.  This  is 
because  both  surface  and  localized  defects  are  fast  recombi¬ 
nation  centers,  with  the  subsequent  release  of  the  carrier  en¬ 
ergy.  This  leads  to  the  formation  of  thermal  spikes  at,  and  a 
possible  dissolution  of,  localized  defects. 

As  has  already  been  mentioned,  the  “band-gap  energy 
source  term"  in  the  heat  equation  requires  knowledge  of  the 
electron  concentration,  n(x.i ).  For  constant  r,  D,  and  v. ,  we 


have  solved  the  electron  diffusion  equation  exactly  for  a 
square-shaped  laser  pulse  (no  drift  current ).  The  solution  for 
n(x,f )  is  a  closed-form  expression  consisting  of  complemen¬ 
tary  error  functions.  To  the  best  of  our  knowledge,  this  solu¬ 
tion  has  never  before  been  reported  in  the  literature.  The 
existence  of  a  closed-form  expression  for  n(x,r )  is  important 
for  our  purposes  since  it  avoids  numerical  integration  of  the 
electron  diffusion  equation.  The  heat  transport  equation  is 
then  solved  numerically,  with  the  temperature  dependence 
of  the  thermal  conductivity  included. 

In  Sec.  II,  we  describe  the  mathematical  basis  of  the 
model.  In  Sec.  Ill,  we  discuss  the  results  of  the  numerical 
computations  of  the  temperature.  Section  IV  summarizes 
the  results. 


It.  MODEL 

The  photon-to-thermal  conversion  consists  of  photoex¬ 
citation  and  thermalization  of  the  photoexcited  carriers,  and 
heat  transport. 


A.  Concentration  ot  photogonoratod  carriers 


The  instantaneous  concentration  of  the  photogenerated 
carriers  is,  in  the  one-dimensional  approximation  and  in  the 
absence  of  electric  fields,  obtained  by  solving  the  electron 
transport  equation. 


r  dx- 


M) 


where  the  generation  function  g  gives  the  photogeneranon 
rate  of  the  carriers. 


g  =  aN exp(  -  ax) ,  (2) 

where  a  is  the  optical  absorption  coefficient  and  .V  =  .V  (r )  is 
the  photon  flux  incident  at  the  surface.  N  (/ )  =  F(i  )/hv. 
where  F(t )  is  the  absorbed  laser  power  and  hv  is  the  photon 
energy. 

dn&znfx,  i)  -  na,  (.') 

where  n0  is  the  equilibrium  concentration  and  r  and  D  arc 
minority-carrier  lifetime  and  diffusion  constant, 
respectively. 

In  Eq.  ( I ),  a,  r,  and/)  are,  in  general,  functions  of  both  x 
and  t.  The  x  dependence  arises  due  to  the  spatial  inhomo¬ 
geneities  in  the  semiconductor;  the  time  dependence  arises 
from  the  dependence  of  a,  D,  and  r  on  temperature.  Addi¬ 
tional  time  dependence  arises  due  to  changes  in  the  semicon¬ 
ductor  properties  as  the  damage  disappears  during  the  an¬ 
nealing.  For  high  carrier  concentrations,  a  may  be  a 
function  of  the  photocarrier  concentration  n(x.r ).  It  would 
be  quite  complicated  to  fully  account  for  the  effect  of  the 
changing  material  properties.  We  approximate  material 
properties  by  some  "average"  values  and  the  crystal  is  as¬ 
sumed  homogeneous,  with  the  exception  of  the  front  surface 
Taking  into  account  surface  recombination  of  the  earn¬ 
ers.  the  boundary  conditions  are 


Xx  =  0.O  = 


i  (wtoi  «„>  (4> 
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and 

n{x  « / )  *  n0  • 

where  v ,  is  the  surface  recombination  velocity  and  /  is  the 
semiconductor  thickness.  Equations  (4)  and  (5)  correspond 
to  a  finite  flux  of  carriers  at  the  semiconductor  front  surface, 
and  a  constant  carrier  concentration  at  the  slab  back  bound¬ 
ary.  The  initial  condition  is 

n(x,  t  =  0)  *  n„ .  W 

In  general,  if  AF.  a,  r.  and  0  are  functions  of  x  and/,  Eg.  | 


(I)  has  to  be  solved  numerically.  However.  Eq  ( I  >  can  be 
solved  exactly  analytically  for  the  special  case  when  a.  r.  and 
Pare  constants,  and  the  generation  function  A'  (i )  isa  square- 
shaped  laser  pulse. 

AF(/)=Ar„.  0  </</, 

=  0.  otherwise.  (7> 

In  Eg.  (7),  t,  is  the  length  of  the  laser  pulse  and  .V„  is  the 
number  (constant)  of  arriving  photons  at  the  semiconductor 
front  surface.  The  mathematical  methods  are  discussed  in 
the  Appendix.  The  result  is 


d»(x.  i )  = 


g.V0r  I  y  1  +aD/v, 
2(1  -f.  V)  L  ,  ±  I  -L/o.r 
(±  1  +  00/1',) 


exp(  ±  x/L )  erf]  ± 


m 


(aD/v 


(40/  )'  2 

-y-  exp  {^±ax  ~  —  (1  -  o‘£  :)j  erfc  (  +  x/(40f ')'•’  +  o(0r ')'  ’) 

- exp  (  —  !:  +  (/•„, +x)Me4.(g%  — 

-1X1  -v]t/D)  \  t  0/1  \0/  (40/)'  JJ, 


(I  -  00/1' 


where  /„  =  max(0,/ - /,),  the  diffusion  length  L  =  (Dr)'  \  I 
and  erf(z)  and  erfc(z)  are  error  and  complementary  error 
functions,  respectively.  Equation  (8)  is  evaluated  between 
the  indicated  upper  and  lower  boundaries.  /  and  /„.  respec¬ 
tively.  n(x.i )  appears  as  a  source  term  in  the  thermal  diffu¬ 
sion  equation  discussed  in  Sec.  II  B. 

B.  Thermal  transport 

The  temperature  distribution  in  the  semiconductor 
after  the  application  of  the  laser  pulse  is  obtained  by  solving 
the  thermal  equation. 

dT _ 3  f*  ar\  _  S(x,/)  (9) 

3t  3x  \c,  3x)  c, 

where  c,  and  «r  are  the  specific  heat  and  thermal  conductiv¬ 
ity.  respectively,  and  where  the  source  function  5  is  the  heat 
flux  deposited  in  the  semiconductor  lattice  per  unit  length. 
Equation  (9)  is  solved  subject  to  appropriate  boundary  con¬ 
ditions.  In  many  instances,  these  are  (i)  initial  temperature, 
(ii)  thermal  flux  at  the  front  surface,  and  (iii)  thermal  flux  or 
temperature  of  the  heat  sink  at  the  back  of  the  sample. 

1.  Ughtty  damaged  matonata 

Following  our  previous  discussion,  we  express  S  as  a 
summation  of  three  terms.  S„  S3,  and  S}.  S,  is  the  "excess" 
portion  of  the  electron  energy,  S2  is  the  band-gap  portion  of 
the  electron  energy,  and  S,  is  energy  released  by  carriers 
recombined  at  the  semiconductor  surface. 

The  excess  portion  of  the  electron  energy  is  released 
during  quasithennaliation  of  photoexcited  electrons  inside 
the  conduction  band.  Since  the  characteristic  frequency, 
lO'MO"  sec  ',  is  much  faster"  than  any  other  rate  in  the 
problem,  we  approximate 

S,(x,  / )  =>  (Av  -  Et)aN (/ ) exp(  —  ax) ,  (10) 

where  vis  the  laser  frequency  and  £,  is  the  band-gap  energy. 

The  energy-gap  portion  of  the  electron  energy  is  trans¬ 
ferred  to  heat  at  rate  r„  ' .  where  r„  is  the  nonradiative  life- 
time  of  the  electron, 


(8) 


S,(x,  t  )—E,  n(X,  l)  —  ■  (ID 

The  surface  recombination  contributes  a  term 

S,(x, /)  =  £„i\(»t(x. /)-n0)A(jt).  (!2) 

where  6  (x)  is  the  delta  function.  Mathematically,  the  effect 
of  T,(x./ )  will  be  to  introduce  a  finite  derivative  3T / dx  at 
x  =  0;  see  Eq.  (14)  in  what  follows. 

Similarly  as  for  the  electron  diffusion,  we  consider  an 
essentially  uniform  crystal.  It  is  worth  noticing,  however, 
that  the  effect  of  the  presence  of  localized  carrier  "sinks"  in 
the  crystal  on  temperature  could  be  included  through  source 
terms  of  the  type  of  Eq.  (12)  (which  in  turn  introduce  discon¬ 
tinuities  in  ST  /dx).  Physically,  this  corresponds  to  a  local¬ 
ized  deposition  of  the  energy,  and  a  subsequent  formation  of 
thermal  spikes  at  defect  sites.  This  makes  dissolution  of  lo¬ 
calized  defects  feasible.  In  this  sense,  location  x  =  0  of  this 
paper  may  be  viewed  as  a  defect  site,  with  the  recombination 
rate  for  the  carriers  equal  to  u,dn(x,r  )6  (x). 

Neglecting  loss  of  heat  radiated  from  the  front  surface 
(it  has  been  estimated  to  be  typically  six  orders  of  magnitude 
smaller  than  the  incoming  laser  flux),  Eq.  (9)  is  then  solved 
subject  to  the  following  initial  and  boundary  conditions: 

r(x,f*0)*ro,  (13) 

aT(X~  0-'l)  =  -  -8-  i\(n(0.  / )  —  n„) ,  (14) 

dx  K 

r(x  =  /,f)-r„,  (i5) 

where  T0  is  the  temperature  of  the  heat  sink  at  the  back  of  the 
sample.  Equation(l4)  was  obtained  by  substituting  Eq.  (12) 
into  Eq.  (9)  and  integrating  the  resulting  equation  across  the 
front  boundary.  Equation  (9)  has  to  be  solved  numerically 

2.  Heavily  damaged  materials  and  k  andcp  an  constants 

For  heavily  damaged  materials  and  «r,  cf  *  const.  Eq 
(9)  has  a  closed-form  solution.  This  is  because  for  extremely 
short  carrier  lifetimes,  n/,o  <  1  and  rDor  <  I  (i.e..  bulk  re¬ 
combination  is  faster  than  either  the  surface  recombination 
or  carrier  diffusion),  the  source  functions  simplify  to 
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TABLE  I  The  values  of  the  esaienal  aad  laser  parameters  used  m  the  calruiati 


' 

. 

11 

a  (cm  ') 

Damage 

r(s*c> 

O (cm  i« 

v  <cm/sev> 

cr  (J/cm‘> 

Rub> 

Nd  VAG 

1  IfiH'sl 

OaAs 

Good  quality 
to 

amoephutd 

Iff*- 10" 

200-13 

0-10' 

1  7S 

2  s  10* 

N  A 

10 

S« 

Good  quality 
to 

amorpfiucd 

Iff*- 10  " 

30-3 

0-10’ 

16 

2  V  l0‘-2  v 

10*  100-10’ 

too 

S:(x,  t  )mE,  ~£,g  =  E,aF(t )  ex p(  -  ax) 

(16) 


and 


S, (*./)= 0.  (17) 

This  yields  the  source  term 

SmS,  +  S2  — hvaN(i)exp(  — ax)  (18) 

and  the  boundary  condition 


(19) 


<?r(jt  =  0,/) 

3x 

The  solution  to  Eq.  (9)  is  then  obtained  by  taking  the 
limit  l/r-»0  in  Eq.  (A  10)  of  the  Appendix.  This  yields 

Fn 


jr(x.»)  = 


Tan 

Xerfc 

+  4a 


exp(  ±  ax) 


Hrv),.£; 

a(~)  *  i 

+  2«x  erf  ,  (20) 

where  F^mhvNo,  and  the  upper  and  lower  boundaries  are 
/ '  =  t  and  i0  *  max(0.r  -  /, ),  respectively. 


IH.  NUMERICAL  RESULTS 

Equation  (9)  was  solved  numerically  for  ruby  and 
Nd :  YAO  lasers.  GaAs,  and  silicon.  The  material  param¬ 
eters  used  ranged  from  those  appropriate  for  heavily  dam¬ 
aged  materials  to  those  corresponding  to  high-quality  mate¬ 
rials.  The  values  of  these  parameters  are  summarized  in 
Table  I.  The  thermal  conductivity  sr(  T)  was  represented  by  a 
summation  of  exponentials  chosen  to  At  data  measured  for 
pure  OaAs  and  silicon. 11  The  carrier  lifetimes  were  taken  as 
r  -  v  The  laser  pulse  lengths  used  (10  nsec  for  ruby  and 
100  nsec  for  Nd :  YAO  laser)  are  the  optimum  pulse  lengths 
for  heavily  damaged  materials. 

Samples  of  the  numerical  results  are  shown  in  Figs.1- 
12.  The  profiles  are  shown  as  a  function  of  the  normalized 
distance  ax  away  from  the  semiconductor  front  surface, 
with  time  (in  fractions  and  multiples  of  the  pulse  length  /, )  as 
a  parameter.  The  symbols  shown  (in  cp  units  unless  speci¬ 
fied  otherwise)  have  their  usual  meaning  The  normalized 
temperature  ( T  -  T,MF,/*,aY'  »  dimension  leu,  «r„  is  a 


quantity  (equal  to  unity)  which  has  the  dimension  of  thermal 
conductivity.  F0  is  the  actual  absorbed  power  (in  W/cnr). 
The  profiles  shown  correspond  to  4  rw,  ~900- 1 300  ‘C 
Oust  below  the  melting),  the  exact  number  can  be  calculated 
from  parameters  specified  in  each  figure.  Before  inspecting 
the  profiles,  it  is  convenient  to  qualitatively  summarize  the 
main  results  of  these  calculations.  This  will  aid  understand¬ 
ing  when  inspecting  the  profiles. 


For  heavily  damaged  materials,  the  effects  of  earner 
diffusion  and  surface  recombination  on  heat  generation  can 
be  neglected.  The  simulated  profiles  approached  those  cal¬ 
culated  using  Eq.  (20)  [Eq.  (20)  assumes  x  m  const,  n\D 
<  I,  and  rDa*  <  I).  As  compared  to  Eq.  (20),  the  effect  of 
k  =  ic{T)on  the  profiles  was  to  slightly  increase  the  surface 
temperature  and  narrow  the  profile  width.  This  is  because 
•KT)  decreases  with  increasing  temperature. 

1.  Surfaca  tampan  tun  and  prom  width 

For  the  pulse  length  shorter  than  the  thermal  diffusion 
time,  /,  <  t,h,  where  t,h  =  cp/xa~.  the  profiles  are  optical  ab¬ 
sorption  depth  limited.  The  profile  width  ~o  ',  and  JT„1JV 
grows  linearly  with  both  the  pulse  length  and  optical  absorp¬ 
tion  coefficient,  4  rm<>  ~F„tra/cr .  For  the  pulse  length 
longer  than  rfk .  tr  >  t,h ,  the  profile  is  thermal  diffusion  limit¬ 
ed.  Its  width  is  asbe t,/cf)' and  4  Tm„  is  independent  of  a 
and  proportional  to  the  square  root  of  the  pulse  length, 

^FJt,/Kcpy-'  (the  surface  temperature  approaches 
saturation  for  extremely  long  pulses).  In  either  case.  4  TBU, 
is  approximately  linear  with  pulse  peak  power  (4  rn,„  is  ex¬ 
actly  linear  with  F0  for  k  *  const).  The  heating  is  most  effi¬ 
cient  (least  laser  energy  required)  when  tr  ~ttk .  The  optical 
depth-limited  profiles  are  generally  narrower  than  the  ther¬ 
mal  diffbsion-limited  profiles. 

2.  Tampantura  dacay 

After  the  end  of  the  laser  pulse,  the  temperature  drops 
and  the  profile  broadens  as  the  heat  propagates  into  the  bulk. 
The  decay  is  much  faster  for  strong  absorption  than  for  weak 
absorption.  This  is  because  the  thermal  diffusion  time 
c,  /«*J  is  much  longer  for  small  a  than  for  large  a.  For  ex¬ 
ample, fora  ■  2  X 10' cm  ',  the  surface  temperature  typical- 
lydroppedtoone-halfofitsmaximum valueat/>  100/,  For 
a  *  2  x  I04  cm  ',  the  half-time  was  only  t^ftr 
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For  comparison  purposes,  Eq.  (9)  was,  for  heavily  dam¬ 
aged  materials,  also  solved  for  the  “Gaussian"  pulse  shapes. 
We  have  observed  really  no  significant  difference  between 
the  shapes  of  the  thermal  profiles  created  by  square  and 
Gaussian  pulse  shapes.  The  resulting  temperature  distribu¬ 
tions  were  affected  much  more  by  the  condition  a  '  versus 
(*ir/cr),,!  than  by  the  pulse  shape.  For  the  Gaussian  pulse, 
slightly  higher  laser  powers  were  required  to  achieve  the 
same  surface  temperature  as  for  the  square-shaped  pulse. 
The  Gaussian  profiles  tend  to  decay  more  slowly  than  the 
square-shaped  pulse  profiles. 

ff.  Lightly  damaged  matarWa 

The  main  effect  of  carrier  dilhision  on  temperature  pro¬ 
files  consists  of  the  following: 

(i)  Higher  laser  powers  are  required  to  raise  the  surface 
temperature  of  a  good-quality  material  than  of  a  damaged 
material. 

(ii)  Temperature  profiles  in  good-quality  materials  are 
generally  wider  than  in  damaged  materials. 

(iii)  The  time  (Oat  which  the  surface  temperature 
reaches  maximum,  and  time  (r,/: )  required  for  the  surface 
temperature  to  drop  to  one  half  of  its  maximum  value  after 
the  end  of  the  laser  pulse,  are  longer  in  good-quality  materi¬ 
als  than  in  damaged  materials. 

(iv)  The  effect  of  carrier  diffusion  on  temperature  pro¬ 
files  is  much  stronger  for  Nd  :  YAG  laser  than  for  ruby  laser 
irradiation. 

(v)  The  effect  of  the  “sinking”  of  the  carriers  is  to  pro¬ 
duce  thermal  spikes  at  the  surface  and  lower  temperature  in 
the  material  interior.  For  a  sufficiently  fast  surface  recombi¬ 
nation.  the  profiles  may  develop  a  secondary  local  maximum 
in  the  material  interior. 

Conclusions  (i)-(iv)  »re  the  direct  consequence  of  the 
model  used  for  the  transfer  of  the  photon  energy  to  heat.  The 
excess  and  band-gap  portions  of  the  electron  energy  heat 
volumes  bounded  by  the  depths 


<r«  =  max[a\  (x^/c,)1'1] 

(21) 

r.i). 

(22) 

respectively. 

Thus,  in  good-quality  materials  where  the  carier  diffu¬ 
sion  length  L>Otg,  the  band-gap  portion  of  the  energy  heats 
a  volume  which  is  larger  than  the  one  heated  by  the  excess 
energy.  This  results  in  temperatures  which  are  generally 
lower  (for  the  same  laser  power)  than  those  in  damaged  ma¬ 
terials.  Furthermore,  since  both  the  electron  lifetime  (r)  and 
the  thermal  diffusion  time  associated  with  the  band-gap  por¬ 
tion  of  the  energy 

*wa-Vira/ '«*•  <23> 

may  be  quite  long,  the  heating  times  0m  and  f,,2)  may  be 
substantially  longer  than  t, .  Clearly,  the  effect  of  carrier  dif¬ 
fusion  on  the  temperature  will  be  greater  for  Nd :  YAG  than 
for  ruby  laser  irradiation.  This  is  because  Nd :  YAG  fre¬ 
quency  is  very  close  to  the  silicon  band-gap  frequency;  the 
excess  portion  of  the  electron  energy  is  s»l%.  For  ruby  irra¬ 
diation  of  silicon,  the  excess  portion  of  the  electron  energy  is 
large,  os  35<3( ,  and  is  not  affected  by  electron  diffusion. 

asm  j  Apt*  PTiy* .  Vot  St.  No.  8, May  1900 


Conclusion  (v)  is  the  unique  consequence  of  the  surface 
recombination  of  the  carriers.  For  t>,  —  0,  the  tempera)  ure  is 
smooth  and  monotonically  decreasing  with  the  distance 
away  from  the  surface.  This  is  because  both  the  thermal 
source  S,  (excess  energy)  and  S,  (bang-gap  energy  distribu¬ 
tion)  are  also  monotonically  decreasing,  and  the  surface  re¬ 
combination  sourced,  =  0.  Fore,  #0,  however,  the  earners 
preferentially  recombine,  and  release  their  energy,  at  the  sur¬ 
face.  This  leads  to  the  formation  of  thermal  spikes.  Further¬ 
more,  since  the  carrier  distribution  develops  a  "sink"  at  the 
surface,  the  source  Sj(x)  will  be  nonmonotonic.  If  the  maxi¬ 
mum  of  S,  is  suffictefitly  strong  so  that  also  the  total 
S  =  5,  +  Jj  is  nonmonotonic,  then  the  resulting  thermal 
distribution  may  develop  a  local  interior  maximum.  (We 
note  that  all  peaks  and  maxima  will  eventually  disappear 
with  time  because  of  carrier  and  thermal  diffusion  ) 

Figures  I  and  2  show  carrier  distributions  for  silicon, 
with  time  (r  /l.)  and  v,  as  parameters.  In  Fig.  1,  the  electron 
lifetime  is  10"*  sec,  in  Fig.  2,  r  —  I  O'1  sec.  The  effect  of  i\  is 
clearly  seen.  The  carrier  distribution  (per  unit  photon  flux)  is 
high  and  monotonic  for  u,  «=  0,  lower  an  nonmonotonic  for 
v,  #0.  For  short  r  (Fig.  I),  the  maximum  is  sharp  and  close 
to  the  surface.  The  maximum  moves  into  the  material  interi¬ 
or  with  increasing  o,  and  higher  t/if.  Because  of  the  short  r 
( r<tf ),  the  distribution  decays  extremely  rapidly  after  the 
end  of  the  pulse  (( >  tf).  For  longer  r  (r>f,,  see  Fig.  2),  the 
maximum  is  much  more  pronounced.  The  concentration  re¬ 
mains  high  even  when  t>tf. 


FIG  I  Fhotocarmr  concentration  in  rub) -irradiated  mIkiio.i.  10‘xc. 
»■  KlO'cm’.r*  10  *  tec.  0  »  J  cm’  nee  ’  For  t  •  0.  nt  *.i  I  k  non¬ 
monotonic  Since  r«  .  nte.f  1  (level v  rnpwlli  sfler  (  i 
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FIG  2.  Photocarhcr  concentration  in  Nd  :  Y AG-irradiated  tilicon,  I, 

=  It)-' ice.  a  -  lO'cm'.r  -  10  tec.  D  -  3  cm- tec'  Fore.  #0.«K<,r  I  in 
nonmonotonic  Since  ra/,.  n0rj  )  is  Mill  high  at » * 


o  i  ?  )  • 

An 

FIG.  4.  Temperature  for  ruby  irradiation  of  heavily  damaged  GaAs 
r  =  10  "  tec.  D  =  I5cm:  tec  '.  p.  -  2  *  I04  cm  sec  a  =  2  -  104  cm  .  i 
=  10 "  tec,  F„^lx  10’  W/cm\ 


Inspection  of  Figs.  3-12  reveals  the  qualitative  features 
discussed  above. 


1.  Ruby  irradiation  of  GaAs 

Figures  3-3  show  temperature  distributions  for  ruby 
laser  irradiation  of  GaAs.  Figures  3  and  4  show  the  effect  of 
v ,.  The  principal  effect  of  large  t>,  in  GaAs  is  to  lower  the 
surface  temperature.  For  the  same  laser  input,  the  surface 
temperature  obtained  for  w,  =  2  X 10*  cm  sec  ’  is  higher  (by 
~50%)  than  the  temperature  obtained  for  v,  =  2  X 10’ 
cm  sec* 1 .  Figures  4  and  3  compare  temperature  distributions 
for  heavily  damaged  and  good-quality  GaAs.  The  principal 
effect  of  the  electron  diffusion  is  to  increase  the  laser  power 
required  for  heating  good-quality  material  in  Fig.  3.  approxi¬ 
mately  by  (I  -  £,/fiv)'~5  (for  the  same  a).  Its  effect  on 
the  width  of  the  profiles  is  negligible;  the  widths  of  profiles 
shown  in  Figs.  4  end  5  are  practically  identical.  This  is  be¬ 
cause  in  good-quality  material,  the  band-gap  energy  (80%  of 
the  photon  energy)  is  distributed  into  the  depth  L=s3  x  10  ' 
cm.  which  is  much  larger  (by  nearly  two  orders  of  magni¬ 
tude)  than  the  depth  for  the  excess  energy  (20%), 
o’~5x  I0*5  cm.  This  makes  the  band-gap  portion  of  the 
energy  ineffective  for  the  heating.  We  should  note  that  for 
intermediately  damaged  GaAs,  more  of  the  available  power 
will  be  used  for  the  heating,  and  fora  '^L,  the  profiles  will 
be  somewhat  wider  than  in  either  the  heavily  damaged  or 
high-quality  materials 


Ft®  i  *  Ruby  kradtstron  of  sdicon 

r  *  10  kc.  D  —  1 5 CHI- ICC '.r.  »  2 x  10’ cm  see a  -  2 x  10* cm  t, 

to  ‘  <«..  F„  *  :  •  io'  w/em-'.  Temperature  profiles  in  amorphous  silicon  are  very 
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FIG  5  Temperature  for  ruby  irradiation  of  GaAt:  r  =  10  *  tec.  fl  =  200 
cm-'  tec  f.  «  2  X  10*  cm  tec  a  »  2  x  10*  cm  f,  «  10 "  tec. 

F„  -  8  x  10’  W/cm 


much  the  same  as  those  in  GaAs.  This  is  because  of  the  simi¬ 
lar  input  parameters  (a  =  2x  10*  cm  '.  r<  10-*  sec)  used  in 
the  two  cases.  Figures  6  and  7  show  temperature  distribu¬ 
tions  in  crystalline  silicon.  The  main  differences  between 
amorphous  and  crystalline  silicon  are  higher  power  require¬ 
ments  and  long  heating  times  for  crystalline  materials.  The 
higher  powers  are  required  because  of  (i)  weaker  optical  ab¬ 
sorption  (o= r  10'  cm1)  and  (ii)  inefficient  utilization  of  the 
band-gap  portion  of  the  energy  for  the  heating  in  good-qual¬ 
ity  materials.  The  heating  times  are  long  because  of  the  influ¬ 
ence  of  the  band-gap  energy  on  temperature  (for  r  >  tf )  at 
long  times.  At  short  times,  the  profile  widths  in  Figs.  6  and  7 
are  ~a'\  as  only  the  excess  energy  heats  the  material  at 
/<  r. 

Figures  6  and  7  correspond  to  a  lightly  damaged 
(r  =  10-1  sec,  D  =  3cm:  sec),  and  good  quality  ( r  =  10' sec, 
D  =■  30  cm'  sec  ’).  silicon,  respectively.  The  main  difference 
between  the  two  profiles  is  that  at  time  lOOr, ;  the  profile 
width  in  Fig.  7  is  narrower  than  the  one  shown  in  Fig.  6.  This 
is  because  at  very  long  L,  L  >  o  ',  in  Fig.  7,  the  heating  is 
entirely  governed  by  the  excess  portion  of  the  electron  ener¬ 
gy;  the  band-gap  portion  affects  the  profile  only  at  very  long 
times,  fir. 

3.  Nd:  YAG  wndmdon  of  silicon 

Several  representative  profiles  are  shown  in  Figs.  8-12. 
For  Nd  :  YAG  in  silicon,  E,  ^0.99 hv  and  the  profiles  are 
therefore  band-gap  recombination  dominated.  Figures  8  and 
9  compare  profiles  for  two  values  of  if,  for  heavily  damaged 
material  (o  =  10'  cm  r  -  lO-**  sec);  r,  *  2  x  10’  cm  sec  1 


FIG.  6.  Temperature  for  ruby  irradiation  of  silicon  r  =  10  hv.  /> 

cm’ tec  '.  v.  =  2x  10' cm  tec  \  o  =  2x  I0‘  cm  '. ».  -  10 "  sec. 

F„  *  2.5x  Iff  W/cm‘. 


FIG  1  Temperature  for  ruby  irradiation  of  good-quality  silicon  -  in 
tec. D  »  Mem'  sec  '.r  ■»  ’  .  10‘cmsev  -  2  •  10'cm  '.r  .  IO’hv 
F„s=  *  -  l(T  Wrcm 
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FIG.  S.  Temperature  for  Nd  :  YAG  irradiation  of  damaged  silicon: 
r  «  IO’sec.0  m  Jcm:  tcc'.v.  =  2x  10’ cm  MC  '.a  «  10' cm  '.I,  =  10  ’ 
sec.  f„«  2.5x10’  W/cm: 


and  2x  104  cm  sec  1  in  Fip.  8  and  9,  respectively.  The  half¬ 
widths  are  ~a  1  >  L  in  both  cases.  The  primary  effect  of 
slow  v,  (v,  <  L )  in  Fig.  9  is  to  suppress  the  interior  maximum 
appearing  in  Fig.  8  for  large  v, .  There  should  be  a  sharp  peak 
at  the  surface  in  Fig.  8.  The  peak  is  unresolved  because  of  the 
numerical  limitations.  The  peak  corresponds  to  thermal 
spikes  (discussed  earlier),  created  at  x  =  0  as  the  result  of  the 
preferential  deposition  of  the  carrier  energy  at  the  surface. 

As  has  already  been  discussed,  the  occurrence  of  ther¬ 
mal  spikes  and  local  interior  maxima  (the  local  maximum  is 
always  lower  than  the  spike)  is  the  consequence  of  fast  sur¬ 
face  recombination.  By  inspecting  the  simulated  profiles,  a 
nonmonotonic  T (x)  occurs  for  v,  sufficiently  large,  v,r>L 
[so  that  maximum  in  n(x,r )  is  large  and  separated  from  the 
surface],  r  not  too  short  (so  that  the  energy  does  not  relax 
instantaneously),  a’ 1  not  too  small,  a'ZL(%o  that  the  band- 
gap  energy  is  not  dissipated  into  too  large  a  volume),  and  the 
excess  portion  of  the  energy  not  too  large  compared  to  the 
band-gap  energy  [so  that  the  source  term  S,(x)  does  not  sup¬ 
press  the  peak  in  S2(x)].  If  these  conditions  art  satisfied,  T(x) 
becomes  nonmonotonic,  with  a  sharp  spike  at  the  surface 
and  a  local  maximum  at  the  distance  0  <x  <  L.  The  magni¬ 
tude  of  the  spike  can  be  estimated  by  using  the  energy  con¬ 
servation  arguments.  For  example,  noticing  that  the  energy 
carried  by  carriers  recombined  at  the  surface  must  be  equal 
to  the  heat  contained  in  the  spike,  we  obtain  that  in  the  limit 
of rapid  t',.  short  r,  and  large  a,  the  height  of  the  spike  is  on 
the  order  of  aorFofu././xc,)''1,  where/,  it  the  char¬ 
acteristic  length  associated  with  the  position  of  the  maxi- 


F1G  4  Temperature  for  Nd  YAG  irradiation  of  damafed  mIkihi 
r  *  tO  "sec.  D  •  }cm;  lec  ’.i>.  *lxlO,cm«c,.o.  10' cm  ',  i 
see. F„>  2  5x10’  W/cm’ 


FIG  10  Temperature  for  Nd  YAG  irradiation  of  uhcon  - --  Id 
Dm  Jem’ mc  r.  m  Jx  10* cm  wc  a  -  10’ cm  '.  t  -  10  « 
F„  -■  tv  to’  W/cm' 
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FIG  1 1  Temperature  for  Nd  :  YAG  irradiation  of  silicon  r  -  10  *  sec. 
D  -  3  cnr  ICC  a.  =  2x  10*  cm  icc  a  -  10'cm  \  /,  «  10  ’  sec. 

F„  =  2xlff‘W/cm;. 


mum  in  An(x,t ).  For  typical  parameters  of  interest,  this 
yields  4  T^v  ^several  103  *K  and  spike  characteristic  width 
tp  10  *  cm.  The  local  maximum  will  initially 

move  deeper  into  the  bulk  (with  the  diffusing  carriers)  until 
it  eventually  disappears,  with  the  characteristic  time 
t  =  maxfr,  f„,  (c,  /Att)L  Jj. 

Figure  10  corresponds  to  light-to-intcrmediately  dam¬ 
aged  material  (a~10'  cm1,  r~lCT7  sec).  The  development 
of  non  monoton  icd  T  (x,t )  is  clearly  seen  (again,  there  should 
be  a  spike  at  x  —  0).  The  temperature  becomes  monotonic  at 
times  ti  f,h ,  ~(c,/4*)Z.  J~10*  sec.  Because  of  the  long  r, 
the  required  power  is  about  twice  higher  than  the  one  for  the 
heavily  damaged  material  in  Fig.  9  (for  otherwise  the  same 
parameters).  Figures  1 1  and  12  show  the  effect  of  the  elec¬ 
tron  diffusivity  D  on  the  profiles.  The  main  effect  of  increas¬ 
ing  D  (at  constant  r,  r  *  10"*  sec)  is  to  suppress  the  interior 
maximum.  Because  of  the  long  r,  the  heating  times  are  very 
long,  the  temperature  is  high  even  at  time  1 2  I  (XV, . 

We  have  also  simulated  profiles  for  a  —  100  and  500 
cm  '.  On  the  scaled  distances  ax,  the  shapes  of  the  profiles 
were  qualitatively  similar  to  those  calculated  for  a  =  10* 
cm  '.  The  required  powers  were  also  much  higher  (approxi¬ 
mately  by  the  ratio  of  a'i  and  the  heating  times  accordingly 
longer  [as  determined  by  the  thermal  time  (c,  /4*r )l  see 

Eq  (23)] 

IV.  DISCUSSION 

In  order  to  facilitate  a  more  comprehensive  under - 


FIG  12  Temperature  for  Nd  YAG  irradiation  of  silicon  r  =  10  *  tec. 
D  =  30 cnr  tec  r.  »  2x  10*  cm  tec  o  »  10‘cm  '.  /,  «■  10  ‘  tec. 
F,-25xl (f  W/em' 


standing  of  our  results.  Figs.  1 3  and  14  show  the  derived 
dependence  of  (i)  pulse  energy  (£3 )  required  to  raise  the  sur¬ 
face  temperature  by  4  7~  1400  *C,  (ii)  time  at  which  the  sur¬ 
face  temperature  reaches  its  maximum  (t„ ),  and  (iii)  time  at 
which  the  temperature  drops  to  one-half  of  its  maximum 
value  (f,,).  on  the  quality  of  the  processed  material  (silicon 
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FIG  13  Lawrnwrsx^lmiuindlnhnithrMlKonturfacriutibclKu  ihc 
melons  temperature  Nd  YAG  t.  »  10  tec.  ruhy  i.  *  lU  x. 
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FIG.  14.  The  dependence  of  (i)  lime  at  which  the  surface  temperature 
reaches  the  maximum  (/..):  (n)  time  when  the  surface  temperature  decays  to 
one-half  of  its  maximum  value  «,  .).  on  material  quality  (silicon). 


only).  The  quality  is  characterized  in  terms  of  the  electron 
diffusion  length  L  (the  corresponding  values  of  r  and  D  and 
a  are  also  indicated).  Q  (Q^Ft, )  is  the  actually  absorbed 
energy;  the  incident  energy  is  nearly  twice  as  high  as  the 
absorbed  energy  because  of  material  reflection.  We  empha¬ 
size  once  again  that  the  values  of  the  material  parameters  are 
some  “average"  representative  values  (this  is  true  particular¬ 
ly  for  the  values  of  a  for  Nd  :  YAG  irradiation  of  silicon), 
and  thus  the  derived  quantities  are  very  approximate  The 
results  can  be  generalized  to  different  laser  parameters. 

A.  Ruby  In  silicon 

The  required  energy  increases  with  L  increasing,  from 
about  0~ 0.25  J  for  amorphous  to  0=4  J  for  good-quality 
material.  Only  the  excess  portion  of  the  energy  is  operativt- 
for  heating  of  the  quality  material  (£,//iv=65%),  and, 
therefore,  Q  saturates  at  long  L,  this  is  seen  in  Fig.  1 3.  With 
the  a' s  and  t,  used,  the  heated  volume  is  thermal  diffusion 
dominated  for  amorphous  material  and  optical  absorption 
depth  limited  for  good-quality  material.  The  surface  tem¬ 
perature  reaches  maximum  at  tm  =  tf.  The  effect  of  the 
band-gap  portion  of  the  energy  is  to  slow  down  the  tempera¬ 
ture  decay  in  lightly  damaged  materials.  At  intermediate 
values  of  L,  L  £  a  ',  the  temperature  remains  high  up  to  the 
time  (c,/4w)L 2  (i.e.,  t =  r).  At  very  long  L  (I>o  '),  how¬ 
ever.  the  band-gap  energy  is  dissipated  into  an  extremely 
large  volume  (depth=several  L  ),  so  that  its  effect  on  tem¬ 
perature  becomes  insignificant.  The  temperature  starts  to 
decay  more  rapidly.  The  decrease  in  (,,,  is  seen  in  Fig.  14  for 
L2  5x  10" ’em. 

B.  Nd :  VAO  In  •Mean 

The  heating  is  nearly  entirely  band-gap  energy  domi- 
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nated.  For  amorphous  and  damaged  materials,  the  optical 
absorption  depth  a'1  >  (<rlk  X )  and  the  profile  width  is  there¬ 
fore  ~a'\  unaffected  by  the  diffusion.  The  surface  tempera¬ 
ture  is  J  rMa  =0a/c, .  In  good-quality  material,  L  becomes 
i.  (O'2  cm.  This  is  larger  than  the  probable  value  of  a  1  in  the 
material,  fl  '~2x  10"J  cm  for  lightly  damaged,  and 
a'~l(y!cmfor  good-quality  material.  Therefore,  the  depth 
of  the  heated  volume  is  =£,  and  the  power  required  to  heat 
the  material  increases  accordingly.  In  Fig.  1 3,  Q  increases  at 
long  L  with  L  approximately  linearly. 

Because  the  heating  is  band-gap  energy  dominated, 
rm„  is  obtained  in  good-quality  materials  at  times  on  the 
order  of  several  r  (tm  in  amorphous  and  heavily  dam¬ 
aged  material).  The  decay  of  the  surface  temperature  is  de¬ 
termined  by  the  thermal  diffusion  time  r,k<  (/„,~flh  , 

=  max[(c,/4*)L 2,  (c,/4*Xl/<*2).  *,]•  Since  the  thickness  of 
the  heated  region  increases  with  L  increasing,  there  is  no 
saturation  of  flo  at  long  L  which  has  been  observed  for  ruby 
irradiation  of  silicon. 

In  summary,  we  have  investigated  in  this  work  the  in¬ 
fluence  of  the  diffusion  of  the  photoexcited  carriers  on  laser 
heating  of  strongly  absorbing  semiconductors.  The  diffusion 
was  found  to  increase  the  powers  required  to  heat  good-qual¬ 
ity  materials  and  to  prolong  the  heating  times.  Surface  re¬ 
combination  was  found  to  create  thermal  spikes  at  the  sur¬ 
face.  The  profiles  could  become  nonmonotonic  as  the 
consequence  of  the  interplay  of  surface  recombination  and 
diffusion  of  the  photoexcited  carriers.  The  strength  of  the 
spike  was  estimated  from  energy  arguments. 

In  practice,  the  existence  of  the  thermal  spikes  may  be 
of  importance  for  annealing  of  localized  defects.  For  a  suffi¬ 
ciently  strong  spike,  it  is  qualitatively  feasible  for  the  heat  to 
dissolve  the  defect.  This  results  in  localized  defect  removal. 
For  practical  utilization,  further  analyses  are  desirable  This 
is  because  of  the  effect  of  both  the  finite  size  'e.g..  finite  thick¬ 
ness  of  the  surface  “dead  layer”),  and  spatial  distribution,  of 
(he  defects  on  the  temperature  profiles.  Further  effects,  not 
considered  here,  include  other  energy  “diffusion”  mecha¬ 
nisms  (such  as  radiative  recombination  of  the  earners),  non¬ 
equilibrium  effects,  and  distributed  damage  such  as  that  cre¬ 
ated  during  ion  implantation.  The  damaged  regions  have 
different  optical,  thermal,  and  electronic  propenies  than  the 
rest  of  the  material  (they  often  absorb  the  radiation  more 
strongly).  Furthermore,  as  the  damage  disappears,  material 
properties  change  during  the  processing.  This  was  not  taken 
into  account  here;  in  the  model  used,  with  the  exception  of 
the  front  surface,  the  material  was  assumed  homogeneous, 
and  its  properties  independent  of  time.  Further  analyses  are 
being  planned  to  account  for  these  effects. 
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APPENDIX  A:  ELECTRON  TRANSPORT  EQUATION: 
SOLUTION  BY  LAPLACE  TRANSFORM  TECHNIQUE 


The  diffusion  equation  we  want  to  solve  is 

=aF{t)txp(-ax)-  +  D  §.  (Al) 

at  r  ox 

We  consider  constant  a,  r,  and  D,  and  a  squre-shaped 
generation  function  F(t ): 

F(t)  =F0,  0  <t<tf 

—  0 ,  otherwise .  (A2) 

The  boundary  conditions  are 

=  t\(/i(0, /)-«„)  (A3) 

for  the  front  surface  and 


G(x,  x',t  —  t ') 

~  I4ri>«  -  nr",  ff.*-11' **(-£) 

and  Ans*n(x,t )  —  n0  is 


(A6) 


An(x, 0  =  dx‘ 

Assuming  the  slab  thickness  l>a  \  the  semiconductor 
can  be  replaced  by  a  semi-infinite  medium,  with  the  result 
that  the  summation  2„  in  Eq.  (A6)  reduces  to  only  one  term, 
namely,  the  term  n  =  0.  Substituting  the  resulting  equation 
into  Eq.  (A7)  and  integrating  it  over  coordinate  x  yields 


J*  G(x,x‘,  t  -  t')aF0exp(-ax) .  (A7) 


n(x  =  /)  =  «„  (A4) 

for  the  back  surface. 

The  initial  condition  is 

n(x,  t  =  0)  =  n„ .  (A5) 

Equation  (Al)  is  a  linear  partial  differential  equation 
for  n(x,t ).  For  v,  =  0,  Eq.  (Al)  can  be  solved  by  using  the 
Green’s  function  technique.  For  v,  #0,  it  is  solved  by  using 
the  Laplace  transform  technique.  In  using  the  Laplace  trans¬ 
forms,  advantage  is  taken  of  the  knowledge  of  the  steady- 
state  solution  to  Eq.  (Al). 

Cam  v,  m  Oi  Groans  function  technique 

For  i>,  =  0,  Green’s  function  is  _ 


An(x,t)=  £  ®*p(  ±  ax)  f  dt ' 

A  (  •  .  I  Jl. 

Xexp  ^  —  —  (1  -  a'L  5)j 
J  ±  x  +  2a  Dt  ’  \ 

A  2 (/>/  ’)’ :  )' 


xerfc 


(A8) 


where  L  —  (Dr)'  -,  erfc(r)  is  the  complementary  error  func¬ 
tion,  and  the  integral  lower  boundary  f0  is 


rn  =  0,  t<tr , 

(A9) 

tn  =  t-t„,  t>tr. 

Integral  (A8)  is  evaluated  by  using  the  partial  fraction 
This  yields 


I)  ,  A  T [«P < ■  - O-J. '>) )  j ; 

+  S’,?1  ’  ”p  (~t~  (l  -  ±  (±  irp)“p|-  °'D' 

(A  10) 

The  second  term  on  the  right-hand  side  of  Eq.  (A  10)  is  evaluated  by  using  the  following  change  of  variables: 
t  =  r2  and  t  —  z2  (All) 

for  the  first  and  second  term  under  the  integral  sign,  respectively.  Performing  the  change  of  the  variables,  each  of  the  two 
terms  transforms  into  an  expression  of  the  type  (o  #0) 

Jexp  ^  -  aV  -  £  e*P<  ±  ±  -jj  +  const . 

Collecting  the  resulting  terms  completes  the  calculation.  The  final  result  is 


(AI2) 


A"(x •'>-  d  ("p  ( -  1  '»■>  7),  £  “*  (“<0' 1  ww) 


(A13) 


Cam  v,  #0:  laplaca  Snaforwt  tachnlqua 

Since  the  steady-state  solution  to  Eq.  (Al)  is  known.  Eq.  (Al)  is  Laplace  transformed  with  respect  to  time.  One  obtains 
(for  the  square-shaped  F{t )] 

exp<  -  «)( I  -  exp(  -  *,)]  -  (1  +  *>  +  -  0  •  <A14) 

where  n{xj)  is  the  Laplac?  transform  of  n{xj )  The  solution  of  Eq.  (A  14)  is 
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4*. L - „> - SLtfg-  - «Pf-i(1+  1 . 

f  I  +  «•  -  alL 1  ^  v,  +  (2>(l/r  +  i))'/J  PV  L  /J 


(A  1 5) 


aD 

[f>(l/r  +  s)]'/l 

The  reverie  Laplace  transform  is  easiest  to  find  by  using  the  well-known  results  for  the  partial  expressions  appearing  in 
Eq.  (A  13),  and  making  time  convolutions  of  the  partial  Laplace  transforms.  That  is. 


where 


re 

foh*FAt),  /i(i)/2(s)— */"(t ) ,  /(j  +  c)—F(,t )  exp(  —  c» ) , 


(A16) 


a  \  a 

The  partial  transform  used  in  evaluating  the  inverse  Laplace  transform  of  Eq.  (A15)  are  of  the  type 


M 


exi 


P (-«) 
J(r  +  b) 


-HI -expj  -b(t-a)])H(t-a\  a> 0 
0 


—  (1  —  exp(  —  bt )], 

D 


a<0 


and 


exi 


■p(-«l,a) 


_L_exp(  _  ^-iexpfab  +  i  v)erfc  +  bt  >'J ) . 


6  +  r,/l  («)' 

In  Eq.  (A  17),  H(t)  is  the  step  function.  H(t)  =  1  for  r  >  0  and  zero  otherwise. 
Using  Eqs.  (A16)-(A18),  the  inverse  transform  of  Eq.  (A1S)  is 

! )  -  —  ar^'T  -  a>)[ 1  —  eip^  —  i-(l  J‘ 


(A17) 


(A  18) 


(4 Dty 


“■(ail 


x(l-exp(-  i-^<l-a*I*))-[l-exp(-  -^-^.(1  -a2L *))  ]/f(r - t‘ ' -  )  (A19) 

Integral  (A  1 9)  is  evaluated  by  using  the  same  changes  of  variables  as  those  used  for  evaluation  of  integrals  in  the  Green's 
function  expression.  Eq.  (A8)  and  the  following  discussion.  It  is  convenient  to  rewrite  Eq.  (AI9)  in  the  form 

X  £/(f  )  { 1  -  g(/  ■)  exp  (  -  i ( 1  -  a JL  *))  -  [  1  -  g(i  ”)  exp  (  -  ( 1  -  aJI J))  ]//  (r  -  r  ’  -  /, )]) , 

(A20) 

(A21) 
(A22) 


where 


/c ;) h -  "■“pfe<’+"'))“fel2(Sp+  "-fen)' 

g<r)*exp(i(l  -a,L*)). 

The  solution  of  Eq.  (A20)  is 

dn(x.O  =  i  -d-a**1))-  ^-^[Ffr^-IGfr'l-GfO)) 


(23) 


where 
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(A24) 


(A25) 


dt '/{!')  +  const, 


G(t)  =  jdt  */(/ *(/  )  +  const 


The  first  term  in  the  integral  on  the  right-hand  side  of  Eq.  ( A24)  is  evaluated  by  using  the  change  of  the  variables) '  =  z2. 
The  resulting  integral  is  of  the  type  of  Eq.  (A12).  The  second  term  under  the  integration  sign  is  ev  aluated  by  first  using  the 
partial  fraction,  and  then  making  the  change  of  the  variables,  t '  *  z  2,  in  the  two  resulting  integral  expressions.  The  integral  is 
of  the  type 


1  ±**>(7  +  j)  erf(" ±  7$ 


(A26) 


Integral  (A25)  is  evaluated  by  using  the  same  procedure.  In  fact,  the  only  difference  between  the  integrals,  Eqs.  (A24)  and 
(A23),  is  the  replacement  of  r1  in  Eq.  (24)  by  a*D  in  Eq.  (A23). 

f(t.  r)g(t,  ry+f(t,  \/Da})  (A27) 


The  mathematics  is  then  straightforward,  though  rather  tedious.  After  having  collected  and  combined  all  terms,  the  final 
result  is 


An(x,  t)- 


aF0r 


2(1  -Z.V) 
+ 


zM 


(4  ory 


(±  1  +Da/v,) 
<*.-»  (1  -Da/v,) 
2(1  -  L  W) 


(aD/u,  -  1X1  -  o\r/D) 


m 

+  a(Di 


■>)'" 


(A28) 


It  is  easy  to  verify  that  Eq.  (A28)  reduces  to  Eq.  (A  1 3)  in  the  limit  of  v,  =  0.  In  the  steady  state,  Eq.  ( A28)  reduces  to  the 
well-known  expression 


lim  4«i(x,  oo )  = 

v-« 


arF0 
1  -  a2L 2 


exp(  -  ax)  - 


u,  +aD 

v,  +  L/r 


exp 


(A29) 
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'"Other  mechanisms  include  photon-phonon  interaction  and  free  absorp¬ 
tion.  We  have  made  estimates  of  the  free^arrier  absorption  and  found  that 
for  laser  intensities  of  interest  the  effect  is  small  al  the  ruby  laser  wave¬ 
length.  However,  it  may  become  large  at  the  Nd  :  YAG  laser  wavelength 
The  photon-phonon  and  free-elect ron  absorptions  were  not  taken  into 
account  is  this  work 

'  'The  pulse  lengths  (t,)  are  on  the  order  of  10'*- Kk  ’  sec.  The  thermal  diffu¬ 
sion  time  (time  it  takes  the  thermal  energy  to  diffuse  a  distance  a  '  s»s> 
from  the  semiconductor  surface)  is  sc,/«r  slO'-IO^  sec  for 
a  m  IC-IO’  cm  ',  respectively.  The  lime  it  lakes  the  electron  (provided 
that  electron  lifetime  is  long  enough  that  the  electron  does  not  recombine 
earlier,  L>o'')  to  traverse  a  distance  cr'  is  r  =  \/Da'Z  (O'*  sec  for  the 
typical  value  of  DS  10  cm1  sec1. 

I3S.M.  See,  Physics  of  Semiconductor  Dntct!  (Wiky-lnterscience,  New 
York,  1969). 

"It  is  well  known  that,  for  Nd  :  YAO  irradiation  a  increases  rapidly  with 
temperature.  For  single  crystal  silicon,  as  10  cm  1  at  300  *K  and  in- 
cremes  by  approximately  two  orders  of  magnitude  at  the  melting  tempera- 
lure.  This  was  not  taken  into  account;  the  values  used  are  some  represen¬ 
tative  avenge  values. 
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APPENDIX  B 


LASER  ANNEALING  OF  LOW  FLUENCE 
ION  IMPLANTED  SILICON 
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Later  annealing  of  low-fluence  ion-implanted  ailicon 

S.  Prusatn 

TRW  Semiconductors,  Lawndale.  California  90260 

W.  vonderOhe 

TR  W  Defense  and  Space  Systems  Group,  One  Space  Park,  Redondo  Beach.  California  90278 
(Received  29  October  1979;  accepted  for  publication  25  February  1980) 

Ruby  laser  pulses  of  1.5-2. 5  J  cm1  were  shown  to  be  effective  in  dissolving  defect  nuclei  in  1 100) 
and  1 1 1 1 1  silicon  surfaces  implanted  with  10'*— 10,s  B  cm  2  and  (1-5)  X  10u  P  cm  2.  No  melting 
and  only  partial  electrical  activation  appear  to  occur.  A  sequential  moderate  thermal  anneal,  10 
min  at  1000  *C,  results  in  full  electrical  activation.  A  |  Ill  |  silicon  surface,  made  amorphous  by  an 
ion  implantation  of  5  X 10'4  P  cm 2  was  electrically  activated  and  completely  cleared  of  defect 
nuclei  without  apparent  melting  by  preceding  the  laser  anneal  by  a  two-step  thermal  anneal  (30 
min  at  550 ‘C,  10  min  at  1000  *C).  The  mechanism  of  pulse  laser  annealing  of  nonamorphous  ion- 
implanted  surfaces  is  discussed. 

PACS  numbers:  78.S0.Ge,  6l.70.Tm,  79.20  Ds,  6l.80.Jh 


I.  INTRODUCTION 

Ion  implantation  has  been  found  to  result  in  the  intro¬ 
duction  of  electrically  active  impurities  in  an  extremely  uni¬ 
form  and  reproducible  manner.  Concomitant  with  this  de¬ 
sirable  characteristic,  ion  implantation  leaves  a  residue  of 
lattice  defects.  These  primary  defects  consist  of  impurity 
atoms  on  nonlattice  positions,  and  the  presence  of  many  sili¬ 
con  vacancies  and  interstitials.  A  relatively  moderate  ther¬ 
mal  anneal,  1000  *C  for  10  min,  results  in  relocation  of  the 
impurity  atoms  to  lattice  positions,  as  shown  by  their  electri¬ 
cal  activation.  There  is,  however,  a  considerable  rmidiium  of 
lattice  damage  remaining.  These  secondary  defects'  have 
been  studied  by  TEM  for  both  light  and  heavy  implantations 
and  have  been  shown  to  consist  primarily  of  dislocation 
loops,  with  and  without  enclosed  stacking  faults.2  For  heavy 
implantations  the  Rutherford  ion-backscat taring  and  ion- 
channel  technique  has  been  used  to  show  the  presence  of 
residual  defects.2  A  method  developed  by  the  author4  uti¬ 
lizes  a  two-step  technique.  By  growing  a  thermal  oxide,  ex¬ 
trinsic  dislocation  loops  are  found  to  expand  to  a  size  of 
several  microns.  Using  chemical  etching  techniques  the  in¬ 
tersection  of  these  ternary  defects  with  the  surface  can  be 
observed  by  optical  or  scanning  electron  microscopy.  This 
method  was  applied  to  boron  and  phosphorus  over  the  range 
lO'MO”  cm 2.  A  thermal  anneal  of  1  h  at  1 100  *C  in  dry  N, 
was  effective  for  10M  cm'2,  but  not  for  10' 3  cm'2  tuna  cue  in 
dissolving  the  defect  nuclei  introduced  by  the  ion  isaplants- 
tion,  to  timt  i  mtupqwm  oxidation  (Mid  to  mult  n  tkt 
dcvdopHKBt  Of  CXpHMltd  ^g)ocMiOQ  hovtVCf, 

causes  a  significant  redistribution  of  the  implanted 
impurities. 

In  this  study  a  ^-switched  ruby  laser  was  used.  A  re¬ 
view  of  the  literature,  as  wen  as  a  recent  review  article,  indi¬ 
cates  that  this  laser  has  only  been  effectively  used  when  the 
implant  ilucnccs  were  high  enough  to  convert  the  surface  to 
an  amorphous  condition.3  Under  such  conditions  it  was 
found  that  the  laser  beam  melts  the  crystal  to  a  depth  some¬ 
what  greater  than  that  of  the  implanted  profile.  This  molten 
layer  recrystallizcs  from  the  undamaged  substrate  generat¬ 


ing  an  epitaxially  grown  defect-free  single-crystalline  silicon 
layer  containing  the  implanted  impurity  atoms  in  normal 
lattice  sites.3  Thus  both  defect  elimination  and  electrical  ac¬ 
tivation  of  the  dopant  atoms  are  simultaneously  carried  out. 
During  the  time  the  implanted  region  is  molten,  the  impurity 
atoms  diffuse  rapidly  in  the  liquid,  resulting  in  a  significant 
redistribution  of  the  implanted  dopant  atoms.  It  has  been 
inggrstad  that  the  trine  available  for  annealing  with  Q- 
s witched  ruby  I  seen  is  too  short  for  an  unmolten  damaged 
layer  to  be  annealed  with  a  single  pulse  and  that  melting 
followed  by  liquid  phase  epitaxy  is  essential  for  effective 

ContinuouB  At  laser  annealing  was  applied  by  Gat  et 
ml  *  to  sskeon  surfaces  implanted  with  fluences  of  As  heavy 
enough,  >X  I0l4cm  2at  100  keV,  to  convert  the  surface  to  an 
amorphous  condition.  Melting  did  not  occur  and  thus,  also, 
no  redistribution  of  the  implanted  As.  Fewer  defects  re¬ 
named  than  were  found  for  a  single-step  thermal  anneal  ( 30 
nun  at  1000  T),  but  complete  elimination  of  defect  nuclei 
did  not  occur. 

In  this  study  we  chose  B-implant  fluences  of  I0*4. 

Sx  I014,  and  10”  cm  2.  and  P -implant  fluences  of  I014  and 
Sxl0,4cm  1  tato  both  1 100)  and  1 1 1 1 1  silicon  surfaces 
One  implantation,  that  of  5  X I014  P  cm 1  into  a  1 1 1 1 1  sur¬ 
face.  resulted  in  amorphism.  For  the  rest  of  the  wafers,  the 
surfaces  rstsiasd  their  aiagk  crystallinity 

For  these  lower  flamer  implants  it  was  found  that  laser 
energies  which  were  described  ia  the  Htcrsture  as  resulting  in 
melting  of  the  surface  and  completely  activating  the  im¬ 
planted  impurity  atoms  bene  foiled  to  beat  the  silicon  sur¬ 
faces  sufociently  to  permit  more  than  a  fraction  of  the  im¬ 
planted  dopant  moms  to  migrate  to  substitutional  positions. 
Nevertheless  the  temperature  distribution  within  the  im¬ 
planted  layer  was  such  that  complete  solution  of  implant 
defect  nuclei  occurred.  By  combining  a  laser  anneal  with  a 
short  thermal  anneal  it  was  possible  to  achieve  both  defect 
diminatioo  as  well  as  complete  electrical  activation.  Such  a 
process  would  avoid  impurity  redistribution,  permitting  the 
dopant  distribution  to  remain  as  implanted 
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TABLE  1.  Ei pennies taj  conditions. 


Wafer  numbers 

Type 

Oriental  ion 

Implant 
specie* ' 

Implant 
flucnce 
(I014  cm  ’) 

IDA.  II3B 

II 

000) 

n» 

1.0 

II4A.  II4B 

P 

(100) 

31, 

1.0 

IDA.  IDB 

H 

(100) 

II, 

50 

IISA.  USB 

P 

000) 

3Ih 

5.0 

II7A,  II7B 

a 

000) 

II. 

10.0 

213A. 2I3B 

IV 

(III) 

II. 

1.0 

2 DA,  2I4B 

p 

(III) 

31, 

1.0 

2I5A.2I5B 

A 

(III) 

II, 

5.0 

2 ISA,  2ISB 

p 

(III) 

Mr 

5.0 

2I7A,  2I7B 

A 

(III) 

11. 

10.0 

'Boros  implaataliooa  «i  M  hV;  ptapfcorus  implmmirmi  it  X)  trV 


H.  EXPERIMENTAL  INVESTIGATION 

Twenty  wmfen  woe  subjected  to  implsnUtkm  with  bo¬ 
ron  and  phosphorus  ions.  The  orientations,  conductivity 
types,  and  resistivities  of  the  wafers,  together  with  the 
fluences  and  implant  energies  are  given  in  Table  I.  The  im¬ 
plant  energies  were  chosen  to  give  similar  distributions  for 
the  B  and  P  ions.  The  wafers  were  commercially  polished 
and  had  resistivities  in  the  range  of  1-4  H  cm.  After  implan¬ 
tation  the  wafers  were  examined  for  the  presence  of  amor- 
phism.  Thu  was  determined  by  the  presence  c^a  characteris- 
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tic  “milky"  appearance  of  the  implanted  area  as  compared  to 
the  adjacent  areas  protected  from  implantation  by  the  wafer 
holder.2  Wafers  2I6A  and  216B  exhibited  such  amorphism. 

After  implantation  the  series  A  wafers  were  subjected 
to  an  identical  laser  annealing  process  and  the  series  B  wafers 
were  halved,  each  half  being  subjected  to  a  separate  laser 
annealing  process.  Figure  1  describes  the  experimental 
procedure. 

The  laser  annealing  was  carried  out  with  a  ^-switched 
ruby  laser  operated  in  the  multimode  condition.  The  pulse 
duration  was  13  nsec  and  the  beam  size  was  approximately  1 
cm  in  diameter.  The  beam  was  projected  through  a  rectangu¬ 
lar  aperture  4.4  x  7.9  mm.  For  the  A  wafers,  large  areas, 
approximately  one- fourth  of  the  surface  of  each  2-in.-diam 
wafer,  were  pulsed  with  a  single-laser  energy.  This  was  ac¬ 
complished  by  using  pulses  of  the  same  energy,  as  the  wafer 
was  repeatedly  translated  a  distance  equal  to  the  width  or 
height  of  the  projected  beam.  Figure  2  illustrates  the  laser 
pulse  distribution.  The  laser  energies  utilized  were  1 .0, 1 .3, 
2.0,  and  2.3  J  cm  2.  After  laser  annealing,  all  surfaces  were 
examined  by  high-power  optical  microscopy.  The  ripple  ap¬ 
pearance  characteristic  trf' surface  melting7  was  not  found.  In 
the  center  of  each  wafer  there  remained  a  small  reference 
area  which  was  not  subjected  to  laser  annealing.  The  laser 
annealed  areas  were  sufficiently  large  to  permit  a  sheet  resis¬ 
tance  value  to  be  obtained  directly  without  introducing  cor¬ 
rection  factors. 

Table  II  lists  the  results  of  the  sheet  resistance  measure¬ 
ments  following  the  laser  annealing.  A  minimum  of  six  read¬ 
ings  were  taken  for  each  annealed  area  and  the  average  was 
taken.  For  the  two  lower  energies  there  was  considerable 
scatter  while  for  the  two  higher  energies  all  readings  were 
within  ±  5%  of  the  mean.  This  would  appear  possibly  to  be 
due  to  inhomogeneities  within  the  laser  beam.  For  the  higher 
energies,  lateral  thermal  conduction  could  smooth  out  the 
annealing  effect. 

After  measuring  the  sheet  resistances  of  the  laser  an¬ 
nealed  areas  all  wafers  were  subjected  to  a  thermal  anneal  of 
10  min  at  1000  *C  in  a  N,  atmosphere.  The  results  of  sheet 
resistance  measurements  repeated  on  the  same  laser  an¬ 
nealed  areas  are  found  in  Table  III. 
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FIG.  2.  DiMnbMioaaflMerpuho.  Area  I  <1.01  cm  ’).  arts  110  5  J  cm  !). 
a It*  III  (2.0  J  cm1),  uta  IV  <2.5  I  cm  ') 

As  seen  in  Fig.  I,  the  B  series  wafers  were  halved,  one- 
half  designated  as  B1  and  the  other  as  B2.  The  B1  halves 
were  subjected  to  laser  annealing  at  nine  discrete  areas.  Each 
area  was  4.4  X  7.9  nun,  the  size  of  the  aperture  through 
which  the  beam  was  projected.  The  laser  energies  used  were 
0.3, 1.0, 1.3, 2.0,  and  2-3  J  cm'2  with  a  single  pulse  as  well  as 
0.3(2X),0.3(4x),  1.0  (2  X),  and  1.0  Jem2  <3x).Each 
discrete  laser  annealed  area  was  surrounded  by  a  nonan- 
nealed  area  and  after  oxidation  and  Sirtl  etching  (he  laser 
pulse  boundary  (LPB)  was  examined  by  optical  microscopy 
and  compared  to  the  ion-implant  boundary  (I1B).  between 
the  area  of  the  wafer  covered  by  the  wafer  holder  and  the 
area  of  the  wafer  subjected  to  implantation  without  laser 
annealing. 

Wafers  1 16b  and  216B  received  identical  implantations 
of  5 xlO14  Pent'2,  but  only  2I6B  with  a  (lilt  orientation 
exhibited  amorphism.  The  wafer  halves  1 16B2  and  216B2 
were  subjected  to  a  preanneal  treatment  of  I  h  at  SXFC  to 
permit  slow  transformation  of  the  amorphous  layer  to  single 
crystallinity.  This  was  conllrmed  by  die  disappearance  of  the 
milky  surface  on  216B2.  All  B2  halves  were  then  annealed 
for  10  min  at  1000  *C.  Sheet  resistivity  measurements  were 
carried  out  and  are  given  in  Table  III  under  the  heading  of 
zero-laser-beam  energy  / 

The  B2  halves  were  then  subjected  to  laser  annealing  at 
five  discrete  areas,  using  single-laser  pulses  of  0.3, 1.0, 1.3, 


2.0.  and  2.3  J  cm  2.  Each  discrete  laser  annealed  area. 

4.4  x  7.9  mm,  was  surrounded  by  a  nonannealed  area  After 
steam  oxidation  and  Sirtl  etching,  the  surfaces  were  exam¬ 
ined  by  optica)  microscopy. 

III.  EXPERIMENTAL  RESULTS 

Examination  of  the  results  listed  in  Table  II  shows  that 
in  every  case  there  is  s  significant  increase  in  electrical  activ¬ 
ity  of  the  implanted  ions  in  going  from  1.0  to  1 .3  i  cm  2  and 
from  1.3  to  2.0  J  cm'2.  In  going  from  2.0  to  2.5  J  cm'2,  the 
change  in  Ku  appears  to  be  much  smaller.  At  first  it  might 
appear  that  this  is  a  result  of  an  approach  to  complete  activa¬ 
tion.  Following  the  laser  anneal  with  the  moderate  thermal 
anneal  of  10  min  at  1000  *C  results  in  a  very  significant  de¬ 
crease  in  sheet  resistance  for  all  specimens  for  all  laser  an¬ 
neals,  as  is  shown  by  comparing  Table  III  with  Table  II.  This 
is  illustrated  in  Fig.  3,  in  which  the  sheet  resistances  for  laser 
annealing  alone  and  for  combined  laser  and  thermal  anneal¬ 
ing  are  compared. 

The  question  arises  whether  this  drop  in  sheet  resis¬ 
tance  is  due  to  the  activation  of  the  implanted  impurity  or  to 
the  redistribution  of  the  implanted  atoms  over  a  larger  vol¬ 
ume,  resulting  in  a  higher  average  mobility  for  the  carriers. 
In  the  appendix  the  expected  decrease  in  sheet  resistivity  for 
the  B-implanted  wafers,  due  to  the  redistribution  occurring 
fen  an  anneal  of  10  min  at  1000  *C,  is  calculated  to  be  7. 5%. 
This  was  checked  by  subjecting  a  wafer  implanted  with 
I X I0'4  B  cm'2  at  50  keV  to  a  succession  of  five  2-min  an¬ 
neals  at  1000  *C.  The  sheet  resistivities  measured  after  each 
anneal  were  found  to  remain  within  the  calculated  range.  We 
therefore  conclude,  for  the  fluences  studied,  that  laser  an¬ 
nealing  by  itself  is  incapable  of  bringing  about  full  electrical 
activation  of  the  implanted  ions. 

A  measure  of  (he  effectiveness  of  laser  pulses  in  elimi¬ 
nating  implant  damage  is  the  comparison  of  the  implanted 
and  laser  annealed  surfaces  with  that  of  an  unimplanted  and 
unpulsed  surface  after  oxidation  and  Sirtl  etching.  Where 
them  are  identical  we  can  reasonably  assume  that  the  laser 
annealing  eliminated  the  defect  nuclei  introduced  by  the  ion 
implantation.  For  the  B1  specimens,  single-laser  pulses  of 
0.5  J  cm  2  failed  to  show  any  reduction  in  defect  density.  In 
every  case  a  laser  pulse  of  2.0  or  2.5  J  cm'2  resulted  in  com- 
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TABLE  III  Sfcctt  mutant*  in  tl  AD  of  taxr  ■■■*»lwl  area*  (dIIowmi*  iktrmal  aaaral  of  10  MIN  al  1000  *C  m  N .  ambtcm 


Laacf  beam  energy 

(l/cm  ) 


SpKMMa 

0 

10 

15 

20 

2  5 

II) 
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65X 
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W) 
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WO 
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ISS 
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1)1 
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1  IT 
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99 
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M3 
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I4S 
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109* 
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pfete  chminarion  of  the  ion-implant  defects,  la  most  of  the 
Bl  specimens  (213, 214. 114. 115, 116.  end  117)  all  detects 
were  eliminated  by  a  1.5  J  cm"2  pulse,  for  the  others  (113, 
215, 216  217)  a  slight  defect  residue  remained  after  laser 
puking.  For  a  laser  pulse  of  1 .0  J  cm'2  there  ares  a  transition¬ 
al  structure  between  no  annealing  and  complete  annealing. 
This  structure  was  alto  found  at  the  laser  pdee  boundary 
(LFB)  of  the  higher-energy  laser  annealed  stem 

Figure  4  illustrates  the  laser  pulse  boundary  of  speci¬ 
men  1I3BI  for  1.3 1  cm1.  The  aonpulsed  surface  is  charac- 
by  §  light  of  oaadslkwi  baht 

(OSF).  the  puhed  ares  is  relatively  dear,  and  the  transitional 
npon  m  CMnctcraeo  oy  «  nig pct  oonccmranom  ov  uor 
correaptmdmg  to  that  found  far  the  lower  lamrenetgy  of  1.0 
J  cm'2.  Figura  5  and  6  illustrate  the  LPB  for  specimen 
U7B1  far  1.51  cm  *.  In  Pig.  5  we  see  the  laser  annealed  area 


ssee  J  Aow  «ws.,  Vot.  si.  no  7,  July  iseo 


free  of  OSF  on  the  right  and  the  transition  zone  containing  a 
heavy  concentration  of  OSF  with  a  light  background  on  the 
left.  In  Pig.  6  the  transitinn  sons  is  on  the  right  and  the 
unannealed  area  containing  a  heavy  background  together 
with  a  heavy  concentration  of  OSF  is  found  on  the  left. 

In  this  investigation  we  also  studied  the  effect  of  multi¬ 
ple  laser  pulses  of  0.5  and  1.0  Jem'2  with  a  1 -min  interval. 

Since  a  laser  pulse  of  1 .0  J  cm'2  results  in  a  significant 
reduction  of  the  defect  concentration  for  various  implanta¬ 
tions  while  1 .5  J  cm~2  was  elective  in  eliminating  all  defects 
in  moat  cases,  it  appeared  reasonable  that  multiple  pulses  of 
1.0  J  car2  might  be  equivalent  to  a  single  pulse  of  higher 
energy.  It  was  found,  however,  that  multiple  pulses  of  0.5 
and  1.0J  cm2  did  not  reduce  the  defect  concentration  to  any 
significant  extent  below  the  level  found  for  a  single  pulse  of 
the  tame  energy. 

Following  the  thermal  anneal  of  the  A  series  wafers, 
these  too  were  subjected  to  steam  oxidation  and  Sirtl  etch¬ 
ing.  Microscopic  examination  showed  results  similar  to 
those  obtained  with  the  Bl  series,  confirming  our  prior  ex¬ 
perience  that  an  anneal  of  !000*C  for  10  min  has  little  or  no 
effect  in  dissolving  ion-implant  defect  nuclei.4 

The  aaneahag  procedure,  10  min  at  1000*Cindry  N}, 
utilized  for  the  B2  wafer  halves  prior  to  laser  pulsing  con¬ 
verts  the  ioa-implant  damage  into  secondary  defects.1  These 
consist  of  small  dislocation  loops,  probably  resulting  from 
the  agglomeration  of  the  primary  defects,  the  vacancies  and 
the  interstitial  atom  generated  in  the  implantation  process. 
The  secondary  defects  appear  to  be  much  more  stable  than 
the  primary  defects  m  shown  by  their  greater  resistance  to 
solution  by  laser  annealing  Of  the  ten  specimens,  only  one. 
216B2,  exhfoitedcooapieledsfect  removal  with  laaer  anneal¬ 
ing.  For  all  the  others  there  was  some  reduction  in  defect 
concentration,  but  even  for  the  highest  laser  pulse  energy 
need,  2.5 1  cm*2,  there  was  a  very  significant  residue  of 
defects. 

The  1 1 11 1  oriented  wafers,  2I6A  and  2I6B,  which  re¬ 
ceived  5  X 1014  Pcm'2  isapianution  were  the  only  ones  of  the 
20  studied  which  exhibited  amorphism  after  implant*  >  ion. 
The  ( I00|  oriented  wafers,  I I6A  and  I 16B,  which  rwcived 
the  same  implantation,  did  not  .  Specimens  1 1>B2  and  216B2 
received  identical  two-step  thermal  anneals  prior  to  laser 
pulsing.  The  first  step,  40  min  at  $50  *C,  was  chosen  to  pet  - 
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FIG  5  Lasei  pulse  boundary  foi  specimen  I I7B1.  junction  of  laser  an 
nealrd  and  iransmon  area  pulsed  area  on  the  right.  I  x  10''  B  cm  .19 
J  cm  ’  lasei  energy.  Sinl  etch  10  Sec.  magnification  250  » 
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FIG  4  Laser  pulse  boundary  for  specimen  1 1 3BI.  pulsed  area  on  the  right. 
1  X  10'4  B  ern  \  1.5  J  cm  :  laser  energ>.  Sirtl  etch  10  Sec.  magnification 
100X 


FIG  6  Laser  pulse  boundary  for  specimen  1 1 7BI.  junction  of  unannealed 
and  transition  area  Pulsed  area  on  the  right.  I  x  10,s  cm  \  15  J  cm  ;  laser 
energ).  Sirt I  etch  10  Sec.  magnification  250  >  . 


mit  slow  regrowth  of  the  amorphous  layer  into  single-crystal 
silicon.  Examination  of  specimen  2 16B2  after  this  step  did 
indeed  show  that  the  implanted  surface  had  lost  its  milky 
appearance  and  could  not  be  distinguished  from  the  adjacent 
nonimplanted  surfaces.  Sheet  resistivity  measurements  gave 
values  of  392  Cl  /□  for  1 1 6B2  and  264  Cl  AD  for  2 1 6B2,  agree¬ 
ing  with  work  of  Miyao  ei  al  *  that  only  partial  activation  is 
obtained  for  550  "C  anneals 

The  second  step,  10  min  at  1000  *C.  was  a  standard 
thermal  anneal  used  to  produce  electrical  activation  of  the 
implanted  impurities.  The  two-step  process  has  been  shown 
to  reduce  significantly  the  implant  damage  for  heavier  im¬ 
plant  fluences  v 

We  found  that  the  use  of  a  two-step  thermal  anneal  for 
this  amorphous  surface  permits  us  to  dissolve  the  ion-im¬ 
plant  defect  nuclei  with  a  laser  pulse  energy  of  2  0  J  cm 


FIG  7  Scanning  electron  micrograph  of  neutron  irradiated  silicon  surface 
espooed  to  a  rub>  laser  pulse  of  2  0  J  cm  magnification  .1000  x 
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This  last  result  suggests  that  the  technique  of  com'c;ned 
laser  and  thermal  annealing  which  was  developed  for  low- 
fluence  implants  as  described  here  might  also  be  effectively 
applied  to  heavy  implants,  where  the  surface  is  converted  to 
an  amorphous  condition  and  therefore  prone  to  melting  and 
concomitant  impurity  redistribution  when  subjected  to 
pulsed  laser  annealing  By  introducing  a  two-step  thermal 
anneal,  the  amorphous  surface  is  solid  phase  epitaxially  re¬ 
grown  without  impurity  redistribution.  Subsequent  laser  an¬ 
nealing  results  in  the  dissolution  of  the  implant  defects  with¬ 
out  melting  and  the  resultant  impurity  redistribution. 
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IV.  DISCUSSION 

The  results  of  this  study  suggest  that  the  energy  release 
in  the  interaction  of  a  pulsed  ruby  laser  with  a  single-crystal¬ 
line  surface  containing  discrete  implant  defects  occurs  pri¬ 
marily  at  the  defects,  creating  localized  hot  spots  which  re¬ 
sult  in  the  dissolution  of  the  defect  and  in  the  heating  and 
electrical  activation  of  the  area  adjacent  to  the  defect.  The 
partial  activation  of  the  implanted  impurities  is  explained  by 
the  dejecta  bring  concentrated  closer  the  w*f^*r*  the 
impurities. 10  This  phenomenon  was  dearly  illustrated  in  a 
separate  unpublished  study  of  the  interaction  of  Q-swhch 
ruby  laser  pulses,  identical  to  the  ones  used  in  this  study, 
with  silicon  which  has  been  subjected  to  neutral  bombard¬ 
ment.  In  that  study  we  found  clear  signs  of  the  occurrence  of 
localized  melting  at  the  silicon  surfaces  in  the  vicinity  of 
defects.  Fig.  7.  We  did  not  detect  such  meh  spots  for  ion- 
impbntcd  and  laser  annealed  silicon. 

The  results  obtained  suggest  that  melting  does  not  oc¬ 
cur  in  the  laser  annealing  of  a  single-crystalline  implanted 
silicon  surface.  This  belief  is  based  on  the  occurrence  of  par¬ 
tial  electrical  activation  by  laser  annealing,  the  partial  solu¬ 
tion  of  defects  with  lower  laser  energies,  and  the  absence  of  a 
rippled  surface  characteristic  of  melting. 

The  interaction  of  a  laser  pulse  with  a  semiconducting 
crystal  is  well  understood.  For  each  photon  absorbed  we 
generate  a  hole-electron  pair.  These  have  a  kinetic  energy  as 
well  as  a  potential  energy  equal  to  the  band-gap  energy.  The 
kinetic  energy  it  transferred  almost  immediately  to  the  crys¬ 
tal  lattice  as  heat,  while  the  potential  energy  is  released  at 
such  a  time  and  place  as  recombination  of  the  hole-electron 
pair  occurs. 

We  recognize  that  the  intensity  of  a  laser  pulse  I  in  pene¬ 
trating  a  distance  x  from  the  surface  of  a  homogeneous  crys¬ 
tal  is  given  by  the  relation 

where  /„  is  the  initial  intensity  and  a  is  the  absorptivity. 

The  concentration  of  excess  hole-electron  pairs  gener¬ 
ated  is  proportional  to  the  rate  of  absorption. 

[a]  »[/>]»  —  k  m  kal (x). 

dx 

The  surface  volume  contains  not  only  the  implanted 
impurities  but  also  a  high  concentration  of  crystal  defects 
which  can  act  as  recombination  centers  for  the  hole-electron 
pairs  generated  by  the  laser  pulse.  This  permits  the  genera¬ 
tion  of  high  temperatures  at  defect  sites  by  the  diffkisioa  of 
excess  electrons  and  holes  from  adjacent  areas.  A  require- 
meat  for  this  model  to  be  applicable  is  that  thf  rate  of  diff\»* 
non  of  botes  md  cfoctroas  Co  tbs  detect  isoosbiMftos  ocn* 
ten  exceed  by  a  significant  amount  the  rate  of  thermal 
cSdihseioo  firoioc  tSss  ctctsct ,  ^fbis  coo^td  bos  bee®  ^Sevc&oped  in 
ooitscdeebbte  ® 

Another  model  views  the  defect  site  as  important  in  the 
generation  of  hole-electron  pairs. 

The  abaorptivttya  is  a  Auction  of  the  crystalline  perfec¬ 
tion  of  the  silicon.  The  value  for  amorphous  silicon  exceeds 
that  for  crystalline  silicon  by  an  order  of  magnitude.12  This 
suggests  that  the  absorption  of  laser  energy  at  a  defect  site  is 


much  higher  than  in  adjacent  undamaged  silicon.  According 
to  this  view,  we  would  generate  high  concentrations  of  holes 
and  electrons  around  individual  defect  sites.  We  recognize 
that  the  recombination  rate  is  proportional  to  the  hole-elec¬ 
tron  concentration  product,  i.e.,  R  =  r[n][p]  An  increase  in 
the  electron  and  hole  concentrations  of  an  order  of  magni¬ 
tude  results  in  an  increase  in  the  recombination  rate  of  two 
orders  of  magnitude.  The  rate  of  diffusion  of  excess  holes  and 
electrons  away  from  the  defect  site  where  they  were  generat¬ 
ed  will  be  proportional  to  their  concentration  gradient  and 
hence  to  their  concentration.  This  is  significantly  leas  than 
the  rate  of  recombination,  insuring  that  recombination  and 
hence  thermal  release  occurs  at  the  defect  sites  where  the 
bole-electron  pairs  were  generated. 

On  the  basis  of  available  evidence  we  are  unable  to  de¬ 
cide  whether  either  or  both  of  these  suggested  models  are 
valid. 

In  our  results  we  found  that  repeated  lower-energy 
pulses,  i.e.,  0.S  and  1.0  J  cm  2  resulted  in  no  greater  solution 
of  defect  nuclei  than  occurs  with  a  single  pulse.  It  would 
appear  that  the  first  pulse  dissolves  tome  of  the  primary  de¬ 
fect  nuclei  while  converting  others  due  to  released  thermal 
energy  into  secondary  defects.  As  shown  in  our  second  set  of 
experiments,  secondary  defects  are  much  more  resistant  to 
dissolution  by  laser  pulsing,  so  that  repeated  pulsing  contrib¬ 
utes  little  to  defect  removal. 

V.  CONCLUSIONS 

(1)  Pulsed  laser  annealing  of  lower-fluence  ion  implants 
does  not  result  in  foil  electrical  activation.  (2)  Laser  anneal 
•ingle-pulse  energies  of  1.5  J  cm  2  are  sufficient  to  dissolve 
implant  defect  nuclei  in  lower-fluence  ion  implants  without 
causing  melting.  (3)  Following  a  pulsed  laser  anneal  with 
t hemal  annealing  pemits  us  to  obtain  both  dissolution  of  im¬ 
plant  defect  nuclei  as  well  as  complete  electrical  activation 
without  redistribution  of  low-fluence  ion  implanted  impuri¬ 
ties.  (4)  A  two-step  thermal  anneal  followed  by  a  pulsed  laser 
anneal  permits  us  to  obtain  both  dissolution  of  implant  de¬ 
fect  nuclei,  as  well  as  complete  electrical  activation  without 
impurity  redistribution,  for  implant  fluences  sufficiently 
high  to  cause  amorphism  at  the  surface. 

This  work  was  supported  in  part  by  the  Defense  Ad¬ 
vanced  Research  Projects  Agency. 

APPENOIX:  CALCULATION  OF  EXPECTED  DROP  IN 
SHEET  RESISTANCE  RESULTING  PROM  IMPURITY 
REDISTRIBUTION  DUE  TO  THERMAL  ANNEAL 

The  mean  deviation  ARp  for  the  50-keV  boron  implants 
is  3.7  x  IfrVm.  After  a  thermal  anneal,  the  mean  deviation 
is  increased  to  ARp'  which  can  be  expressed  as 

ARp'  -  [4 Rp2  +  20r),',t 

where  D  is  the  diffusion  coefficient  and  t  is  the  diffusion  time. 

For  the  thermal  diffusion  of  implanted  boron,  Ricco  « 
al. 11  found  that  the  diffusion  coefficients  of  Kendall  and 
Dekries14  gave  a  good  fit.  Using  their  value  for  1000  *C,  we 
obtain  a  value  of  ARp'  for  a  10  min  diffusion  of  7.38x  10 2 
pm. 
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Nicholas1'  has  published  an  empirical  expression  giv¬ 
ing  the  sheet  resistance  An  as  a  function  of  the  fluencefand 
(he  mean  deviation, 

*u  =  1 X  10"(  F )  '"(AKpX  ia4)  ° 

The  percentage  decrease  in  sheet  resistance  due  to  thermal 
anneal  can  therefore  be  calculated  as  7.5%. 
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1.0  INTRODUCTION 


The  purpose  of  this  report  Is  to  establish  whether  free  carrier 
contribution  to  optical  absorption  In  semiconductors  will  significantly 
modify  absorption  characteristics  of  Intrinsic  semiconductors.  The  work 
which  Is  described  in  this  report  has  been  stimulated  by  an  Idea  of  using 
lasers  for  processing  of  surface  regions  of  semiconductor  materials. 

Laser-Induced  surface  processing  (LISP)  Is  based  on  the  notion  that, 
when  Irradiated  by  a  laser  whose  radiation  Is  strongly  absorbed  In  the 
processed  material,  the  surface  of  the  irradiated  material  Is  elevated  to 
high  temperatures  In  controllable  manners.  In  the  usual  mode  of  operation, 
LISP  uses  a  high-power  laser  whose  frequency  Is  higher  than  the  semi¬ 
conductor  bad  gap  frequency,  thus  causing  In  It  band-to-band  transitions. 
High  densities  of  free  carriers  may  result  due  to  photoexcitation  of 
electrons  into  the  conduction  band.  One  then  asks  whether  these  high 
densities  will  contribute  significantly  to  the  optical  properties  of  the 
absorbing  material,  specifically,  to  Its  absorption  characteristics. 
Significant  changes  In  the  optical  properties  of  the  semiconductor  would 
In  turn  affect  absorption  of  the  laser  beam,  l.e.,  the  process  through 
which  the  carriers  have  been  created. 

It  Is  well  known  that  high  concentrations  of  free  carriers  In  metals 
may  cause  significant  absorption  In  visible  or  Infrared,  due  to  collective 
motion  of  the  carriers.  The  absorption  peaks  up  just  below  the  plasma 
frequency.  In  semiconductors,  high  carrier  concentrations  may  be  obtained 
by  heavy  doping  of  the  semiconductors,  thermal  excitation  at  elevated 
temperatures,  carrier  photo-generation  and  current  Injection,  among  other 
things.  The  relative  Importance  of  the  contribution  of  the  free  carriers 
to  optical  absorption  will  depend  on  the  radiation  frequency.  In  the 
frequency  region  below  the  band  gap  energy,  free  carrier  absorption  may  be 
the  major  mechanism  of  optical  absorption.  In  the  frequency  region  above 
the  band  gap  energy,  band-to-band  transitions  may  be  expected  to  be 
stronger  than  collective  excitations  unless  the  carrier  concentration 
is  very  high.  Thus,  the  strength  of  the  absorption  Is  expected  to  depend 
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on  the  number  of  free  carriers,  carrier  lifetimes  and  frequency.  It  Is  the 
carrier  contribution  to  optical  absorption  In  the  frequency  region  above 
the  band  gap  frequency  which  is  of  Interest  In  LISP. 

The  Increase  In  the  number  of  free  carriers  In  strongly  absorbing 
semiconductors  during  LISP  Is  due  to  (1)  carrier  photogeneration  and 
(2)  thermal  excitation  of  a  large  number  of  free  carriers.  As  has  already 
been  noted,  the  primary  mechanisms  for  optical  absorption  In  strongly 
absorbing  semiconductors  are  band-to-band  transitions,  creating  (for  long 
pulse  length  and  long  carrier  lifetimes)  large  Instantaneous 
concentrations  of  photoelectrons.  The  source  of  the  thermally  excited 
electrons  Is  heat;  thermal  excitations  occur  after  carriers  photo-excited 
Into  the  conduction  band  relax  back  Into  the  valence  band,  transferring 
optical  energy  to  thermal  vibrations  of  the  lattice. 

It  may  be  of  Interest  to  note  here  that  thermal  runaway  In 
semiconductor  laser  windows  Is  caused  by  thermal  excitations.  These 
windows  are  fabricated  from  Intrinsic  semiconductors  transparent  to  the 
laser  frequency.  However,  for  high  power  lasers,  even  a  very  low  residual 
absorption  causes  the  window  to  heat  up,  thus  thermally  exciting 
carriers  Into  the  conduction  band  and  rapidly  Increasing  absorption  with 
temperature. 

This  report  will  be  organized  as  follows.  The  number  of  photo- 
excited  and  thermally-generated  carriers  during  LISP  will  be  calculated 
In  Section  2.0.  The  diffusion  equation  will  be  used  for  estimating  values 
of  instantaneous  densities  of  the  photogenerated  carriers.  A  theory  of 
optical  absorption  through  collective  motion  of  free  electrons  Is  developed 
In  Section  3.0.  In  Section  4.0  we  calculate  the  absorption  coefficient 
due  to  the  free  electrons  as  a  function  of  frequency  and  further  discuss 
our  results. 
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2.0  CARRIER  CONCENTRATION 


2.1  Transport  of  Photogenerated  Carriers.  The  Instantaneous 
concentration  of  photogenerated  carriers  n(x,t)  In  a  semiconductor  In  the 
absence  of  an  accelerating  electric  field  Is  obtained  In  the  low  carrier 
concentration  approximation  by  solving  the  following  diffusion  equation. 


(1) 


where  the  generation  function  g  gives  the  photogeneration  rate  of  the 
carriers, 

9  -  a  Fe_0lX ,  (2) 

where  a  Is  optical  absorption  coefficient  and  F  =  F(t)  gives  the  rate  of 
the  photon  flux  Incident  at  a  unit  surface  of  the  semiconductor. 

An  =  n(x,t)  -  nQ,  (3) 

where  nQ  Is  the  equilibrium  carrier  concentration,  t  and  D  are  minority 
carrier  lifetimes  and  diffusion  constant,  respectively. 


Assuming  that  g,  x  and  D  are  Indenpendent  of  carrier  concentration, 
Eq.  (1)  can  be  solved  by  using  Green  function  techniques.  For  the 
boundary  conditions. 


j(x  *  0,t)  s  -  o  |5 


=  0 


x=0 


(4) 


and 


n(x  -  t)  *  nQ  , 


(5) 


the  Green  function  Is 


G(x,x' ,t-t' ) 


1 

/4irD(t-t'  ) 


(6) 


♦"  / 

st 

n«-°° 


(x-x1  -  2nt)2  (x+x'-2n&)2  \  (t-t1 ) 

40(t-t' )  +e'  4D(t-t' )  /’  x 
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where  i  Is  the  semiconductor  thickness.  The  boundary  conditions  (4)  and  (5) 
describe  the  physical  situation  when  there  Is  no  flux  of  carriers  across 
the  semiconductor  front  surface,  and  carrier  concentration  at  the  slab  back 
boundary  Is  constant  (nQ  is  the  equilibrium  minority  carrier  concentration). 

The  solution  of  Eq.  (1)  is  then  given  by 

»  ^  t 

n(x,t)  *  f dx  J dt'  G(x,x',  t-t')  g(x',t')  ,  (6) 

0  0 


where  the  generation  function  g(x,t)  is  given  by  Eq.  (2). 

For  the  slab  thickness  i  such  that 

l  »  1/a  ,  (7) 


the  semiconductor  can  be  replaced  by  a  semi -Infinite  medium,  with  the  result 
that  the  summation  I  In  Eq.  (6)  Is  reduced  to  one  term,  namely,  term  n=0. 
Carrier  concentration  is  then  equal  to 


oo  t 


n(x,t)  *  J dlt'  G(x,x' ,  t-t')  oF(t' 


-  aX 


)l 


(8) 


0  0 


where  the  Green  function  is 


G(x,x' ,t-t' )  * 


1 


-4^4 j(5) 


ArrD(t-t’) 


(■ 


+  e 


Our  reason  for  solving  the  transport  equation  Is  to  find  out  whether 
the  Instantaneous  concentration  of  photoexclted  carriers  Is  sufficiently 
large  so  that  the  strength  of  free  carrier- Induced  optical  absorption 
becomes  comparable  with  band-to-band  absorption.  In  order  to  estimate  the 
strength  of  the  effect,  we  have  evaluated  the  concentration  of  carriers 
photoexclted  at  the  semiconductor  front  surface  for  a  square  shape  laser 
pulse  (See  Figure  1).  Carrier  concentration  at  x  *  0  gives  an  upper 
estimate  for  concentrations  Inside  the  semiconductor  (x  >  0). 
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FIGURE  1.  Photogeneration  and 


For  the  square-shape  laser  pulse,  the  generation  function  F(t)  is 
equal  to 

F(t)  =  Fq  0  <  t  <  t  ,  (10) 

F(t)  =0  tp  <  t  <  - 

where  t  is  the  length  of  the  laser  pulse  and  FQ  is  a  (constant)  number 
of  arriving  photons  at  the  semiconductor  front  surface.  Using  Eqs.  (9) 
and  (10),  Eq.  (8)  then  reduces  to 

n(x=0,t)  -  J" dt'  e_t  ^7  "  ^  erf  c  (a/oF  ),  (11) 

t_tP 

where  erf  c  (z)  is  the  complementary  error  function. 

Clearly,  the  number  of  photogenerated  carriers  is  largest  at  the  end 
of  the  laser  pulse.  Since,  for  our  purposes,  we  are  interested  in  the 
maximum  values  of  the  concentration  of  the  photo-excited  carriers,  in  what 
follows,  we  shall  be  concerned  with  the  integral, 

S  i  2 

n(0,tp)  -  aFo  J  dt  e"t(x  ’  *  erf  c  (a/Dt~).  (12) 

0 

Inspecting  Eq.  (12),  we  find  that  there  exist  various  physical  regimes  of 
our  Interest  characterized  by  relative  length  of  three  time  constants: 
laser  pulse  length  (tp),  carrier  lifetime  (t)  and  "diffusion"  time  (1/Da2). 

In  order  to  obtain  approximate  expressions  for  the  concentrations,  in  what 

follows,  we  shall  therefore  evaluate  Eq.  (12)  in  the  following  limits: 
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2.1.1  Steady  State  or  a  very  long  pulse  limit:  tp  >>  t,  t  >>  1/Da2. 
Setting  the  Integral  upper  boundary  to  Infinity,  t  -*•  »,  Eq.  (10) 

yields 

n^°’V  =  T+aL  ’ 

where  the  diffusion  length  L  =  /Dt  .  That  Is,  photocarrier  concentration 
Is  limited  by  the  length  of  the  carrier  lifetime  or  the  length  of  the 
diffusion  time,  whichever  Is  smaller.  Here,  the  "diffusion  time"  is 
defined  as  the  length  of  time  It  takes  the  carrier  to  diffuse  away  from 
the  surface  region  (of  the  depth  1/a),  where  the  carrier  was  originally 
created.  Into  the  semiconductor  interior. 


2.1.2  Short  pulse. 

2. 1.2.1  Pulse  shorter  than  diffusion  time:  t_  <  1/Da2 

-  p 

Expanding  the  error  function  for  small  arguments,  Eq.  (12)  yields 


n(0,tp) 


aFo 


-t„(l  -  Ott2) 
e  p  T  -1 

D»2  -  I 


04) 


T 

In  order  to  better  understand  the  physical  content  of  this  result,  we 
further  reduce  Eq.  (14)  in  the  following  limits: 


(a)  t  <  tp(tp  <  l/g2P). 


n(0,tp)  -  «F0t 


(15) 


For  the  carrier  lifetime  shortest  of  the  three  time  constant, 
carrier  concentration  is  carrier  lifetime-limited. 


(b)  tp  <  t  (t  <  l/g2P  or  t  <  1/g2P  >  tp). 

n(0,tp)  =  aF0  tp.  (16) 

For  the  pulse  length  shorter  than  the  carrier  lifetime,  the 
concentration  is  pulse  length-limited. 
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2. 1.2. 2  Pulse  longer  than  diffusion  time:  t  >  1/Da2 

In  this  limit,  the  integral  in  Eq.  (12)  is  broken  into  two  intervals, 
namely,  t  <  1/Da2  and  t  >  1/Da2.  In  the  first  Interval,  the  error  function 
is  expanded  for  small  arguments.  In  the  second  Interval,  the  error 
function  is  expanded  for  large  arguments.  We  obtain 


A  further  Insight  Into  Eq.  (17)  Is  obtained  by  reducing  It  into  simpler 
expressions  for  the  following  limiting  cases: 

(a)  tp  <  t  (tp  >  1/Da2). 


Equation  (17)  reduce  to 

n(0,tp>  .  F0 


08) 


The  result,  Eq. 
the  following  form. 


(8),  can  be  better  understood  by  rewriting  It  into 


(18' ) 


Notice  that  denominator  of  Eq.  (18')  has  dimension  of  velocity.  That  Is, 
for  tp  <  t,  carrier  concentration  Is  limited  by  carrier  diffusion  within 
the  time  Interval  equal  to  the  length  of  the  laser  pulse. 
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(b)  1/Da2  <  t  <  tp. 

Equation  (17)  reduces  to 

»<«-V  ■  f0  >TF  •  091 

which  can  be  rewritten  as 

n(0,t  )  =  (19* ) 

p  /Dt 

T 

That  Is,  by  analogy  with  Eq.  (18*),  for  t  <  tp,  the  concentration  Is 
limited  by  carrier  diffusion  within  the  time  Interval  equal  to  the  carrier 
lifetime. 

(c)  t  <  1/Da2  <  tp. 

Equation  (17)  reduces  to 

n(0,tp)  «  aFQT  (20) 

That  Is,  carrier  concentration  Is  limited  by  carrier  lifetime.  (It  should 
be  noted  here  that  Eqs.  (19)  and  (20)  In  fact  reproduced  the  steady  state 
result,  viz. ,  Eq.  (13)). 

The  results,  Eqs.  (13)  through  (20),  are  summarized  in  Table  1.  It  Is 
seen  that  the  following  "rule-of-thumb"  applies.  The  surface  concentration 
of  the  photoexclted  carriers  Is  limited  by  the  length  of  the  shortest 
time  Interval  In  the  problem,  t  or  t.  If  these  are  the  shortest  time 
Intervals;  or.  In  the  case  of  fast  carrier  diffusion,  by  diffusion  of 
carriers  within  the  time  Interval  t  or  t,  whichever  Is  shorter.  Comments 
In  the  third  column  of  Table  1  refer  to  real  physical  cases  to  which  the 
respective  limiting  expressions  apply.  The  values  of  the  representative 
parameters  (l.e.,  time  constants,  mobilities,  dlffuslvltles)  used  for 
Identifying  the  various  physical  regimes  are  also  listed  In  Table  1.  In 
the  last  column  of  Table  1  are  listed  plasma  frequencies  corresponding  to 
the  maximum  concentrations  (column  four).  The  plasma  frequencies  were 
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obtained  by  evaluating  the  expression 


where  c  Is  the  dielectric  constant  and  nr  Is  electron  effective  mass.  The 

e 

maximum  concentrations  were  evaluated  for  the  photon  flux 
25  -2  -1 

FQ  •  3  x  10  cms  ,  where  FQ  ■  PQ/hv.  PQ  Is  the  peak  power  of  the  ruby 
laser  (x  ■  .6943  pm)  and  was  taken  PQ  =  107  w/cm2.  The  pulse  width  was 
30  nanoseconds.  The  pulse  was  square-shaped. 

2.2  Surface  Recombination  Velocity 

The  effect  of  the  presence  of  the  surface  In  the  vicinity  of  the  active 
region  Is  to  Introduce  additional  carrier-loss  mechanism  Into  the  problem. 

It  Is  customary  to  define  an  effective  surface  recombination  velocity  s  of 
minority  carriers  by  the  expression: 


0 


an 

ax 


x=0 


s  An(0) 


s(n(0)  -  nQ), 


(22) 


where  nQ  Is  the  equilibrium  concentration  of  minority  carriers.  Thus,  If 
one  wishes  to  account  for  the  effect  of  the  surface  recombination  on 
carrier  concentration,  Eq.  (22)  will  replace  the  boundary  condition, 

Eq.  (4),  of  the  previous  section. 


The  transport  equation,  Eq  (1),  with  the  boundary  conditions,  Equations 
(5)  and  (22)  may  be  solved  using  again  the  Green  function  techniques.  We 
have  not  constructed  the  complete  Green  function  for  this  problem.  However, 
we  have  solved  Eq.  (1)  directly  for  the  case  of  the  steady  state.  In  the 
steady  state,  solution  to  Eq.  (1)  with  the  boundary  conditions  Eq.  (5) 
and  (22)  and  with  the  source  term  Eq.  (2)  applied  over  a  long  period  of 
time  Is 


An  » 


a  Ft 

1  -  a2 l2 


(  (s  +  aO) 
'  S  +  L/t 


■x/L 


(23) 


w*»ere  all  parameters  were  defined  in  Section  2.1. 


Taking  the  limit  of  large  s,  s  >  aD  and  s  >  L/t,  Eq.  (23)  reduces  at 
x  ■  0  to  the  simple  expression: 


An{0)  *  sTTtVlcO 


Thus,  In  the  case  of  a  rapid  surface  recombination,  concentration  of 
photoexclted  carriers  Is  limited  by  surface  recombination  velocity,  s,  as 
expected.  Intuitively,  one  expects  a  similar  result  also  for  the 
transient  case:  For  a  finite  pulse  length,  carrier  concentration  will  be 
surface  recombination-limited  If  the  surface  velocity  Is  higher  than  any  of 
the  other  characteristic  velocities  In  the  problem,  l.e..  If 


s 


>( 


1  L  1 

ot’  t*  at  * 
P 


aD) 


(25) 


In  the  other  cases,  l.e.,  for  the  surface  recombination  slower  than 
any  of  the  characteristic  velocities,  carrier  concentrations  will  be 
limited  by  the  corresponding  velocities,  viz.,  Eqs.  (15)  through  (20). 


2.3  Thermal  Excitations 

Heating  the  semiconductor  to  high  temperatures  will  have  a  threefold 
effect  on  the  smalconductor  properties:  (1)  Increase  the  Intrinsic 
carrier  densities,  (2)  Increase  of  extrinsic  carrier  densities,  and 
(3)  temperature  dependence  of  carrier  lifetimes,  mobilities  and 
dlffuslvltles. 


2.3.1  Intrinsic  Carrier  Densities.  Intrinsic  carrier  densities  of 
Ge,  SI  and  6aAs  as  a  function  of  reciprocal  temperature  are  shown  In 
Figure  2. 
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Figure  2. 


Intrinsic  carrier  densities  of  Ge.  si  >  and  GaAs  as  a  function  of 
reciprocal  temperature. 
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2.3.2  Extrinsic  Carrier  Densities.  In  the  limit  of  low  temperature 
(kT  <  Eg  -  E^)  or  a  high  doping  concentration,  the  ratio  of  Ionized 
dopants  at  temperatures  and  T2  Is 


n(T2) 


exp 


E  -E, 
g  d 

/L. 

_L\ 

2k 

VT2 

V 

* 

(26) 

where  -  E^  Is  the  position  of  the  dopant  level  below  the  bottom  of  the 
conduction  band.  In  the  limit  of  high  temperature,  or  low  doping,  the  Ionized 
dopant  concentration  Is  approximately  constant,  equal  to  the  total  dopant 
concentration. 


n  »  N 

Thus,  at  elevated  temperatures,  the  maximum  concentration  of  extrinsic 
electrons  Is  equal  to  the  dopant  concentration,  yielding  the  following 


plasma  frequencies 

Up  ■  (4tt  ne2/em)^2: 

GaAs 

up  *  6.2  x  104  4T 

(27a) 

S± 

a  *  2  x  io4 

p 

(27b) 

2.3.3  Temperature  Dependence  of  Physical  Constants.  Temperature 
dependence  of  semiconductor  parameters  such  as  lifetimes,  mobilities  and 
dlffuslvltles  Is  rather  complex.  In  the  context  of  the  present 
Investigation,  changes  In  temperature  will  modify  Intervals  of  physical 
regimes  defined  In  Table  I  of  Section  2.2.  Due  to  their  complexity,  and 
their  comparatively  small  effect,  we  shall  neglect  these  temperature 
changes  In  the  present  problem. 
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3.0  EFFECTS  OF  FREE  ELECTRONS  ON  SEMICONDUCTOR  ABSORPTION 

3.1  Propagation  of  electromagnetic  wave  In  a  lossy  medium. 

Attenuation  of  electromagnetic  waves  In  semiconductors  follows  from 
solving  Maxwell  Equations  In  a  conducting  medium.  Solving  the  Maxwell  Equa¬ 
tions  for  the  electromagnetic  wave. 

E  *  E0  exp  [1(kfx  -  ut)]  ,  (28) 

propagating  In  a  lossy  medium,  one  obtains  for  the  complex  (frequency-dependant) 
wave  vector  k, 

E*  *  T  %  ,  (29) 

c 

where  u  is  the  angular  frequency  and  the  generalized  complex  dielectric  constant 
t  Is 

e  *  «  +  1  ~  •  (30) 

Alternatively,  since  the  generalized  complex  conductivity  Is 


fe  (1  -  e) 


Equation  (31)  can  be  rewritten  as 

7=l+i$I_o  _  (32) 

In  Eqs.  (29)  -  (32),  e.o  and  e,o  are  complex  and  real  quantities,  respec¬ 
tively,  and  are  functions  of  angular  frequency  ».  Using  Eq.  30,  the  gener¬ 
alized  refractive  Index  n  , 

n  =  /l  ■  n  +  Ik  ,  (33) 


can  be  written  as 


.2  2  . 

n  -  k  ■  c 


K"  ZT  "  0  • 
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where  n  Is  the  real  part  of  the  generalized  refractive  Index  n  and  <  (the 
Imaginary  part  of  rf)  Is  related  to  the  extinction  coefficient  a  (extinction 
of  energy  carried  by  the  electromagnetic  wave),  by 

a  «  2ie  |  .  (36) 

Hence,  the  absorption  coefficient  6  (describing  attenuation  of  the  amplitude 
of  the  electromagnetic  wave)  Is  equal  to 

=  (37) 

and  the  electromagnetic  wave  Is 

E  =  E  e1 (kx-wt)  -  .  (38) 

o 

We  note  that  Eqs.  (34)  and  (35)  can  be  rewritten  alternatively  as 


giving  n  and  k  directly  In  terms  of  e  and  a.  Once  the  dielectric  constant  and 
conductivity  are  known,  the  above  Equations  completely  describe  dispersion 
and  attenuation  of  electromagnetic  waves  propagating  In  a  lossy  medium. 

3.2  Electron  Susceptibility  and  Conductivity 

In  a  radiation  field,  the  current  of  the  free  carriers  has  a  component  In 
phase  and  a  component  out  of  phase  with  the  electric  field.  The  In-phase 
current  contributes  to  the  conductivity  and  the  out-of-phase  component  con¬ 
tributes  to  the  susceptibility.  In  order  to  estimate  the  free  carrier  con¬ 
tribution  to  o  and  e.  It  Is  convenient  to  employ  a  model  of  a  damped 
oscillator.  In  this  model,  carrier  motion  Is  described  by  the  following 
equation  of  motion. 


(41) 


1  »  ? 

m(x  +  t  x  +  ui_x)  =  eE^  e 
o  o 

where  the  damping  constant  i  is  electron  “relaxation"  time  and  uQ  is  the 
restoring  force.  Since,  for  free  carriers,  there  is  no  local  restriction  on 
their  motion.  In  what  follows,  the  restoring  force  will  be  neglected. 
Putting  u>Q  *  0,  and  using  the  definitions  polarization  P  =  nex  and  current 
j  ■  nex,  we  obtain 


e 


eo  +  4tf 


U>(  1  -itjJT  ) 


) 


(42) 


a 


ne2t 


l(l  -  Iut) 


(43) 


2  4itne^ 

where  u  Is  the  plasma  frequency,  u  *  (-— — ),  and  n  Is  the  free  electron 

p  p  tin 

density.  It  Is  straightforward  to  verify  that  Eqs.  (42)  and  (43)  satisfy 
the  relations  of  Section  III .1  -  Rewriting  Eqs.  (42)  and  (43)  as 

2  2 


and 


,  “p  T  * 
e  =(1  -  22  )Eo 

1  +  U)  T 


ne2t _ 

0  m(l  +  u^t*) 


(44) 


(45) 


and  substituting  them  Into  Eqs.  (39)  and  (40,  the  absorption  coefficient  can  now 
be  evaluated. 

The  final  remark  will  concern  time  constant  t.  t  Is  the  relaxation  time, 
which.  In  the  context  of  the  damped  oscillator  model,  is  related  to  carrier 
mobility. 


(46) 


where  m  and  e  are  electron  effective  mass  and  charge,  respectively. 
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4.0  ABSORPTION  COEFFICIENT 


The  refraction  and  absorption  coefficients  associated  with  free 
electrons  are  calculated  by  substituting  Eqs.  (44)  and  (45)  into  Eqs.  (39) 
and  (40).  The  results  are  sketched  in  Figure  3  through  5. 

In  Figure  3,  we  plot  the  qualitative  behavior  of  both  the  real  and 
imaginary  parts  of  the  generalized  refractive  index  versus  frequency  for 
the  case  of  t  ■  •.  This  is  a  rather  illustrative  case.  It  shows  that,  in 
the  absence  of  any  loss  mechanism,  there  is  a  total  reflection  at  frequenci 
u)  <  Up.  At  u  <  Up,  the  wave  is  evanescent  with  the  skin  depth  6(6  h  b"1 ) 
approaching  infinity  as  frequency  u  approaches  up.  In  the  frequency 
region  above  the  plasma  frequency ,u  >up,  there  is  zero  absorption 
(k  =  0),  the  wave  being  partially  reflected  and  partially  transmitted. 

There  is  a  practically  complete  transmission  at  very  high  frequencies. 


In  Figure  4,  we  sketch  the  qualitative  behavior  of  real  and  Imaginary 


parts  of  the  generalized  refractive  index,  n  and  *,  respectively,  versus 
frequency  for  both,  the  case  Up>  t-1  and  wp  <  t'1 .  Also  noted  in  Figure  4 
are  leading  terms  of  expansions  of  n  and  ic  versus  u  in  various  frequency 


regimes;  that  Is,  regimes  defined  by  relative  magnitudes  of  parameters  u. 

Up  and  t  .  The  leading  terms  were  obtained  by  expanding  expressions  for 
n  and  k  in  terms  of  smallness  parameters  appropriate  to  the  various  regimes 


It  Is  seen  that  finite  values  of  n  and  <  are  obtained  at  all  frequencies 
Indicating  partial  transmission,  reflection  and  absorption  of  electro¬ 
magnetic  wave  in  all  the  frequency  regions. 


In  Figure  5,  we  sketch  qualitatively  the  Imaginary  part  of  the  refractive 
Index,  k,  and  the  absorption  coefficient  B;  Figure  5.1  shows  the  case  of 
Up  <  t"1  and  Figure  4.2  the  case  up  >  t"^.  By  comparing  Figure  5.1  and 
5.2,  It  Is  seen  that  the  qualitative  behavior  of  b  in  the  two  cases  is  the 
same.  The  absorption  coefficient  starts  from  zero  at  u  *  0  and  grows 
approximately  with  frequency  as  (u)^2  until  u  approaches  frequency 
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(GOODMETALS:  r",«Wp/  Xp  3  1.3pm) 


FIGURE  3.  Real  (n)  and  Imaginary  (k)  parts  of  the  generalized  refractive 
Index  In  the  absence  of  any  absorption  mechanism.  (Good  metal 
approximation:  t-1  «  w  a  *  1.3  ym) 
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(o)  (INTRINSIC  SEMICONDUCTOR) 


FIGURE  5. 


lalltatlve  behavior  of  the  Imaginary  part  of  the  Index  of 
•fraction  (<)  and  the  optical  absorption  coefficient  ($)  due 
>  the  free  electron  effect,  (a)  t*1  >  up  (Intrinsic  semi- 
jnductors ) ;  (b)  up  >  t'1  (extrinsic  semiconductors). 
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(b)  (INTRINSIC  SEMICONDUCTOR) 

FIGURE  4.  Qualitative  plot  of  the  real  (n)  and  Imaginary  (k)  parts  of 
the  generalized  refractive  Index,  (a)  uD  >  t*1  (extrinsic 
semiconductor);  (b)  uD  <  t-1  (Intrinsic  semiconductor). 


«)  ■  a)  or  t  ,  whichever  is  smaller.  Absorption  then  levels  off  (attains 
P  _1 

maximum  value)  until  the  frequency  reaches  t  or  wp,  whichever  is  larger. 

Thereafter,  e  decays  Inversely  proportional  to  the  second  power  of  the 

frequency.8  ^  of2. 

Thus,  we  can  see  that  there  exist  several  frequency  regions  in  which 

the  absorption  coefficient  behaves  differently.  It  is  easy  to  verify  that, 

for  semiconductors,  unless  the  semiconductor  is  intrinsic,  the  limit  of 

t”1  <  a)  usually  applies.  Using  Eq.  (46),  one  finds  that  typical  values 
P  _-i  3  _12 

of  t  range  from  4x10  to4x10  seconds  (SI,  InAs,  GaAs  and  Ge); 

-1  14  -3 

therefore,  t  <  <op  for  electron  concentrations  above,  typically,  10  cm  . 

Hence,  for  most  semiconductors,  the  absorption  coefficient  8  *  w„/c  In 

-1  2  2** 
the  frequency  region  of  t  <  «d  <  up  and  8  =  iup  /2ctu>  for  frequencies 

u  >  (Up.  A  list  of  calculated  values  of  8  at  the  ruby  frequency 

(x  *  .694  pm)  and  at  7  pm  for  silicon,  InAs,  GaAs  and  germanium  for 

17  18  19  -3 

doping  concentrations  of  10  ,10  and  10  cm  Is  given  In  Table  2. 

In  Table  2,  the  effect  of  frequency  and  carrier  concentration  on  values 

of  8  Is  clearly  seen.  Absorption  peaks  up  In  the  vicinity  of  u>  *  Wp. 

With  frequencies  Increasing  Into  the  near  Infrared  region  of  the 

spectrum,  the  absorption  decreases  until.  In  the  visible  region  of  the 

spectrum.  It  become  very  small.  We  note  that  this  region  also 

coincides  with  the  region  of  the  onset  of  strong  band-to-band  transitions , 

hence.  In  semiconductors,  the  effect  of  free  carriers  on  the  optical 

properties  becomes  important  really  only  at  frequencies  lower  than  the 

bandgap  frequency.  Whereas  at  frequencies  u  <  un,  the  absorption 

4-1  “  10  -3 

coefficient  8  *  w./c  »  10  cm  for  electron  concentration  of  n  *  10  cm  , 
P 

the  value  of  the  absorption  coefficient  at  the  band  gap  frequency  Is  lower 
than  this  maximum  value  by  several  orders  of  magnitude. 

We  note,  finally,  that  in  the  near  Infrared  region,  where  the 
frequency  u  >  up.  It  holds  that 
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TABLE  2,  ABSORPTION  COEFFICIENT  ASSOCIATED  WITH  FREE  ELECTRONS. 


16ir 2a2/w2  <  e 


(47) 


and  hence. 


n  a  ST 


(48) 


In  addition.  If  the  carrier  concentration  is  not  too  high,  the  dielectric 

constant  e  may  be  approximated  by  Its  value  In  the  absence  of  free 

2 

electrons,  e  “  e0  -  n0  •  Therefore,  electrons  affect  optical  properties 
only  through  their  contribution  to  conductivity.  Influencing  k  or  b,  but 
not  n:  „ 


6 


n„fa)_ 

fa)  CTT  0  OP 

—  |T  S  -  ■  -  »  ■  *-■  ■ 

c  n_  c  -  2 

O  2c  fa)zT 


(49) 


which  Is  the  expression  written  In  Figure  4  for  the  limit  of 

fa)  >  (r  ,  ft>p). 


In  the  usual  mode  of  Its  operation,  LISP  utilizes  lasers  whose  radiation 

Is  strongly  absorbed  In  the  processed  material,  necessitating  frequency 

regimes  above  the  band  gap  frequencies.  The  ruby  laser  Is  used  for 

Irradiation  of  semiconductor  surfaces  In  our  laboratory.  The  Instantaneous 

concentrations  of  free  electrons  Induced  at  the  semiconductor  surface  by 

the  laser  Irradiation  were  calculated  In  Section  1.0  to  be  on  the  order  or 
19  -3 

less  than  10  cm  .  Hence,  as  has  already  been  noted.  In  the  frequency 
regime  above  the  band  gap  frequency  where  the  semiconductor  Is  strongly 
absorbing  (a  >  10*  cm~^),  optical  absorption  due  to  the  photoexclted 
electrons  Is  less  than  the  band-to~band  absorption,  and  hence,  the 
effect  of  free  electrons  on  the  absorption  coefficient  can  be  neglected. 
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SUMMARY  OF  FINAL  REPORT,  PART  II 


-V  U 

The  objective  of  Tertc  2  was  to  evaluate  the  laser  blow-off  Ion 
Implantation  source.  This  part  reports  the  results  of  this  evaluation. 


1.0  SUMMARY 


A  new,  potentially  less  expensive  and  relatively  simple  technique 
for  Ion  Implantation  of  extended  areas  of  semiconductor  material  was 
experimentally  developed  under  contract  to  DARPA.  ^Thls  new  method  Involves 
a  source  of  Implant  Ions  produced  by  laser  ^blow-off*  from  a  piece  of 
dopant  material.  Energy  from  a  short-pulse  laser  Is  focused  onto  the 
surface  of  the  dopant  material,  forming  a  small  bubble  of  dopant  plasma 
which  'blows-off*.  This  plasma  bubble  Is  allowed  to  expand  Into  vacuum 
until  Its  density  and  area  become  suitable  for  electrostatic  acceleration  to 
the  necessary  Ion  Implant  energy.  Key  features  Include  the  easy  formation 
of  Intense  Ion  beams  of  vlrtuallyiany  dopant  material,  with  beam  areas 
of  several  Inch  diameter,  and  wlthouivthe  necessity  for  slower  and  more 
expensive  rasterlng  of  Ion  beams  from  rophlstlcated  Ion  beam  machines. 
Perhaps  most  significant,  although  also  more  speculative,  was  the  potential 
for  self-annealing  of  Implant  damage  during,  the  Implant  process.  The 

possibility  arose  from  the  unique  situation  the  short  pulse  of  high 

\ 

Intensity,  high  energy  Ion  beam,  which  rap1dly\depos1ts  energy  exactly 
In  the  Implant  region  to  momentarily  raise  the  material  temperature.  If 
this  possibility  could  have  been  realized  In  practice,  an  entire  and 
sometimes  troublesome  separate  annealing  process  normally  associated  with 
Ion  Implantation  might  have  been  eliminated.  \ 


Recent  theoretical  studies  during  the  present  contract  Indicate  that 
for  the  pulse  durations  and  magnitudes  likely  to  be  obtained,  only  small 
temperature  rises  are  expected.  Thus  self-annealing  due  to  thermal  effects 
alone  Is  not  likely.  However,  the  very  Interesting  and  potentially 
Important  finding  was  made  that  the  much  higher  current  densities  used  here 
(pulsed)  produce  amorphous  silicon  at  total  doses  far  smaller  than  required 
with  conventional  Implantation.  The  effect  Is  quite  striking,  but  no 
satisfactory  explanation  has  been  given  as  yet.  The  production  of 
amorphlsm  at  low  dose  may  prove  very  useful  In  annealing  some  low-dose 
structures . 

During  the  contract,  an  apparatus  was  built  and  tested  for  Implanting 
Ions  by  this  method.  Implants  of  boron  and  antimony  Ions  Into  silicon 
substrates  were  made  at  energies  up  to  65  keV.  The  upper  energy  was 
limited  by  electrical  breakdown  across  vacuum  feedthroughs  and  electrode 
spacing,  a  problem  which  appears  to  be  solvable. 

The  quality  of  the  implants  was  assayed  with  Ion  Microprobe  Mass 
Analysis  (IMMA)  equipment  and  with  Spreading  Resistance  Probe  techniques 
(SRP).  The  IMMA  measurements  are  Ideal  for  determining  precisely  the 
Implant  Ion  depth-profiles  and  surface  areal  uniformity,  and  provide 
meaningful  comparison  between  the  new  and  the  conventional  methods.  In 
fact,  IMMA  was  more  appropriate  as  a  diagnostic  than  actual  semiconductor 
device  fabrication  tests  during  early  phases,  since  the  success  or 
failure  of  device  fabrication  may  say  little  about  large  and  small  scale 
uniformity  and  overall  depth  profiles  which  were  our  Initial  concern. 

Early  major  emphasis  was  therefore  placed  on  the  IMMA  evaluation 
technique. 
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During  the  later  phase,  major  emphasis  was  shifted  from  IMMA  as 
diagnostic  toward  the  quality  of  actual  semiconductor  materials  and 
devices  fabricated  with  the  laser  blow-off  Implantation  technique.  This 
shift  was  particularly  Important  In  studying  possible  self-annealing  with 
the  new  method  since  IMMA  could  say  little  about  semi-conductor  structure 
and  activation  properties.  In  the  final  analysis  the  quality  of  the 
fabricated  semiconductor  device  operation  Is  clearly  the  key  measure  of 
success.  Spreading  resistance  probe  measurements  were  utilized  to 
determine  resistivity  and  Implant  uniformity.  Implant  uniformity  was  found 
to  be  excellent.  A  groove-and-staln  technique  was  used  to  determine 
junction  depths. 

Careful  analysis  of  IMMA  comparisons  between  depth  profiles  for 
conventional  and  laser  blow-off  Implants  raised  the  possibility  of  a 
potential  problem  with  the  blow-off  method.  Namely,  multiply-charged 
Ions  are  present  In  the  blow-off  plasma,  and  result  In  Implant  energies 
at  multiples  of  the  acceleration  voltage.  This  could  produce  a  "tall" 
on  the  depth  profile  extending  to  considerable  distances.  The  experiments 
to  determine  the  actual  junction  depths  showed  that  for  many  Implant  purposes, 
the  problem  Is  not  serious. 
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2 .  EXPERIMENTAL  APPARATUS 

A  schematic  diagram  and  picture  of  the  Ion  Implantation  system 

developed  here  Is  shown  as  Figure  la  and  lb.  This  system  Is  an  adaptation 

of  the  high  Intensity  metal  Ion  beam  source  described  In  Reference  1. 

The  Ions  to  be  Implanted  are  created  by  pulsed  laser  Irradiation  of  a  solid 

block  of  the  Implant  material.  In  the  present  case  a  Nd:YAG  laser  Is 

Q-swItched  with  a  Pockels  cell,  producing  an  0.1  joule  laser  pulse  with 

less  than  0.01  microsecond  length  at  a  wavelength  of  1.06  micron.  The 

laser  beam  passes  through  steering  prisms  and  vacuum  system  windows  and 

Is  focused  to  a  1  mm  diameter  spot  on  the  Implant  material  by  a  50  cm 

focal  length  lens.  Tests  thus  far  have  utilized  solid  boron  and  antimony 

targets.  This  laser  target  Is  located  In  an  electrical  field-free 

region,  with  the  normal  to  Its  surface  directed  approximately  toward  the 

silicon  wafer  to  be  Implanted.  The  laser  beam  strikes  the  surface  from 

an  angle  of  about  25°  to  the  normal.  The  position  and  angle  of  the  boron 

target  are  adjustable  to  optimize  the  Intensity  and  uniformity  of  the  Ion 

beam  In  the  direction  of  the  silicon  wafer,  since  these  parameters  are 

1  2 

complicated  functions  of  the  laser  beam  energy  and  angle.  ' 

The  result  of  the  rapid  deposition  of  laser  energy  on  the  Implant 
material  surface  Is  to  produce  a  high  temperature  plasma  "bubble"  of 
Implant  material  Ions  and  electrons,  which  "blow-off"  from  the  surface. 

In  the  present  configuration  this  plasma  bubble  expands  freely  until  It 
reaches  the  Ion  acceleration  region.  This  distance  Is  presently  32  cm. 

In  the  first  measurements  no  acceleration  was  utilized,  and  the  silicon 
wafers  and  an  array  of  current  measuring  Faraday  cups  were  placed  at  this 
position.  This  set-up  allowed  the  measurement  of  Implant  Ion  fluence  and 
areal  uniformity.  This  configuration  may  prove  to  be  valuable  In  Its  own 
right  as  a  method  for  depositing  a  thin  layer  of  certain  materials  for 
subsequent  laser  diffusion  processing,  since  atoms  of  essentially  any 
solid  dopant  can  be  deposited  under  vacuum  conditions.  The  Ion 
energies  are  a  few  hundred  electron  volts  In  this  case,  and  thus  barely 
penetrate  the  surface. 
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YAG  LASER  M  STEERING  PRISMS 


Figure  la.  Schematic  Diagram  of  Laser  Blow-Off  Ion  Implantation  System 


In  the  case  where  high  energy  implantation  is  to  be  used,  acceleration 
of  the  implant  ions  is  necessary  as  they  reach  the  acceleration  region. 

In  the  present  configuration,  this  acceleration  occurs  between  the  grids 
shown.  The  first  grid  is  held  at  the  same  potential  as  the  field-free 
plasma  expansion  region.  The  second  is  held  at  a  high  negative  voltage 
equal  to  the  ion  kinetic  energy  in  eV  desired.  The  silicon  wafer  is  held 
at  the  same  potential  as  the  second  grid  at  a  distance  approximately  one 
centimeter  beyond.  This  "drift"  distance  is  intended  to  wash  out  any 
non-uniformity  in  implant  caused  by  grid  shadows  or  focusing  effects. 

The  drift  distance  cannot  be  too  long  because  space  charge  spreading  becomes 
severe  because  the  plasma  electrons  are  rejected  by  the  Ion  acceleration 
region.  The  grids  are  presently  100  mesh  stainless  steel. 

For  inital  tests,  two  types  of  silicon  wafers  were  obtained  from 
Monsanto.  Both  types  are  3  inches  in  diameter  and  0.015  inches  thick,  are 
oriented  100,  and  have  resistivities  of  10  ohm-centimeters.  One  type  is 
N  doped  with  phosphorus  and  the  other  type  is  P  doped  with  boron.  In 
order  to  economize  on  wafer  material,  the  wafers  are  cut  into  cm-size  chips 
for  use  when  the  measurements  permit. 

The  accelerating  voltage  is  supplied  by  a  Universal  Voltronics, 

160  kV  power  supply,  which  charges  an  0.05  microfarad  capacitor  across  the 
accelerating  grids.  In  the  event  of  arc  breakdown,  the  grids  would  have 

to  dissipate  60  joules  of  energy  at  50  kV.  If  an  accele, ating  region 

2  2 
area  of  100  cm  and  an  implant  Ion  density  of  1  micro- coulomb/ cm  /pulse 

Is  considered,  the  accelerating  voltage  drops  from  50  kV  to  48  kV  during  the 

pulse.  A  larger  capacitance  value  would  reduce  the  voltage  change,  but 

the  present  value  is  safer  during  development. 

The  vacuum  system  consists  of  fairly  standard  components  Including 
a  10  Inch  oil  diffusion  pump  with  liquid  nitrogen  cooled  baffle  and 
mechanical  forepump.  This  system  should  be  adequate  for  development 
purposes,  but  final  definitive  quantitative  device  tests  may  necessitate 
an  ultra-high  vacuum  system  to  minimize  surface  contamination. 
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2.1  The  Laser  Blow-Off  Plasma 

In  order  to  examine  the  characteristics  of  the  expanding 
plasma  plume  with  regard  to  ion  energy,  current  and  number  den¬ 
sity,  and  angular  distribution,  an  array  of  5  Faraday  cups  was 
placed  symmetrically  along  an  arc  perpendicular  to  the  normal 
to  the  boron  laser  target  in  a  plane  containing  the  normal.  The 
cups  were  41  cm  from  the  target  and  were  spaced  4.2°  apart  with 
respect  to  the  target,  that  is,  the  5  cups  observed  the  target 
at  angles  of  8.4°,  4.2°  and  0°  each  side  of  the  normal.  The 
distance  between  centers  of  the  outer  cups  was  12  cm.  The  Fara¬ 
day  cups  detected  ions  within  the  expanding  plasma  plume  as  the 
ions  reached  the  cup  collectors.  These  cups  were  designed  and 

tested  for  high  current  operation  at  the  plasma  densities  en- 

? 

countered.  The  effective  detection  area  was  1  cm  .  These  cups 
were  of  planar  construction  so  that  the  ion  transit  time  within 
the  cup  was  negligibly  small  compared  to  the  transit  time  from 
the  laser  target  to  the  cup.  Elements  within  the  cup  were  electri¬ 
cally  biased  to  suppress  secondary  electron  emission.  The  Faraday 
cups  could  be  operated  in  either  the  current  sensing  mode  or  the 
total  charge  collection  mode,  with  the  signals  displayed  directly 
on  an  oscilloscope  in  either  case.  The  present  results  with 

boron  ions  are  consistent  with  earlier  theoretical  and  experi- 
2 

mental  work  done  in  this  laboratory  on  laser  blow-off  plumes. 

This  work  involved  a  model  based  on  the  following  physical  as¬ 
sumptions. 

It  was  assumed  that  the  deposited  laser  energy  vaporizes 
and  ionizes  Nq  particles  from  a  small  region  of  the  source  mate¬ 
rial  surface.  This  hot  plasma  "bubble”  of  NQ  ions  begins  to 
expand  against  the  substrate  and  imparts  momentum  to  it,  thereby 

giving  the  center  of  mass  of  the  N  ions  a  velocity  v  in  the 

o  1  cm 

laboratory  coordinate  frame.  Because  of  symmetry,  v  is  normal 

cm 

to  the  surface.  After  the  center  of  mass  of  the  N  ions  has 

o 

moved  a  distance  from  the  surface,  the  number  density  of  ions 
remains  high  enough  so  that  most  of  the  Nq  ions  equilibrate  at 
a  temperature  T  relative  to  their  center  of  mass,  which  continues 
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to  move  at  vcm.  That  is,  an  observer  moving  with  the  center  of 
mass  would  observe  an  essentially  spherically  symmetric  plasma 
"bubble"  in  equilibrium  at  T,  and  freely  expanding.  As  the  number 
density  continues  to  decrease,  the  ion  velocity  distribution 
corresponding  to  equilibrium  at  temperature  T  becomes  "frozen" 
in  the  center-of-mass  system.  A  laboratory  observer  at  a  larger 
distance  from  the  source  laser  spot  observes  ions  arriving  with 
velocity  essentially  along  a  line  from  source  spot  to  observer, 
but  with  a  plume- like  angular  distribution  due  to  the  center-of- 
mass  velocity  vcm  superimposed  on  the  purely  thermal  velocity 
distribution. 


The  theory  thus  consists  simply  of  calculating  the  number 

density  p  and  current  density  j"  of  ions  to  be  observed  in  the 

laboratory  frame  resulting  from  Nq  ions  expanding  with  spherical 

symmetry  and  having  equilibrium  velocity  distribution  corresponding 

to  T  in  their  center-of-mass  system,  and  with  a  superimposed 

center-of-mass  velocity  v  . 

cm 


Figure  2  shows  the  coordinate  system  for  the  theoretical 
model.  The  present  adaptation  differs  in  that  a  solid  target 
is  struck  by  the  laser  beam  on  the  front  surface.  The  current 
density  of  ions  j  in  ions  per  unit  area  per  unit  time  is  given 

Kr.e.x.n,  .(i|i)(|)3/2(il)1/2^.xp(^).(»nj4 
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gives  the  relation  defining  the  "blow-off"  parameter  X.  Small  X 

implies  small  vcm  and  hence  a  broad  angular  distribution,  while 

large  X  implies  a  narrow  "plume".  The  other  parameters  above 

include  time  t,  ion  mass  m,  Boltzmann's  constant  k,  temperature 

T,  and  time  of  J  t„„  .  Figure  3  shows  these  results  graphi- 

cally  in  terms  of  versus  n»  and  therefore  represents  an 

+max 

expected  picture  of  j  versus  t  on  an  oscilloscope  trace  of  a 
current  sensing  Faraday  cup.  Picture  A  of  Figure  4  shows  an 
actual  trace  of  a  boron  ion  pulse  from  the  present  measurement. 
The  points  on  Figure  3  are  the  same  data  plotted  there.  The 

2 

skewness  away  from  the  theoretical  curves  has  been  noted  before 
for  ions,  while  the  agreement  is  quite  good  for  neutral  atom 
blow-off.  The  X  values  depend  on  laser  energy  density  and 
the  target  materials,  and  are  thus  subject  to  some  control.  The 
general  agreement  with  theory  is  sufficiently  good  to  predict 
currents  and  densities  as  functions  of  time,  distance  and  angles. 
Of  particular  importance  is  the  sensitivity  of  the  angular  dis¬ 
tribution  to  X,  as  seen  next.  The  total  number  of  ions  passing 
across  unit  area  normal  to  r  per  pulse  is  given  by 

f  7  dt  -  -2-j  exp(-3X2/2) 

o  4,rr 

x  |\  ♦  /ir(l/2  +  x2)U  +  erf  x)  ex2J  r 


where 

X  »  /3/2 

(X  cos  6) 

2 

X 

and 

erf  x  ■  /? 

f  exp(-u2)  du 

o 


Figure  5  shows  these  theoretical  angular  distributions  as  functions 
of  X.  Also  included  are  previous  experimental  points  obtained 
for  aluminum  ions  and  present  points  obtained  for  boron  ions. 
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It  is  obvious  that  small  values  of  X  are  necessary  if  one  is  to 
implant  large  areas  with  high  uniformity  without  "rastering" 
the  wafer  between  shots.  The  ion  arrival  time  profiles  may  also 
be  very  important  in  self -annealing  aspects,  since  the  thermal 
pulse  will  reflect  the  arrival  time  pulse.  Figures  4  and  5  from 
Reference  2  show  examples  of  widely  differing  angular  distributions 
depending  on  laser  beam  conditions.  The  presently  obtained  angu¬ 
lar  distributions  are  adequate  for  present  purposes  and  will  be 
optimized  further,  but  it  is  clear  that  variations  in  ion  fluence 
as  a  function  of  angle  put  eventual  limits  on  the  size  of  wafers 
which  can  be  implanted  with  high  uniformity  over  the  surface. 

Measurements  have  been  made  to  determine  the  actual  ion  cur¬ 
rent  and  areal  implant  ion  density  by  combining  Faraday  cup  array 
current  measurements  with  implant  ion  deposition  without  accel¬ 
eration.  The  currents  and  charge  collected  are  displayed  directly 
on  an  oscilloscope  for  the  two  outer  and  the  middle  Faraday  cups 
in  the  array.  The  apertures  of  the  other  two  cups  are  covered 
with  silicon  wafer  chips.  A  number  of  laser  shots  are  then  made 
to  deposit  boron  on  the  silicon  chips,  at  the  same  time  summing 
the  charge  collected  on  either  side  of  the  chips.  Table  1  shows 
the  results  of  one  such  run.  In  this  particular  case,  cup  A  was 
on  the  side  of  the  normal  toward  the  laser  beam,  and  the  results 
show  the  asymmetry  toward  the  laser  beam  which  was  reported 
earlier.^  The  uniformity  of  ±  2%  over  the  3  cm  interval  between 
cups  A  and  M  is  acceptable  for  now,  but  the  -  81  variation  between 
cups  B  and  M  is  not.  This  run  shows  a  particular  strong  depedence 
on  angle,  a  dependence  itself  dependent  on  laser  conditions  as 
noted  above. 
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Table  1.  Boron  on  Silicon  Experimental  Results 


3.  CHARACTERIZATION  OF  SILICON  SURFACES  ION  IMPLANTED  BY  THE  LASER  BLOW- 

OFF  TECHNIQUE 

Conventional  commercial  Ion  Implantation  equipment  utilizes  a  collimated 
beam  which  Is  rastered  vertically  and  horizontally  across  the  silicon  wafer 
surface.  This  results  In  a  very  uniform  distribution  of  the  Implanted 
atoms,  one  of  the  desirable  features  of  Ion  Implantation.  In  addition,  the 
use  of  an  analyzing  magnet  Insures  that  the  only  species  Implanted  Is  a 
single-charged  Ion  of  one  of  the  Isotopes. 

The  laser  blow-off  technique  utilizes  a  series  of  plasma  blow-offs  of 
Implant  material  from  the  source  surface  which  are  then  accelerated  to  the 
desired  energy  before  Implantation.  Here  one  expects  to  Implant  all  Isotopes 
of  the  source  and  to  Implant  both  doubly  as  well  as  singly  charged  Ions. 

The  doubly  charged  Ions  will  be  accelerated  to  a  higher  velocity,  equivalent 
to  twice  the  energy  for  singly  charged  Ions. 

Silicon  wafers  were  Implanted  here  by  a  conventional  commercial  Ion 
Implanter  as  standards  of  comparison  for  evaluation  and  characterization 
of  the  blow-off  Implanted  surfaces.  The  Implantation  of  such  standard 
wafers  Is  uniform  and  reproducible  and  the  Impurity  distribution  has  been 
extensively  studied  and  described. 

To  demonstrate  the  characterization  techniques,  five  wafers  were 
prepared.  Three  of  these  were  Implanted  with  a  conventional  Implanter.  They 
were  2-1nch  n-type,  (100)  silicon  wafers  with  a  resistivity  of  2  ncm  and 
were  Implanted  with  1013,  1014,  and  1015B11  at  75  keV.  The  two  blow-off 
prepared  wafers  were  3  Inches  In  diameter.  One  was  10  ficm  n-type  and  was 
Implanted  with  boron;  the  other  was  10  ncm  p-type  and  was  Implanted  with 
antimony. 


3.1  The  as-implanted  Surfaces 

There  Is  a  limited  amount  of  Information  that  can  be  obtained  from  the 
as-implanted  surface.  The  Implanted  Ions  are  on  non-lattice  positions. 
Thus  they  do  not  have  electric  properties  which  can  be  measured.  Normally 
the  as-implanted  surface,  due  to  the  effect  of  the  defects  on  carrier 
mobility,  exhibits  a  very  high  resistivity.  It  requires  a  heat  treatment 
to  permit  the  Implanted  Ions  to  diffuse  to  lattice  positions. 
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The  defect  structure  Is  also  difficult  to  evaluate  In  the  as-implanted 
surface.  The  use  of  transmission  electron  microscopy  (TEM)  also  requires 
the  application  of  a  thermal  anneal.  This  causes  the  agglomeration  of  the 
point  defects,  l.e.,  vacancies  and  Interstitial  atoms,  so  that  a  resultant 
defect  Is  obtained  that  Is  large  enough  to  be  detected  by  TEM. 

An  Important  factor  characterizing  an  Ion  Implanted  surface  Is  whether 

that  surface  becomes  amorphous  as  a  result  of  the  Implantation.  This  Is 

usually  Indicated  by  the  development  of  a  "milky"  appearance.  Such 

amorphlsm  Is  found  In  conventionally  Implanted  surfaces  subjected  to  high 

fluences  of  heavy  Ions.  In  the  laser  blow-off  implants,  amorphlsm  was 

noted  at  surprisingly  low  fluences  at  Implant  energies  of  20  keV,  boron 

14  2 

Into  silicon.  Amorphlsm  was  noted  at  fluences  as  low  as  5  x  10  B/cm  . 

At  Implant  energies  above  50  keV,  only  antimony  produced  amorphlsm  at 
this  fluence.  This  surprising  result  Is  attributed  to  the  very  high 
Instantaneous  current  densities  occurring  during  laser  blow-off  Implantation. 

3.2  Uniformity  of  Implantation 

The  presence  of  group  III  and  group  V  Implanted  Ions  In  silicon  can 
be  easily  detected  by  utilizing  their  electrical  properties.  As  mentioned 
earlier.  It  Is  necessary  to  anneal  the  Implanted  surfaces  to  cause  the 
Implanted  Ions  to  diffuse  to  lattice  positions  where  they  exhibit 
electrical  properties.  A  two  step  thermal  anneal  was  used  to  activate 
the  Implanted  Ions.  The  first  step  was  a  1  hour  anneal  in  dry  Ng  at  550®C. 
This  treatment  was  specifically  developed  for  Implantations  where  the  sur¬ 
face  becomes  amorphous.  In  such  cases  It  permits  the  slow  solid  state 
epitaxial  regrowth  of  the  amorphous  surface  layer.  The  second  step,  a  10 
minute  anneal  at  1000*C,  Is  chosen  to  give  complete  electrical  activation 
with  a  minimum  of  Impurity  redistribution  by  diffusion. 

With  the  Ions  activated  one  can  measure  the  uniformity  of  Implanta¬ 
tion  by  taking  a  traverse  across  the  surface  with  a  spreading  resistance 
probe  (SRP).  Figure  6a  Illustrates  a  SRP  traverse  taken  across  a  diameter 
from  a  starting  position  near  the  edge.  We  note  that  near  the  edge  the 
resistivity  Is  higher,  corresponding  to  a  significantly  lower  boron 
concentration.  Within  the  Inner  2.5  Inch  circle  the  Implant  distribution 
Is  fairly  constant. 
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Figure  6b  Illustrates  an  SRP  traverse  taken  In  the  center  of  the  wafer 
with  a  step,  ax,  of  only  10  um.  The  aim  here  Is  to  determine  whether  the 
accelerator  grid  affects  the  Implantation  pattern.  It  may  be  noted  that 
there  Is  a  very  slight  undulation  In  the  SRP  traverse  with  a  period  of 
250  um.  This  corresponds  to  the  grid  mesh  of  100  per  Inch. 

3.3  Implant  Fluence 

The  Implant  fluence  Is  determined  from  the  product  of  the  number  of 

pulses,  the  duration  of  each  pulse,  and  the  current  density  measured  by 

a  Faraday  cup.  For  the  first  experimental  B  Implantation,  assuming 

14  -2 

singly  charged  Ions,  this  was  calculated  to  be  5.4  x  10  B  on 

K.  H.  Nicholas  (3)  has  developed  a  useful  algorithm  which  permits  one 
to  predict  the  sheet  resistances  of  boron-implanted  layers  prepared  by 
conventional  Ion  Implanters.  One  can  use  this  algorithm  to  estimate  the 
fluence  obtained  with  the  experimental  Implanter.  The  algorithm  derives 
the  sheet  resistance  R  In  a/G  from  the  expression 

R  3  1  x  1011  D**7  (10"4  ARp)  ~’3 

.2 

where  D  Is  the  fluence,  Bern  ,  and  ARp  Is  the  projected  standard  deviation, 

um.  Choosing  the  value  for  ARp  of  .057  um  from  the  tables  of  Bernard 

Smith,  Harwell,  for  a  50  keV  Implantation  and  using  the  value  of  330.7 

obtained  by  4-pt  probe  measurement,  one  can  calculate  a  fluence  of 
14  -2 

2.3  x  10  Bern  .  By  comparison  a  similar  calculation  for  the  conven¬ 
tionally  Implanted  B  fluences  of  1  x  10^3,  1  x  10^4  and  1  x  10^5B^cm~2 
with  Ro's  of  1995,  493,  and  92.25  n/a  give  values  of  1.64  x  1013, 

1.21  x  1014  and  1.32  x  1015Bcm"2. 

The  fluence  calculated  for  the  experimental  Implantation  Is  less 
than  half  of  that  obtained  on  the  basis  of  the  Faraday  cup.  One 
possibility  Is  that  a  significant  fraction  of  the  Ions  were  multiply 
charged  resulting  In  an  actual  fluence  lower  than  that  calculated  on  the 
basis  of  singly  charged  Ions.  This  becomes  evident  on  the  basis  of 
Implant  penetration  measurements  to  be  described  next. 
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3.4  Implant  Penetration 

The  penetration  of  the  Implanted  ions  Is  determined  by  three  factors, 
the  atomic  or  molecular  weight,  the  charge  number  on  the  Ion,  and  the 
accelerating  voltage.  Because  of  their  closeness  In  atomic  weight.  Iso¬ 
topes  will  exhibit  only  slight  differences  in  penetration.  Multiple  atomic 
Ions,  If  such  occur  will  penetrate  a  significantly  lesser  distance  into 
the  silicon  while  multiply  charged  Ions  will  penetrate  a  significantly 
greater  distance.  An  accurate  model  (LSS)  exists  and  the  distributions 
for  these  conditions  can  be  predicted. 

The  Ion  microprobe  represents  a  means  of  obtaining  an  impurity  dis¬ 
tribution  profile.  Comparing  such  a  profile  for  an  experimental  laser 
blowoff  Implantation  with  one  obtained  for  a  conventional  Implantation 
of  the  same  fluence  and  energy  enables  the  evaluation  of  the 
contributions  of  multiple  atomic  species  and  multiply  charged  ions. 

Implanted  Impurities  are  often  used  to  generate  pn  junctions  In  device 
preparation,  and  the  depth  of  this  junction  below  the  silicon  surface  is 
an  Important  factor  In  device  behavior.  Of  the  factors  mentioned  above, 
one,  the  presence  of  multiply  charged  Ions,  plays  a  significant  role  In 
affecting  junction  penetration,  x.,  and  Its  presence  can  be  easily 

J 

detected  from  measurements  of  junction  penetration.  Such  measurements 
were  used  to  characterize  the  experimental  Implantations.  This  x^  value 
was  compared  with  the  values  obtained  from  conventional  Implantations 
performed  with  the  same  energy  and  fluence. 

While  techniques  for  measuring  Xj  have  been  available  for  several 
decades,  they  have  never  been  satisfactorily  applied  to  the  very  shallow 
junctions  obtained  with  ion  Implantation.  A  recent  paper  by  C.  P.  Wu  and 
his  associates  at  RCA  Laboratories,  Princeton, "Techniques  for  Lapping 
and  Staining  Ion-Implanted  Layers,"  J.  Electrochem.  Soc.  126,  1982  (1979) 
describes  a  series  of  precautions  that  are  required  to  adapt  present  angle- 
lap  techniques  to  the  measurement  of  such  shallow  junctions.  We  have 
developed  a  new  technique  which  we  feel  gives  the  requisite  accuracy  with¬ 
out  the  Investment  of  excessive  time  and  care  that  prior  techniques  appear 
to  require.  This  technique,  which  we  Intend  to  publish  In  the  scientific 
literature  Is  described  under  the  section  on  Groove  and  Stain  Technique. 
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3.5  Implant  Damage 

Together  with  the  desired  Introduction  of  Impurity  atoms  Into  the 
silicon  lattice  there  also  occurs  undesired  lattice  damage.  The  thermal 
anneal  which  was  used  here  permits  the  Implanted  atoms  to  migrate  to  lattice 
positions  and  thus  become  electrically  active.  It  has  been  shown  (4,  5) 
that  a  considerably  more  drastic  anneal,  i.e.  higher  temperature  at  a 
longer  time  Is  required  to  remove  the  concomitant  lattice  damage. 

A  single  technique  that  has  been  developed  for  studying  implant  defects 
consists  of  exposing  the  ion  implanted  surface  to  a  thermal  oxidation  such 
as  1  hour  In  wet  oxygen  at  1100°C.  This  causes  the  implant  defects  to 
expand  to  stacking  faults  (SF)  which  are  several  orders  of  magnitude  larger 
and  can  thus  be  detected  by  optical  microscopy.  This  is  accomplished  by 
stripping  off  the  grown  Si 02  with  an  HF  rinse  and  subjecting  the  bare  sil¬ 
icon  surface  to  a  short  Sirtl  etch.  Figure  7  illustrates  the  stacking 

14 

fault  density  found  for  a  1  x  10  B  cm  fluence  implanted  by  a 
conventional  ion  implanter. 

It  is  our  intent  to  evaluate  the  implant  damage  developed  by  our 
experimental  technique  in  two  ways.  The  first  is  to  compare  the  SF  con- 
centrationsobtained  for  the  experimental  technique  with  that  obtained  for 
a  conventional  implantation  of  the  same  fluence  and  energy.  The  second  is 
to  determine  the  thermal  annealing  conditions  necessary  to  dissolve  the 
defects  generated  by  the  experimental  implantation  and  to  compare  it  to 
the  annealing  conditions  required  for  a  similar  conventional  implantation. 
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Ion  Implant  Boundary,  10  B  cm,  50  keV  (100)  wet 
Oxidized  at  1100°C  for  1 h ,  Slrtl  etched  30  sec. 
Magnl flcatlon  405X 


3.6  Groove  and  Stain  Technique 

The  measurement  of  shallow  junctions  characteristic  of  Ion  implantation 
raises  special  problems.  We  have  developed  a  technique  to  permit  accurate 
reproducible  measurement  of  such  junctions,  using  a  groove  and  stain 
technique  together  with  a  profllometer.  There  are  two  variants  to  this 
procedure,  one  a  fairly  time-consuming  one,  and  the  other  a  short 
convenient  one.  In  this  study  the  latter  was  used. 

Our  specimens  in  this  study  were  three  n-type  wafers  of  1-4  flcm  resistivity 
conventionally  implanted  with  1  x  10^,  1  x  10^4  and  1  x  10^  cm  ^  at 
75  keV,  and  a  10  flcm  n-type  wafer  experimentally  implanted  with  the  laser 
blow-off  technique  with  approximately  2  x  10^4  B  (natural  isotopic 
abundance)  at  50  keV.  All  wafers  were  subjected  to  a  two  step  anneal, 
lh  at  550°C,  and  10  minutes  at  1000°c. 

The  process  began  by  grooving  the  wafer  surfaces.  For  this  procedure, 
lower  than  conventional  pressures  were  used  so  that  the  groove  depth 
exceeded  the  junction  depth  only  in  part  of  the  groove  length.  In  other 
grooved  areas  the  groove  was  shallower  than  the  junction.  After  staining, 
the  llgher  n-type  region  appeared  only  where  the  groove  extended  below  the 
pn  junction.  Figure  8  illustrates  a  stained  groove  length,  part  of  which 
was  deeper  than  the  junction  and  part  of  which  was  shallower.  As  Section 
A-A,  the  groove  depth  was  exactly  equal  to  the  junction  depth.  Measuring 
the  groove  depth  at  this  point  permits  one  to  obtain  the  junction  depth. 

Figure  9  exhibits  the  profiles  measured  at  these  points  giving  values  for 
the  junction  depths. 

The  following  results  were  obtained: 


FLUENCE 

ENERGY 

IMPLANT  TYPE 

Xj,  A 

TT^B11 

75  keV 

Convent. 

4800 

1  x  1014  B11 

75  keV 

Convent. 

5200 

1  x  1015  B11 

75  keV 

Convent. 

5500 

*  2  x  1014B 

50  keV 

Laser  Blowoff 

6000 
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FIGURE  8.  Microphotograph  of  grooved  and  stained  wafer.  Mag  100X 
Profile  measured  at  A-A 
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FIGURE  9  Profiles  for  four  implantations.  Horizontal  scale  50.1 

Vertical  scale  5000:1 


4.  ION  MICROPROBE  COMPARISONS 

4.1  Preparation  of  Standards  by  Conventional  Ion  Implantation 

In  order  to  make  Ion  microprobe  comparisons  of  silicon  Implanted  with 
boron  by  conventional  techniques  to  silicon  Implanted  with  the  laser 
blow-off  techniques,  three  wafers  were  Implanted  with  a  conventional 
machine. 

4.1.1  Equipment  Oecriptlon; 

The  wafers  were  Implanted  with  a  Varlan/Extrion  Ion  Implanter  model 
200-20A2F  with  a  30CN  end  station.  The  200-20A2F  Is  a  high  throughput, 
multi-purpose  production  Ion  Implanter  with  an  energy  range  from 
10  to  200  KeV.  It  utilizes  a  filament  type  Ion  source  to  achieve  beam 
currents  of  150  pA  for  boron  and  up  to  600  pA  for  phosphorous  and 
arsenic.  An  electrostatic  deflection  system  Is  used  to  scan  the  beam. 
Figure  10  shows  a  schematic  diagram  of  this  conventional  machine. 

4.1.2  Process : 

The  wafers  were  Implanted  with  B^+  at  the  fluences  listed  below: 


Dose 

Beam  Current 

Wafer  1: 

Wafer  2: 

Wafer  3: 

1  x  1013/cm2 

1  x  1014/cmZ 

1  x  1015/cm2 

3.0  pA 

26.0  pA 

125.0  pA 

The  Incident  energy  remained  constant  at  75  KeV. 


An  implant  angle  of  7°  was  used. 
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4.2  Preparation  of  Annealed  Ion  Implant  Standards 

In  order  to  measure  the  amount  of  boron  redistribution  due  to  annealing, 
and  to  compare  any  differences  for  the  case  of  laser  blow-off  Implants, 
silicon  wafers  (10  ohm-cm  phosphorous-doped  n-type)  Implanted  with  boron 
by  conventional  means  were  annealed  according  to  the  following 
procedures. 

4.2.1  Cleaning  procedure: 

The  wafers  are  loaded  In  a  quartz  carrier  and  cleaned  In  pyrex  beakers. 

SN  -  The  wafers  are  Immersed  In  3  successive  baths  of  sulfuric 
acid  at  150°C  for  5  minutes  minimum  In  each  bath. 

-  The  wafers  are  then  Immersed  In  a  nitric  acid  bath  at  55°C  for 
a  minimum  of  5  minutes. 

-  The  wafers  are  then  rinsed  In  flowing  deionized  water  for  a 
minimum  of  1  minute. 

0E  -  The  wafers  are  Immersed  In  7:1  buffered  oxide  etch  for 

exactly  3  seconds  followed  by  Immediately  Immersing  In  flowing 
deionized  water  for  at  least  1  minute. 

PA  -  The  wafers  are  Immersed  In  a  1:1  boiling  mixture  of  freshly 
prepared  H202  and  NH^OH  for  exactly  25  minutes. 

30' R  -  The  wafers  are  rinsed  In  flowing  deionized  water  for  exactly 

30  minutes. 

4.2.2  Annealing  Procedure: 

The  wafers  are  annealed  In  a  diffusion  furnace  tube  at  900°C  for 
1  hour  in  an  ambient  of  flowing  nitrogen. 
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4.3  Procedure  for  Profiling  Boron  Implantation  In  Silicon  Wafers 
using  the  Ion  Microprobe  (IMMA) 

4.3.1  Introduction 

This  section  outlines  the  procedure  for  determining  profiles  of  boron 
Implanted  In  silicon  wafers.  Only  qualified  personnel  Instructed  In  IMMA 
operation  in  conjunction  with  the  system  manuals  supplied  by  the 
manufacturer  make  measurements  using  this  procedure.  The  procedure 
describes  operator  functions  but  avoids  descriptions  of  Instrument 
functions. 

4.3.2  Samples 

4. 3. 2.1  Sample  Preparation.  The  samples  are  prepared  and  Identified 
by  requestor.  After  samples  are  prepared  they  are  stored  In  clean 
sample  holders  and  transported  to  the  IMA  facility. 

4. 3. 2. 2  Sample  Constraints.  Due  to  the  design  of  the  Instrument 
sample  chamber,  all  IMMA  samples  must  conform  to  the  following 
specifications: 

1)  The  difference  In  thickness  of  the  sample  on  a  given  sample 
holder  shall  not  exceed  +  0.009  Inches.  The  difference  between 
the  highest  point  on  the  sample  holder  and  the  lowest  point  that 
can  be  analyzed  Is  <  0.020  Inches. 

2)  The  surface  of  the  sample  must  be  conductive.  A  conductive  path 
must  exist  between  the  sample  surface  being  analyzed  and  the 
sample  holder. 

3)  The  sample  must  be  a  solid. 

4)  The  sample  cannot  exceed  1-1/4  inch  In  diameter. 
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4.3.3  Sample  Receipt 


The  samples  or  sample  Identification  are  logged  In  the  IMMA  sample 
log  book. 

4.3.4  IMMA  Sample  Holder  and  Fixture  Cleaning  Procedure 

1)  Ultrasonic  clean  in  acetone  to  remove  conductive  paint  on 
sample  fixture  surface. 

2)  Ultrasonic  clean  in  trichlorethylene. 

3)  Parts  are  thoroughly  dried  by  placing  in  a  clean  oven  at  ^200°C 
for  1-2  hours. 

4)  After  drying,  sample  fixture  is  stored  In  a  clean  dry  area  such 
as  a  desiccator  containing  silica  gel  or  seal  in  a  clean 
plastic  bag. 

4.3.5  Sample  Fixture  Loading 

The  sample  or  samples  are  attached  to  a  clean  sample  holder  and 
placed  In  a  cleaned  IMMA  sample  fixture.  The  highest  point  of  the 
sample  holder  Is  made  level  with  the  top  surface  of  the  fixture.  The 
sample  fixture  is  placed  in  one  of  the  four  positions  of  the  sample 
select  table.  A  total  of  three  fixtures  can  be  loaded  into  the  sample 
chamber.  The  fourth  position  on  the  sample  table  contains  alignment 
standards . 


4.3.6  Stage  Auto  Yent  and  Pump  Down 

During  either  stage  pump  down  or  venting  operations  it  is  unnecessary 
to  shut  down,  or  make  adjustments  to  any  other  parts  of  the  system. 
Complete  sequencing  of  all  vacuum  operations  Is  Initiated  in  the  auto 
vacuum  mode. 


4.3.6  (Continued) 

It  Is  necessary  to  lower  the  sample  elevator  to  allow  closing  of  the 
gate  valve.  The  average  pump  down  time  Is  30  minutes;  however.  If  the 
stage  Is  open  for  extended  periods,  or  samples  out  gas  heavily,  pumping 
time  extended. 

4.3.7  Operating  Conditions 

4. 3. 7.1  Beam  Alignment.  When  the  chamber  Is  pumped  to  4  x  10  torr, 
the  alignment  standard  Is  raised  to  the  analyze  position.  The  alignment 
Is  performed  at  the  beginning  of  a  series  of  samples  but,  no  less  than 
at  the  beginning  of  each  work  day. 


4. 3. 7. 2  System  Operation  Conditions 

1)  Cold  plate  on  If  required  -  LN2 

2)  Primary  beam  -  02+,  20  kV,  30  nA 

3)  Raster  size  *  100  x  100  microns 

4)  Raster  rate  -  100/sec 

5)  Lines/Raster  -  100 

6)  Function  -  Cr 

7)  Scope  1  mode  -  XLP 

8)  Scope  Input  -  secondary  Ion 

9)  Sample  bias  -  H.V.  position 

10)  Scaler/timer  -  100  seconds 

11)  Set  aux.  drive  -  0.7 

4.3.8  Measurements 

For  depth  profile  measurements  a  flow  diagram  is  given  In  Figure  11. 
The  following  Is  a  step  wise  operation  to  make  said  measurements. 

1)  Turn  Ion  beam  on  and  align  beam,  refer  to  manual  for  Instructions. 
Use  operating  conditions  described  in  4. 3. 7. 2. 


Receive  and 
Set-up 
Sample 


time  scale  Is  for  one  sample.  More  than  one  sample 
can  be  loaded  and  pumped  down  at  the  same  time. 


4.3.8  (Continued) 

2)  Tune  secondary  beam  (refer  to  manual) 

3)  Load  computer  program  DEPTH 

a.  Computer  input  from  operator 

1.  Sample  bias  adjust  -  no 

2.  DAC  adjust  -  yes 

3.  Operating  conditions  -  (operator  selected) 

4.  Data  file  identification  -  (operator  selected) 

5.  Masses  to  be  profiled 

a)  M/E  10  (B)  count  time  =  5  sec 

b)  M/E  11  (8)  count  time  =  5  sec 

c)  M/E  56  (S^)  count  time  -  3  sec 

6.  Masses  to  be  ploted  -  all 

NOTE:  M/E  28  (SI)  was  not  used  for  the  analysis  due  to  after 
glow  causing  high  background  at  M/E  10  and  11. 

7.  Total  profile  time  In  seconds  =  0  (max.  200  cycles). 

b.  Execute  program 

NOTE:  The  computer  drives  magnet  to  selected  mass  (lowest  mass 
first)  counts  for  preset  time.  Moves  to  next  mass,  etc., 
and  repeats  cycle  for  total  profile  time  while  continuously 
sputtering  into  the  sample.  Count  data  Is  plotted  and  stored 
for  future  reduction.  Upon  completion  of  the  profile  the 
operator  can  select  to  repeat  with  same  parameters  on  a 
new  area  of  the  sample  or  select  a  new  sample  for  analysis. 

4)  Data  Reduction 

a.  Load  computer  program  DPLOT. 

b.  Execute  program. 

NOTE:  Operator  selects  stored  data  Input.  The  computer  ratios 
M/E  10  and  11  to  M/E  56  and  plots  data. 
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ION  MICROPROBE  MASS  ANALYZER 


4.4  Results  of  Ion  Microprobe  Comparisons  of  Laser  Blow-Off  Ion 
Implants  with  Conventional  Implants 

Figure  13  shows  the  B  depth  profile  obtained  with  the  Ion  Mlcroprobe 

Mass  Analyzer  (IMMA)  for  the  B  Into  SI  standard  Implanted  with  the 

15  2 

conventional  machine  at  75  keV  and  a  fluence  of  10  B/cm  .  The  results 
are  as  expected.  Except  for  more  random  noise,  similar  results  were 
obtained  at  fluences  of  10  J  and  10  B/cm  ,  with  the  peaks  occurring  at 
the  same  sputter  times.  Intensities  differed  by  factors  of  xlO. 

Figure  14  shows  the  B  depth  profile  obtained  with  IMMA  for  B 

Implanted  Into  an  Identical  SI  wafer  at  50  keV  with  the  laser 

14  2 

blow-off  machine  at  a  fluence  of  2  x  10  B/cm  .  The  dashed  curve  shows 

the  expected  profile  for  a  conventional  Implant  at  50  keV.  The 

++  +++ 

difference  Is  due  to  the  presence  of  B  and  B  In  the  laser  blow-off 
plasma,  which  are  accelerated  to  100  keV  e.id  150  keV,  respectively,  for 
the  50  kV  potential  difference. 

Although  potentially  a  problem  to  some  applications.  It  Is  believed 
that  by  tailoring  the  acceleration  voltages,  this  high  energy  component 
will  not  present  serious  difficulties  In  most  practical  cases. 

See  next  section. 
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•  ORON  IN  SILICON 


SPUTTER  TIME  IN  SECONDS  (DEPTH) 


Figure  13.  Ion  microprobe  depth  profile  of  boron  Implanted  Into  silicon 
by  conventional  Ion  beam  Implanter  at  75  keV  at 
1  x  1015  B/cm2  fluence 


-93- 


INTENSITY  RATIOED  TO  SIUCON  (Sit) 


BORON  IN  SILICON 


Figure  14.  Ion  microprobe  depth  profile  of  boron  Implanted  Into  silicon 
with  laser  blow-off  Ion  Implanter  at  50  keV  and  2  x  1014  B/cm2 
fluence.  Dashed  curve  shows  profile  for  50  keV  conventional 
Implant.  The  difference  Is  due  to  plasma  components  of  B4"*-  and 
B+++  accelerated  to  energies  of  100  and  150  keV.  The  abscissas 
should  be  directly  comparable  to  the  two  preceding  figures. 
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4.5  Nature  of  Non  Mass-Analyzed  Ion  Implantation 

Results  of  IMMA  analysis  of  Ion  Implanted  Impurity  distributions  for 
conventional  and  laser  blow-off  Ion  Implantations  show  very  clearly  that 
for  laser  blow  off  ion  Implantation  a  significant  part  of  the  Implant 
current  Is  made  up  of  multiply  charged  Ions. 

An  approximation  of  the  IMMA  profile  suggests  that  the  Implanted  boron 

beam  contains  60%  singly  charged,  30%  doubly  charged  Ions  and  10%  triply 

charged  Ions.  Figure  15  Illustrates  the  B  distribution  expected  when 
14 

1  x  10  B  Is  Implanted  In  a  conventional  Implanter  Into  silicon  at  50  keV 

with  only  singly  charged  Ions.  Figure  16  Illustrates  the  B  distribution 

14 

expected  when  1  x  10  B  Is  Implanted  at  50  keV  with  60%  singly  charged, 

30%  double  charged  and  10%  triply  charged  Ions.  Note  that  In  the  latter 
case  there  Is  a  significantly  deeper  penetration. 

The  deeper  penetration  of  the  multiply  charged  Ions  does  not  pose  a 
serious  problem  since  the  laser  blow  off  Implanter  permits  the  use  of 
low  Implant  energies.  Thus  the  Implant  energy  can  be  reduced  to  obtain 
the  desired  junction  depth. 

Figures  15  and  16  also  exhibit  the  B  profiles  afer  an  anneal  of 

30  minutes  at  1000°C.  Although  It  appears  that  the  thermal  effect  Is 

greater  for  the  singly  charged  Implantation,  this  Is  not  necessarily 

true.  The  deeper  multiply  charged  Implantation  has  a  shallower  gradient. 

Using  a  lower  energy  for  the  multiply  charged  Implant  results  In  a 

shallower  penetration,  a  steeper  concentration  gradient  and  a  greater 

thermal  diffusion  as  compared  to  results  shown  In  Figure  16.  Figure  17 

14 

Illustrates  the  profiles  obtained  for  multiply  charged  10  B  Implantations 
at  15,  20  and  25  keV.  Figure  18  allows  one  to  compare  the  profile  of  a 
singly  charged  B  Implantation  at  50  keV  with  that  of  a  multiply  charged 
Implantation  of  the  same  fluence  at  20  keV.  We  note  that  both  the 
junction  depth  and  the  concentration  gradient  at  the  junction  are  quite 
similar,  leading  us  to  conclude  that  a  multiply  charged  Implantation  with 
a  correctly  chosen  energy  could  be  substituted  for  a  singly  charged 
conventional  Implantation. 
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FI6URE  15.  Depth  Profile  of  Singly-Charged  Ion  Implantation  as 

Implanted  and  after  Anneal 
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FIGURE  18.  Comparison  of  Slngly-and  Multi  ply- Charged  Ion  Implantation 

Depth  Profiles 
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4.6  RESULTS  OF  EVALUATION  OF  LASER  BLOW-OFF  ION  IMPLANTED  DEVICE 

Laser  blow-off  Ion  Implantation  was  utilized  for  base  production 
on  device  test  wafers.  Boron  was  Implanted  Into  silicon  at  appropriate 
energy  and  dose  so  as  to  replicate  the  conventionally-implanted  boron 
distribution.  The  devices  were  16-bit  multipliers. 

The  devices  were  tested  after  completion  of  fabrication.  It  was 
found  that,  although  none  of  the  devices  (containing  1  O' s  of  thousands 
of  transistors)  were  entirely  functional,  the  Individual  transistors 
showed  transistor  action.  Beta  curves  are  shown  on  Figure  19.  All 
Indications  point  to  a  boron  dose  too  low  by  15  -  20%  as  the  source 
of  the  difficulty.  Resources  did  not  allow  a  repeat  of  the  Implantation 
and  subsequent  wafer  processing. 

These  first  device  tests  Indicate  that  laser  blow-off  Ion 
Implantation  Is  a  viable  method  for  Implanting  semiconductor  material  to 
produce  transistors. 
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FIGURE  19.  Beta  Curves  for  Transistors  Produced  by  Laser  Blow-Off 
Ion  Implantation.  Curves  show  a  beta  of  125;  top  line 
is  10  A. 
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THERMAL  TRANSPORT  IN  THE  TARGET  MATERIAL  DURING  LASER  BLOW-OFF 
IMPLANTATION 


Summary 

The  time  evolution  of  temperature  profiles  was  calculated  for  pulsed 
laser  blow-off  ion  implantation.  It  was  shown  that  for  practical  values 
of  ion  energy  and  ion  current  density  the  temperatures  reached  in  the 
implanted  regions  are  too  low  to  produce  efficient  self  annealing  by 
thermal  effects  alone  of  the  damage  created  during  the  implantation. 

Introduction 

The  laser  blow-off  implantation  technique  utilizes  pulsed  laser 
irradiation  of  a  solid  surface  of  pure  dopant  material  to  produce  a  high 
temperature  plasma  of  dopant  ions.  The  plasma  expands  into  vacuum, 
reducing  the  charge  density  to  the  point  where  an  externally  applied 
electrical  field  can  be  used  to  extract  ions  from  the  plasma.  In 
principle  these  ions  can  be  accelerated  to  an  arbitrarily  high  energy 
and  can  be  used  to  implant  semiconductor  materials.  The  method  offers 
several  advantages  over  the  conventional  ion-implantation  technique.  The 
laser  blow-off  method  is  relatively  simple,  very  fast  (the  pulse  duration 
is  several  microseconds),  permits  large  aerial  coverage,  and  is 
amenable  to  a  very  wide  range  of  dopant  species.  It  has  also  been  suggested 
that  because  of  the  pulsed  nature  of  the  ion  beam  current,  the 
temperatures  reached  in  the  implanted  region  during  the  laser  blow-off 
implantation  could  be  very  high.  This  is  in  contrast  to  the  conventional 
implantation  technique,  where  the  host  crystal  remains  at  low  temperature 
during  the  implantation.  A  deleterious  side  effect  of  ion  implantation 
is  damage  to  the  host  crystal  caused  by  the  high  energy  ions  necessary 
to  penetrate  to  the  desired  depths.  This  damage  must  be  removed.  This 
is  achieved  through  further  processing  of  the  implanted  samples,  such 
as  a  high- temperature  thermal  annealing,  or,  more  recently,  laser- 
induced  annealing.  The  possibility  of  reaching  high  temperatures  In  the 
Implanted  region  during  the  laser  blow-off  Implantation  might  alleviate 
this  requirement  of  further  processing.  If  sufficiently  high,  the 
temperatures  produced  through  the  use  of  high  intensity,  short  duration 
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ion  beam  pulses  might  permit  self-annealing  of  the  implant  damage  during 
the  implantation. 


The  present  work  addresses  the  possibility  of  reaching  high 

temperatures  in  the  implanted  region  during  the  laser  blow-off.  It  presents 

the  results  of  the  thermal  calculations  inside  the  implanted  material.  In 

the  model  used  in  this  work,  it  is  assumed  that  the  heat  is  generated  along 

the  path  of  the  penetrating  ions  at  the  same  rate  at  which  the  ions  lose 

their  energy.  In  the  passage  of  heavy  charged  particles  through  matter, 

the  ion  energy  is  lost  through  ionization,  excitation  (without  ionization) 

of  atoms  in  the  absorber,  capture  and  loss  of  electrons,  and  elastic 

nuclear  scattering.  Thus,  the  assumption  that  the  heat  is  generated  at 

the  same  rate  as  the  ions  lose  their  energy  implies  that  the  transfer  of 

the  energy  through  these  processes  to  the  lattice  proceeds  at  a  rate 

which  is  faster  than  the  inverse  of  the  ion  beam  current  pulse  length. 

This  assumption  is  well  satisfied,  the  typical  pulse  length  and  thermal 

-5  -8 

equilibrization  times  being  on  the  order  of  10  sec  and  10  sec, 
respectively. 

Because  of  the  large  aerial  cross-section  of  the  beam,  the  thermal 
transport  will  be  assumed  one-dimensional  in  this  work.  Further  assumptions 
will  involve  approximating  the  current  beam  pulse  shape  by  a  square 
shape  pulse,  and  approximating  the  values  of  the  material  physical 
parameters  by  constants.  These  assumptions  and  mathematical  methods  are 
discussed  in  more  detail  in  the  Model  Section.  The  results  are 
discussed  in  the  Discussion  Section. 
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Model 


The  temperature  distribution  In  the  semiconductor  In  the  one-dimensional 
approximation  is  obtained  by  solving  the  following  thermal  equation, 

3T  3 / < 3 T  \  $  /  -I  \ 

3t  '  3*  c  ax  "  c„  ’  v  ' 

\  P  /  P 

where  cp  and  k  are  the  "specific"  heat  and  thermal  conductivity, 
respectively,  and  where  the  source  function  S(S  =  S(x,t))  Is  equal  to  the 
absorbed  power  per  unit  volume. 

S(x,t)  =  a(J/e),  (2) 

where  a(a  =  a(x,t))  Is  the  energy  loss  in  the  bombarded  material  per  unit 
length  per  ion,  J(J  =  J(t) )  Is  the  incident  current  density,  and  e  is  the 
charge  of  the  electron.  In  (2),  it  has  implicitly  been  assumed  that  the 
majority  of  the  ions  are  singly  ionized.  Neglecting  loss  of  heat 
radiated  from  the  front  surface.  Equation  (1)  is  then  solved  subject  to 
the  following  initial  and  boundary  conditions: 


T(x,t  =  0)  =  Tq 

(3) 

02 

-1 

Q>  X 

X 

II 

o 

• 

rt 

II 

O 

(4) 

T(x  *  £,  t)  =  Tq 

(5) 

where  t  is  the  sample  thickness  and  Tq  Is  temperature  of  the  heat  sink  at 
the  back  side  of  the  sample. 

In  Equation  (1),  the  material  constants  cp  and  k  are.  In  general, 

functions  of  temperature.  However,  in  the  temperature  region  of  interest, 

cp  Is  practically  constant,  and  k  varies  from  1.2  H/cm®K  at  300°K  to 

.3  W/cm°K  at  1000°K  for  pure  silicon,  and  is  lower  for  doped  silicon.^ 

Hence,  In  what  follows,  the  temperature  dependence  of  c  and  *  will  be 
.  1  3  P 

neglected  (from  our  past  experience  *  the  error  Introduced  by 

approximating  k  by  the  "average"  value  of  *  •  .6  W/cm°K  Is  on  the  order 

of  some  5 X  In  the  value  of  the  calculated  temperature  profiles). 
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Making  the  approximation  that  cp  and  <  are  constant  permits  one  to  solve 
Equation  (1)  anal tyl cal ly  by  using  Green's  function  techniques. 

It  is  easy  to  verify  that  for  the  boundary  conditions  Equations  (3)- 
(5),  the  Green's  function  for  Equation  (1)  is 


G(x,x' ,  t-t') 


1 

'/4rr0(t-t') 


£ 


H)n 


(6) 


(~„-Ax-x'  -  2nt)Z\.  ^(x+x'  -  2n*)2  \  ) 

|exp^  4b{t^F)"  J+  exp  (  WtV)  /  )  f  * 

where  the  diffusion  constant  D  =  k/c .  The  temperature  distribution  Is 
the  convolution  of  the  Green's  function  and  the  source  term  S, 

«  t 

AT(x,t)  =  T(x,t)  -  T  a  f  dx'  f  G(x,x* ,  t-t')  S(x',t')/c  .  (7) 

°  */0  p 


Assuming  the  slab  thickness  I  »  (xR,  xth),  where  xR  Is  the  "range" 
of  the  implanted  Ions,  and  x^  is  a  characteristic  thermal  diffusion 
length,  xth  *  V^D  tp,  where  tp  Is  the  pulse  length,  the  semiconductor  can 
be  replaced  by  a  semi-infinite  medium,  with  the  result  that  the  sumnation 
In  Equation  (6)  reduces  to  only  one  term,  namely  the  term  n  =  0. 
Substituting  the  resulting  equation  into  Equation  (7)  then  yields 


AT(x,t) 


S(x\f) 

■«  4irD(t-t'  ) 


(8) 
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In  Equation  (8),  the  source  term  S  Is  generally  time  and  coordinate- 
dependent.  The  coordinate  dependence  arises  because  of  the  dependence 
of  the  Ion  energy  loss  a  on  depth  below  the  material  surface.  The  time 
dependence  arises  primarily  because  of  the  pulse  shape,  and,  to  a  lesser 
extent,  because  of  a  finite  transit  time  In  the  material  before  the  Ions 


stop  and  transfer  their  energy  to  the  crystalline  lattice.  However,  for 

typical  parameters  of  Interest,  the  stopping  time  has  been  estimated  to 

-12  -11 

be  on  the  order  of  10  to  10  sec.  This  Is  much  shorter  than  the 
typical  pulse  length  t  of  several  microseconds.^  Thus,  the  time 
dependence  of  a  can  be  neglected.  The  shape  of  the  current  pulse  has  been 
described  In  Reference  (1).  The  shape  somewhat  resembles  the  "Gaussian" 


pulse  shape. 


0  <  t  <  t 

m 


otherwise. 


(9) 


where  t  Is  the  time  at  which  the  pulse  reaches  Its  maximum  value  J  . 

m  8  o 

Again,  from  our  past  experience,  we  have  found  that  by  approximating 

the  Gaussian  pulse  shape  by  a  square  shape  pulse. 


0  <  t  <  t 


P 


Otherwl se. 


(10) 


where  the  pulse  length  t  *  2  t^,  Introduces  *  very  small  error  (some  5%) 
Into  our  final  results.  Furthermore,  using  the  square  shape  pulse 
permits  us  to  Integrate  Equation  (8)  analytically.  Therefore,  In  what 
follows,  we  shall  approximate  the  current  pulse  by  a  square  shape  pulse. 
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Finally,  the  energy  loss  per  unit  length  will  be  approximated  by  a 
constant. 


where  E  Is  the  ion  energy  (In  electron  volts),  and  xR  Is  the  projected 
range.  A  better  approximation  to  a  can  be  obtained  through  the  use  of 
the  specific  Ionization  curves  for  passage  of  heavy  charged  particles 
through  matter,  discussed  In  Reference  4.  From  Reference  4,  a 
typical  specific  Ionization  curve  for  a  heavy  charged  particle  start*:  at 
the  surface  at  about  1/2  of  Its  maximum  value,  then  slowly  rises  with 
depth  beneath  the  surface,  reaches  Its  maximum  value  at  some  distance  x 

m 

and  then  falls  abruptly  to  zero.  A  typical  (normalized)  plot  of  the 
specific  Ionization  I(x)  Is  shown  in  Figure  1.  However,  as  has  already 
been  discussed  in  the  Introduction  Section,  other  processes  such  as  host 
atom  excitation,  capture  and  loss  of  electrons,  and  nuclear  scattering, 
also  contribute  to  the  Ion  energy  loss.  Thus,  the  coordinate 
dependence  of  the  energy  loss  o  per  unit  length  will  not  necessarily 
correspond  exactly  to  the  specific  Ionization  curve.  In  the  spirit  of 
the  other  approximations  In  this  calculation,  we  shall  approximate  a  by 
a  constant,  chosen  such  that  the  product  of  a  and  projected  range  xR 
Is  equal  to  the  Ion  energy  eE  (see  Equation  (11)). 


Making  the  above  discussed  approximations,  the  source  term  S(x,t) 
reduces  to  a  constant  for  times  0  <  t  <  tp  and  coordinates  0  <  x  <  xR, 
and  Is  equal  to  zero  otherwise.  Substituting  S(x,t)  Into  Equation  (8) 
and  making  the  change  of  variables  t  -  t*  «  t"  and  x"  *  (x  +  x' )/ V^Dt" 
then  yields 
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A  typical  average  specific  Ionization  along  a  collimated 
beam  of  heavy  particles  In  matter 


AT(x,t)  = 


r  hrR  T,  f  /  dx 

P  R  (...)  Viot7  Q  \  / 


0 

t 


(12) 


cp  ^<i,er,W)  ' 


In  Equation  (12),  the  time  tQ  Is  tQ  =  max  (t-tp,0) 


,t  -  t. 


^‘p 

tit 


03) 


Equation  (12)  can  be  evaluated  by  using  changes  of  variables  and 
partial  fractions.  First,  using  the  change  of  variables  1/Vt"  * 

(xR  +  x)/\Z4Dt'  sign  (xR  +  x),  where  sign  (x)  =  +  1  for  x  £  0,  Equation 
(12)  tranforms  Into 

t(±) 

4T'  5  >  „)  =  iD  X.  s,9"<Vx)  <XSXR >2/  dt'  erf(vjf"/ 


(14) 


where  the  lower  and  upper  Integral  boundaries  are 


(!) 


1 


Application  of  a  partial  fraction  to  Equation  (14)  yields: 


AT‘  *  M  J}ign  (xR±x)(xR-x)< 
(i) 


r  *«*)  ‘*(i) 
l*'  erf^Jt  (i)v^  (u)  fw 

ll  X1 


-1/fl 


The  second  bracketed  term  on  the  right  hand  side  of  Equation  (16)  is 

1  /2 

evaluated  by  using  the  change  of  variables  t'  a  1/x  and  a  partial 
fraction.  This  yields: 

to  x-. 


/dt  -1/t  f  2dx  .-x 

T?'  IT 

t,  1 


e  =  2 


~x 


*/ 


dx  e 


-x 


(17) 


-x 


=  2 


+  yfr  erf(x) 


'1 


Substituting  Equation  (17)  Into  Equation  (16)  finally  yields: 
aT'  =  sf  S  s,9"  (xR±x)(x±xR)2  j  (2»f>  z/TTe 


-1/t' 


t2(0 


(16) 


(O 


(18) 


where  the  lower  and  upper  boundaries  t^  +  ^  and  t2^  are  given  by 
Equation  (15). 
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Results  and  Discussion 


Equation  (18)  was  evaluated  numerically  for  several  values  of  the 

current  pulse  length  tp,  and  several  values  of  the  Ion  energy  E.  The 

Implanted  ions  were  boron.  The  projected  range  x^,  corresponding  to  a 

particular  energy  E,  Is  tabulated  in  Table  1.  The  energy  ranged  from 

20  keV  to  180  keV,  and  the  corresponding  projected  range  was  from 

.0662  ym  to  .4872  yin.  The  pulse  length  ranged  from  t  =  10  ysec  to 

t  =  40  ysec.  The  substrate  material  was  silicon.  The  material 
P  3 

parameters  were  cp  *  1 .6  J/cm  °K,  and  *  »  .6  W/cm°K. 

Samples  of  the  numerical  results  are  shown  In  Figure  2-6.  The 

temperature  profiles  are  shown  as  a  function  of  the  normalized  distance 

x/xR  away  from  the  semiconductor  front  surface,  with  time  (in  fractions 

and  multiples  of  the  pulse  length  t  )  as  a  parameter.  The  temperature 

P  o 

shown  is  a  temperature  rise  per  unit  current  density  JQ(A/cm  ). 

Inspecting  the  calculated  profiles,  the  main  results  of  this 
calculation  can  be  sunmarlzed  as  follows: 

a.  The  temperature  Is  linearly  proportional  to  current. 

b.  The  maximum  temperature  Is  reached  at  the  end  of  the  ion  beam 
pulse,  at  time  t  =  tp. 

c.  The  temperature  profiles  are  thermal  diffusion-limited. 

The  conclusions  are  physically  easily  understood.  The  profiles  are 

thermal  diffusion-dominated  because,  for  the  beam  and  material  parameters 

used,  the  projected  range  Is  much  smaller  than  the  thermal  diffusion 

length,  x^  «  x^,  where  a  typical  value  of  x^  Is  x^  =  V  4r  tp/cp  * 

8  x  10‘3  cm  fOr  t_  *  40  ysec.  From  this  follows  that,  for  example, 

P  o 

for  tp  *  40  nsec,  current  density  1  amp/cm  and  the  ion  energy  E  *  200  keV, 
the  effective  heated  depth  Is  •  8  x  10"3  cm,  the  associated  rise  In  the 
silicon  temperature  Is 
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Table  1. 


LSS  Rang*  Statistics 


for 


Boron  Implanted  In  Silicon 
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SILICON 


-T.)/J 


•  ORON  INTO  SILICON 
TEMPERATURE  NSE  Kl  UNIT  CURRENT  DENSITY  J 
PULSE  LENGTH  TP  -  1.0  E  -  05  mc 
ENERGY  E  -  2.0  E  *  04  #V 
RANGE  R  *  6.6200  E  -  06  cm 


0  230  500  790  1000  1290  1900  1790 


X/R 


Figure  2.  Teuperature  profiles  for  pulsed  laser  blouoff  per  unit  current 
density  J.  Pulse  length  t  ■  10-5  sec.  Ion  energy  e  •  20  keV, 
projected  range  R  ■  6.62  x  1O-6  chi. 
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BORON  INTO  SILICON 
TEMPERATURE  RISE  PER  UNIT  CURRENT  DENSITY  J 
PULSE  LENGTH  TP  =  1.0  E  -  05  Me 


ENERGY  E  =  1.0  E-f  05  «V 
RANGE  R  =  2.9940  E  -  05  cm 


0  100  200  300  400 

X/R 


Figure  3.  Temperature  profiles  for  pulsed  laser  blouoff  per  unit  current 
density  J.  Pulse  length  tp  ■  10‘5  sec.  Ion  energy  E  ■  100  keV, 
projected  range  R  ■  2.994  x  10"5  cm. 
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BORON  INTO  SILICON 
TEMPERATURE  RISE  PER  UNIT  CURRENT  DENSITY  J 
PULSE  LENGTH  TP  =  2.0  E  -  05  mc 
ENERGY  E=  1.0  E+  05  *V 
RANGE  R  *  2.9940  E  -  05  cm 


Figure  4.  Temperature  profiles  for  pulsed  laser  blowoff  per  unit  current 

density  J.  Pulse  length  tp  ■  2  x  10"5  sec.  Ion  energy  E  ■  100  keV, 
projected  range  R  «  2.994  x  IQ*5  cm 
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(T  -  To)  /  J 


Figure  5 


BORON  INTO  SILICON 
TEMPERATURE  RISE  PER  UNIT  CURRENT  DENSITY  J 
PULSE  LENGTH  TP  =  4.0  E  -  05  mc 
ENERGY  E  =  1.0  E*  05  *V 
RANGE  R  =  2.9940  E  -  05  cm 


Temperature  profiles  for  pulsed  laser  blowoff  per  unit  current 
density  J.  Pulse  length  t  ■  4  x  10"®  sec.  Ion  energy  E  ■  100  keV, 
projected  range  R  ■  2.994  x  10’^  cm. 
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BORON  INTO  SILICON 
TEMPERATURE  RISE  PER  UNIT  CURRENT  DENSITY  J 
PULSE  LENGTH  TP  »  4.0  E  -  05  mc 


ENERGY  E  =  1.8E+  05  «V 
RANGE  R  =  4.8720  E- 05  cm 


0  >00  200  300  400  500 

X/R 


Figure  6.  Temperature  profiles  for  pulsed  laser  blowoff  per  unit  current 

density  J.  Pulse  length  t  ■  4  x  10"5  sec.  Ion  energy  E  ■  180  keV, 
projected  range  R  «  4.872  x  10"^  cm. 
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EJ  t 

AT  «  — £  —2—  »  700°K, 

e  Cp  Xth 


(19) 


and  the  time  it  takes  for  the  temperature  to  drop  to  about  1/2  of  its 
maximum  surface  value  is  t^  =  4tR. 


Inspection  of  Figure  2-6  confirms  the  essential  correctness  of  these 

elementary  estimates.  In  Figure  2,  the  maximum  surface  temperature 

per  unit  current  is  AT  =  80°K  for  E  =  20  keV  and  t  =10  nsec.  The 

m  P  -3 

temperature  distribution  halfwidth  is  =  500  xD  =  3  x  10  cm,  which  is  in 

K  ^ 

very  good  agreement  with  x  .  =  Y^/c  tT  =  3  x  10“  cm.  Figure  3  shows 

Cil  r  r 

the  temperature  increase  with  the  ion  energy  E.  Because  of  the  limit 

xth  >:>  XR'  temPerature  increases  with  E  approximately  linearly. 

Figures  4  and  5  show  the  effect  of  Increasing  the  pulse  length  tp.  For 

the  same  energy,  both  the  depth  of  the  heated  volume  and  the  surface 

temperature  increase  approximately  as  the  square  root  of  the  pulse  length 

tp.  Figure  6  shows  the  temperature  profiles  for  t  *  40  nsec  and 

E  =  180  keV.  These  values  represent  probably  the  maximum  achievable 

2 

pulse  length  and  ion  energy,  respectively.  For  JQ  =  1  A/cm  ,  the 

maximum  surface  temperature  is  AT  =  1,400°K,  and  the  surface  remains 

m 

at  this  elevated  temperature  during  a  time  Interval  of,  say,  one-tenth 

of  t  . 

P 


A  practical  annealing  of  damage  In  ion-implanted  semiconductors 
occurs  either  because  of  a  rapid  melting  and  regrowth  of  the  surface 
implanted  region,  or  because  of  a  slow  solid  state  regrowth  of  the 
crystalline  lattice.  For  example.  In  pulsed  laser-induced  heating,  the 
typical  laser  pulse  lengths  are  t  *  10  -  100  nsec,  laser  energy 
Q  2  1  J,  and  the  optical  absorption  depth  xab$  *  xth  *  .5  pm.  Performing 
elementary  estimates  similar  to  those  above,  and  assuming  reflection 
losses  of  50%,  the  temperature  rise  at  the  surface  is  AT  *  Q/e  x 

in  p  ADS 

5,000°K,  which  Indicates  surface  melting.  This  has  been  Indeed 
confirmed  by  more  detailed  calculations  and  experimental  observations. 
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The  temperature  profiles  as  calculated  here  for  the  Ion  bombardment  of 

2 

silicon  were  calculated  for  J  =  1  A/cm  .  In  practice,  this  current 

o  2 

value  Is  difficult  to  achieve,  a  more  probable  value  being  *  .1  A/cm  . 

Clearly,  the  temperatures  obtained  during  the  ion  bombardment  are 

insufficient  to  produce  surface  melting. 

The  solid-state  regrowth  annealing  occurs  during  cw  laser  processing. 

5  2 

There,  typical  laser  parameters  are  laser  power  P  -  2  x  10  W/cm  , 
scanning  velocity  v  =  2.76  cm  s"\  and  the  beam  cross-section  diameter 
d  =  38  pm.  Under  these  conditions,  the  heating  Is  thermal  diffusion- 
limited.  Using  the  same  elementary  analyses  as  above,  the  maximum 
surface  temperature  is  given  by  ATm  =  Pd/vcp  xth  =  £  Vd/vcp.  Assuming 
50%  reflection  losses,  this  gives  AT  =  1,500°K.  This  AT  lasts  during 
an  effective  heating  time  Interval  of  t^gat  s  d/v  *  15  x  10"H  sec,  and 
although  the  amount  of  the  absorbed  energy  is  insufficient  to  produce 
surface  melting,  the  Induced  temperature  is  very  close  to  the  material 
melting  temperature.  Comparing  these  estimates  with  those  for  Ion  beam 
bombardment.  It  Is  clear  that  heating  induced  during  the  Ion 
bombardment  is  too  low  to  produce  efficient  annealing  during  the 
bombardment. 
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